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ABSTRACT 


Significant advantages of Variable Polarity Plasma Arc (VPPA) 
welding process include faster welding, fewer repairs, less joint 
preparation, reduced weldment distortion, and absence of porosity. 
In this report, a mathematical model is presented to analyze the 
VPPA welding process. Results of the mathematical model have been 
compared with the experimental observation accomplished by the GDI 
team. 
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I . INTRODUCTION 


The plasma arc is a concentrated energy source commonly 
used in welding and cutting processes. It is composed of a 
partially ionized gas stream produced by forcing an inert gas 
to flow through an electric arc and emerge from a constricting 
nozzle. With its high energy density and velocity, the 
plasma arc, when impinging on the workpiece, can create a hole 
in the molten pool generated by its heat and can penetrate 
through the workpiece. Depending on the operating 
parameters employed, this hole may either become self- 
k ea ^ n< =T ot remain open as the arc transverses along the 
workpiece. A "keyhole" welding process occurs in the self- 
healing case, in which the molten metal in front of the arc 
flows around the arc column and resolidifies behind the arc. 
On the other hand, a cutting operation is achieved if the hole 
remains open (O'Brien, 1968 ) . 

The plasma welding process is essentially an extension 
of gas tungsten-arc welding (GTAW) ; however, constraining 
the arc to flow through a nozzle produces much higher energy 
density in the arc and much higher gas velocity and momentum. 
The process produces deep penetration welds by forming a 
keyhole in the workpiece by the pressure of the gas flow. The 

melted in front of the advancing keyhole flows around to 
the rear where it solidifies to form the weld bead. This is 
distinct from the essentially surface melting produced by the 
GTAW process (which cannot normally penetrate to a depth 

equal to the width of a weld pool) . 
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In 1955 Linde Air Products introduced a plasma arc torch 

for metal cutting applications, and by 1965 Linde had 

developed an automatic plasma arc welding facility for 

Westinghouse Electric Corporation to fabricate a 120 in. 

diameter, 3/8 in. thick D6AC steel rocket motor case for the 

Titan III-C booster assembly (Miller and Filipski, 1966; 

Privoznik and Miller, 1966) . The plasma arc welding process 

was reported to halve the welding time required. Then in 

1965 Thermal Dynamics Corporation reported the use of direct 

current reversed polarity (workpiece negative) plasma arc to 

join 1/4 in. thick aluminum plate (Cooper et al. , 1965) . At 

the end of the 1960 , s, Van Cleave at the Boeing Company began 

his efforts to combine the plasma arc process with a variable 

polarity feature in which the electrode polarity was 

periodically reversed. Alternating electrode potential for 

aluminum welding had been investigated as early as 1947 

(Herbst, 1947). Difficulties with welding power supplies in 

this application were evident early-on. In addition, when 

variable polarity plasma arc (VPPA) welding was used in the 

U.S. Army Roland Missile Production Program, development 

problems such as arc pressure pulsation were noted (Nunes et 

1981; 1983). As a result of Van Cleave's promising work 

at Boeing, a VPPA research and development project was 

initiated in 1978 at the NASA Marshall Space Flight Center, to 

determine the potential for replacing the GTAW system used in 

the fabrication of the Space Shuttle External Tank (ET) . 

After the original used in the fabrication and process were 
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improved, VPPA welding finally exceeded the expectations 
(Nunes et al . , 1984a; 1984b). 

The Space Shuttle ET (diameter 28 . 6 ft., length 154 ft.) 
is the largest known "drop tank, " carrying 140,000 gallons of 
liquid oxygen and 380,000 gallons of liquid hydrogen. From 
the outset of ET production, the conventional GTAW system 
operating in the direct current electrode negative mode was 
used. With the GTAW system, in spite of the requisite 
careful joint preparation, weld porosity and defects had been 
ever-present problems since the fabrication of the Saturn 
lunar rocket first stage. Therefore, the decision was made 
at the outset to use 100% radiographic inspection. Weld 
porosity and defects had to be ground out and repaired 
systematically . 

The VPPA process significantly reduced porosity and 
defects with less stringent joint preparation and, in 
addition, reduced weldment distortion, and speeded up the 
welding process. 

In this study, the power distribution for an argon 
plasma gas has been analyzed along its course through the VPPA 
welding system. The study includes the following sections: 

(1) Electric potential distribution of straight and 
reverse polarities, and potential intensities for 
the VPPA welding system; 

(2) VPPA welding system electric power input at the 
electrode, within orifice, in the free plasma jet 

column and in the keyhole of workpiece; 
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(3) Power loss at electrode, within the gap between 
electrode and within orifice, standoff column and 
the workpiece; 

(4) Determination of plasma arc enthalpy; 

(5) Determination of crown and root bead widths; 

(6) Determination of crown and root heights; 

(7) Determination of leading edge angle of keyhole. 

(8) Comparison between molding and experimental 
measurements 
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II. VARIABLE POLARITY PLASMA ARC WELDING SYSTEM 
ELECTRIC POTENTIAL CALCULATION 
The significant current carriers in a welding arc are 
electrons and positive ions . Electrons carry the bulk of the 
current, moving rapidly from negative cathode to positive 
anode. The positive ions drift more slowly through the 
interelectrode space. In a fixed polarity arrangement, the 
differential drift rate results in an asymmetrical heating of 
the welding arc's ends. The cathode receives less heat and 
the anode more heat . The straight polarity mode of operation 
entails a negative electrode (cathode) and a positive 
workpiece (anode) . Where the primary object of the weld 
process is to deliver the maximum heat to the workpiece with 

minimum deterioration of the electrode, straight polarity is 
used. 

Reverse polarity, which entails positive electrode and 

negative workpiece, has the advantage that the workpiece is 

subjected to a cleaning process (cathodic cleaning) by the 

impingement of positive ions on the workpiece surface. 

Recent work suggests that heavy positive ion weight is 

irrelevant to the cleaning process which seems to be more like 

capacitor breakdown than hammer impact. In the case of 

plasma arc welding, reverse polarity action appears to 

condition the surface of the aluminum alloy so that the molten 

metal flows easily and controllably under the arc. It is 

conjectured that this fluid control is accomplished through 

breaking up of surface oxide films. Cutoff of reverse 
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polarity during VPPA welding transforms a weld metal flow 
which closes smoothly and soundly behind the keyhold to an 
irregular, intermittent, globular flow. This process 
leaves a rough, lumpy bead sunken below the parent metal 
surface and protruding jaggedly from the root of the weld. 
However, continuous reverse polarity is not necessary to 
provide adequate cathodic cleaning action. 

The variable polarity square— wave with unequal straight 
and reverse polarity time intervals offered a combination of 
the high heating capability of straight polarity with the 
cleaning feature of reverse polarity. Adequate cleaning is 
obtained by incorporating a relatively short (1/10 to 1/5 the 
duration of the straight polarity current) pulse of reverse 
polarity current into the welding current waveform. 

Figure 1 show the VPPA electric power system. During 
the straight polarity mode of operation, it employs a 
negative electrode (cathode) and a positive workpiece 
(anode) ; during the reverse polarity mode of operation it 
employs a positive electrode and a negative workpiece. 
Figure 1 at left illustrates the straight polarity mode; at 
right, reverse polarity mode. 

Calculation of electric potential can be obtained from 
the following formulations: 

(A) For Straight Polarity 

V E = V F - AV & - * v 
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where 0 E = Electrode Work Function 
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4 V a — Anode Drop 
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(B) For Reverse Polarity 
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where A, B = constant varying with workpiece 
materials 

t + = Straight Polarity Time Period 
t_ = Reverse Polarity Time Period 
AV n - = AV A . F . = V A . - V F . = V A . 

AV P = • D * = V D 1 “ V A‘ 

(C) For Potential Intensities 

Electric potential intensities of the plasma fluid in 

the plasma column can be calculated from the electric 

conductivity, <r E , of plasma fluid at the location of 

interest. The value of electric conductivity at the 

location of interest is determined by the local plasma jet 

pressure and temperature. By using the data provided by 

Devoto (1973), Lancaster (1986) has calculated the 

electrical conductivity of argon plasma gas as it varies 

with temperature at one atmosphere pressure. Figure 2 

shows the variation of argon gas electrical conductivity 

with temperature at one atmospheric pressure. In this model 

8 



calculation, we have divided the lengths of plasma column 
along the axial direction into several segments with &L{ as 
the length of segment i. Joule heating, and radiative and 
convective heat transfer losses at the segment i are 
calculated with the assumed value of electric potential 
intensity of that segment, (aV/aL) j ' , and the given inlet 
values of power and enthalpy of plasma fluid at that 
segment. Outlet values of power and enthalpy of plasma 
fluid at that segment are then calculated based on these 
parameters. The value of electric conductivity <r E j at 
segment i is thus determined from the average fluid 
temperature at that segment. The value of average 
temperature is obtained from the corresponding average 
value of plasma fluid enthalpy at that segment which is 
calculated from the average value of inlet and outlet 
plasma fluid enthalpys at that segment. Electrical 
resistance, Rj ; electric potential difference, aVj ; and 
electric field intensity, (aV/aL)j, at segment i can be 
determined from the computed value of electric 
conductivity <r E/ j at that segment. Convergent solution is 
obtained through the trial and error compuations between 
the computed value of (aV/aL) £ and the assumed value of 
(aV/aLJj' at segment i. From these calculations, it is 
clearly shown that the values of electric potential of the 
plasma fluid and heat transfer computations are mutually 
coupled together. Furthermore, the values of plasma 
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electrical conductivity fluctuate with local variations of 
temperature. The physical parameters affecting the 
plasma temperature, such as plasma and shielding gas flow 
rates, radiative and convective heat transfer losses, 
values of main and pilot electric currents, standoff 
distance, etc., all influence the values of electric 
potential and electric field intensity in the plasma fluid. 

The effect of shielding gas flow rate occurs through 
shielding gas mixing into the plasma standoff column L 3 and 
vice versa. The continuous mixing between plasma gas and 
shielding gas induced by the plasma jet stream in plasma 
column results in local variations of species constituent 
concentration and temperature along the axial direction. 
Thus, the effect of jet mixing in column streams between 
plasma and shielding gases must be taken into consideration 
in addition to heat transfer losses in computing electric 
potential along the plasma jet. Figure 3 shows the 
variation of electrical resistance with temperature for 
various argon-helium mixtures. 
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III. VPPA WELDING SYSTEM ELECTRIC POWER INPUT 


VPPA Electric Power Input can be calculated as follows: 

+ l p + av/t + + i.-AV.-t, + I p -AV -t. 

^totil (3—1) 

t + + t. 

where I B + = Straight Polarity Main Current 
I p * = Straight Polarity Pilot Current 
I B ~ * Reverse Polarity Main Current 
Ip" = Reverse Polarity Pilot Current 
Equation (3-1) is the total electric power input of the system which is 
the summation of power input for the plasma arc jet in the gap between 


electrode and orifice, power input within orifice, power input in free 


j e "t column (standoff) , and power input within keyhole of workpiece. 

The power input can be itemized, and calculated as follows: 

(III-A) Power Input at the Electrode 

o (V + Ip + )(AV C - * E )t + + (I a -- Ip") (av a + * E )t_ 

P E — (3-2) 

t + + t. 

(HI— B) Power Input in the Gap Between Electrode and Orifice 
(Joule Heating) 
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(III-C) Power Input Within Orifice 

This includes power input through the work function of metal 
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(3-4) 
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(IXX-D) Power Input in the Free Plasma Jet Column (Stand : 
(Joule Heating) 


(3-5) 
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i m P°rtant to note that shielding gas surrounding th i ' *'.vin 
jet, is induced into the plasma jet column due to the flow gene . •- ’ y 

the plasma jet stream. Power of shielding gas introduced into the 
plasma jet, P, d , shown in Equation (4-29) , shall be included as a p, r t 
of the plasma jet power input, in addition to the power input rl . y ; v,. 
Equation (3-7) . 


(III-E) Power Input in the Keyhole of Workpiece 

This includes power input through the work function of c fc vi 
reverse polarity potential rise 

D + *a ( f R “ ^ v ) 

(3 “ 

and power input through anode and cathode drops 


12 



( 3— 8b) 


I m + ^V 4 t + + I a ^ v ct. 


v, 2 


t + + t. 

(III-F) Total Power Input to the VPPA Welding System 

Total power input of VPPA welding system, shown in Equation (3-1) 
shall be equal to the summation of power input shown in Equation (3-2) 
to (3-8) and plasma jet induced shielding gas power input, P id , namely 
Ptotal “ P E + Pj,l + P*,l + P*,2 + Pj,2 + Pj,3 + Pv,l + Pv,2 + P.d < 3 " 9 ) 
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IV. POWER loss CALCULATION OF VPPA WELDING SYSTEM 
(IV-A) Electrode Heat Conduction Loss 

A one-dimensional formulation of the electrode heat conduction 
loss is given by 



where 

Q e = Electrode Heat Conduction Loss 

k E = Thermal Conductivity Coefficient of Electrode 

A e = Cross-sectional Area of Electrode 

L e = Length of Electrode (from tip to water-cooled collet) 

~ ^*0 = Temperature Difference of Electrode Between Tip and 
Collet Position 

Radiation Heat Loss Between Electrode and Orifice 
Plasma arc heat radiation loss within the gap between electrode 
and orifice can be computed from the following equation (Evans and 
Tankin, 1967; Cram, 1985): 


Qo(W) 


= v 0 (m 3 ) (4 X 10 10 


' T(°K) 
. 15000 



T(°K) 


15000 


U 6 1 


where 

Q g = Heat Radiation Loss Within the Gap Between Electrode 
and Orifice (W) 

V c = Plasma Arc Jet Volume Within Gap Between Electrode 
and Orifice (m 3 ) 

T = Absolute Temperature (°K) 

(rv-C) Power Loss Within the Orifice 


(4-2) 
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Power Loss within the orifice is comprised of the following 
items : 

(IV-C-1) Heat Convection Loss Within the Orifice 

Heat convection loss within the orifice can be computed from the 
following formulation of convective heat transfer model (Hsu and 
Rubinsky, 1987) : 

Nu = 0.2233(X+)-° ,745S + 3.66 ( 4 - 3) 

where 

Xi 

X* = — • 

R 0 (Re • Pr) m 

Here 

^ — A^ci^l Distance, (Distance from the entrance of orifice 
to the half distance of the segment of interest) 

R 0 = Radius of Orifice 

Re = Reynolds Number of Plasma Arc Fluids 
Pr = Prandtl Number of Plasma Arc Fluids 
Nu = Nusselt Number of Plasma Arc Fluids 
Subscript m = mean values 

convection loss within the orifice, thus, can be computed from 
the Nusselt number obtained from Equation (4-3) , namely, 


Q*,c = Nu • A* 

where 



Q* , c = Heat Convection Loss within the Orifice 
A* = Heat Transfer Area of Segment of Interests 


(4-4) 
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D = Diameter of Orifice = 2 R 

Jc = Thermal Conductivity of Plasma Arc 

c p = Constant Pressure Specific Heat of Plasma Arc 

= Average Enthalpy of Plasma Arc through the Orifice 
h, = Average Enthalpy of plasma in Orifice corresponding 

to the Orifice Wall Temperature 

All thermal properties used in this model are determined at the 
temperature corresponding to the average enthalpy of plasma arc. 
(XV-c-2) Direct Power Loss to Surface of Orifice 

Direct power loss to the orifice surface includes P, ( , orifice 
work function induced power input; and part of P, <2 , orifice anode 
drop induced power input, namely 


•M , G 


= P, 


,1 + 


(4-5) 


where , 

Q M/G = Direct power loss from the orifice surface 

* = Fraction power generated by anode drop lost to orifice 

(W)P "'* = Plasma arc power input from power generated at orifice 
anode drop 

(XV-C-3) Electron Collision Power Loss Within Orifice 

Power loss due to free electrons entering the surface of orifice 
can be calculated from the following equation: 


Qx,e * ~ ~ KT 

e \2 

where 

/ • ~ Free electron power loss to the orifice 
e = Electron charge = 1.602 x 10~ 19 coulomb, 

/mole 


(4-6) 
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K * Boltzman constant = 1.38 x 10" 16 ,r Vioi*K 
T * Plasma arc temperature; 

and I is the electric current which can be computed as follows: 


+ yt- 

t f + 1_ 


(4-7) 


Total power loss in orifice, Q„, is the summation of Equations (4-4) , 
(4-5) , and (4-6) , namely, 

Qx = Qx,c + Qx,g + Qx,e (4 ~ 8) 

(XV— D) Power Loss from Standoff Column 

Power loss from the standoff column includes radiative heat 
loss, Q g/r , and convective heat loss, Q 5/C . 

Calculation of radiative heat loss is similar to Equation (4-2) 
for that in the gap between electrode and orifice, namely 

f T ( °K) I 16 /( f T(°K) l 16 ) 


Q s , r (W) = v s ( m3 ) ( 4 x lol °) 


T(°K) 

15000 


15000 


(4-9) 


where V g denotes plasma arc volume within the standoff column; and T is 
the absolute temperature corresponding to the average value of plasma 


arc enthalpy. 

Radiative heat loss is strongly dependent on fluid temperature. 
This and what follows require trial and error computations. 

To compute the convective heat loss, Q S/C , from the standoff 
column, jet stream mixing theory is employed. Turbulent mixing 
between plasma jet and the surrounding shielding gas is shown 
schematically in Figure 4a. r t stands for the inner boundary of jet 
boundary layer, and r 2 expresses the outer boundary of jet boundary 
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layer. Both are functions of axial distance x from the outlet end of 
the orifice. The inner jet core region, surrounded by ri (x) is a non- 
perturbed zone of jet core. The jet boundary layer lies between r, (x) 
and r,(x). jet stream parameters in the jet boundary layer, 
including velocity, temperature and species concentration, vary from 
original (pure plasma gas jet) parameters at r - r, to ambient (pure 
shielding gas) parameters at r = r, . The jet core extends to x - x„ . 

As x > x„ , the jet core disappears and the jet boundary layer occupies 

the entire extent from r = o to r = t • .. . , 

or r 2 . The initial region from x = o to x 

= x H is the region of jet mixing. Boundary thickness at transition x = 

X H is r 2 , H • 

Based on the formulations of Abramovich ( 1963 ), the boundary 

y within the initial region of jet mixing is governed by the 
following equations: 


R o " r t 


cx 


1 - m 


1 + m 


0.416 + 0.134m + 0.021 — 

R. 


(1 + 0.8m - 0.45m 2 ) 


cx fl - m 


\ 1 + mj 


(4-10) 


R o ~ r i 


r 2 - r : 


- 0.416 + 0.134m + 0.021 — fi’d + 0 .8m - 0.45m 2 ) 


(4-11) 

Where R„ denotes the radius of orifice; r t , the radius of inner 
boundary of jet boundary layer; r, , the radius of outer boundary of jet 
boundary layer; and x, the axial distance from the outlet end of 

orifice to the point of interest, m - u H /u 0 , where u H is ambient fluid 
velocity; ;and u„, the original jet flow velocity inside jet core, c 
IS empirical constant which is 0.27 when u »u u . 

0 H * 
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Substituting the values of m = 0 and 


c= 0.27 in Equations (4-10) 


and (4-11), we have 
R„ - r , 


x 

R„ - r 


= 0.27 


0.416 + 0.00567 


i 


r 2 - r t 


= 0.416 + 0.021 


( r 2 - 


(4-12) 


(4-13) 


From equation (4-12) , we can solve for the jet inner boundary r 2 in 
term of the axial distance x which is given by 


r i x 

— = 1 - 0.1123 0.0015 


(x \ 


VRn/ 


(4-14) 


Similarly, from Equation (4-13) , we can determine the boundary layer 
thickness, r 2 - r 2 , i.e.. 


1 

r 2 - r i = - { -19 . 81R 0 + [ 392 . 44R„ 2 + 190.48 (R Q - r x ) R„ ] 1 / 2 } 

(4-15) 

To determine the axial distance of jet initial region x H , we can 
solve Equation (4-14) by substituting r 2 = 0, i.e., 

». = 8 (4 ’ 16) 

In our case, the radius of orifice R 0 = 1.5875 mm, and the axial 
distance of the jet initial region x H = 8(1.5875) = 12.7 mm. In the 
present study, the range of standoff distance considered is between 1 
to 10 mm, so that all standoff values considered in the present study 
lie within the range of the initial region of jet mixing. Let us 
consider velocity, temperature and species concentration profiles 
within the range of jet initial region as follows: 
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U 0 - u 


U 0 - U H 


= (1 - - S ) 2 


T 0 - T 


T o ~ T„ 


= 1 - r, 


(4-17) 


(4-18) 


Pi 


P , + p 

Pi S d 




(4-19) 


where u 0 and T a denote original velocity and temperature of jet core, 
respectively; and u H and T H , the ambient velocity and temperature at 
the outer layer of jet boundary layer, respectively. Dimensionless 
coordinate parameter , - (r 2 - r)/(r 2 - r x ) . Psd and Ppl denote masses 

of shielding gas and plasma gas per unit of volume within the jet 
boundary layer, respectively. 

In the present model, the effect of the mole fraction of 
shielding gas induced by jet mixing into the plasma column has been 
fully taken into consideration. 

It is convenient to employ partial pressure P pl and P Jd , rather 
than density in the calculation. By using equation of state. 
Equation (4-19) can be rewritten as 

P. 


id 


“pi 


5 pl 


1 

\ ’l 


(4-20) 


Since 

P Pl + p id = Ph 

p d H (4-21) 

We can solve partial pressure P pl and P sd from Equations (4-20) and (4- 

21) as follows: 
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(4-22) 


P 


P 


and 


Ph 



(1 ] 


1 


M.d 

i ^ J 



P.d=PH-Ppl (4-23) 
where P p x and P Jd denote partial pressures of plasma and shielding 
gases, respectively; and M p j and M sd stand for molecular weights of 
plasma and shielding gases, respectively. For example, M p L = 39.95 
kg/k-mol for argon gas, and M sd = 4.003 kg/k-mol for helium gas. 
Also, P H is ambient pressure of gases. 

Plasma and shielding gas density, p pl and p sd , can be determined 
from Equation (4-18) for temperature prof ile and Equations (4-22) and 
(4-23) for partial pressure of plasma and shielding gases, i.e.. 


p 


pi 


p Pi M Pi 


(4-24) 


r sd n sd 

Psd (4-25) 

R T 

where R denotes the universal gas constant which is 8315 J/kmol*K. 
Furthermore, plasma and shielding gas enthalpies can be obtained from 
the table at their corresponding values of temperature. Finally, jet 
stream power P Jet can be calculated from the following integration. 


Jet 


(r. 


> o 


2 7rru 


( P p i h p l 


+ />sd h sd> dr 


(4-26) 


It is noted that Ppl = (/> pl ) 0 / (p sd)o " °# h P i = ( h P i)o/ and 

u = u 0 within the region of jet core. Thus, Equation (4-26) can 
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be rewritten as follows: 

P J • t = 1 2u o ( Pp i ) o (kp i ) o + 

f r 2 

+ 2*ru (p pl h pl + />, d h sd )dr (4-27) 

J r i 

At the outlet end of standoff column, the total power of plasma 
jet can be expressed as 

p/ 3,o " P 3,l + p j , 3 - Q S/ r + p *d (4-28) 

where P 3 / j denotes the plasma jet power at the inlet end of standoff 
column (or at the outlet end of orifice) ; Pj^ 3 , the plasma jet power 
due to Joule heating within the standoff column; Q s # r , the radiative 
heat power loss from the standoff column; and P sd , the induced power 

from surrounding shielding gas into plasma jet boundary layer, 
namely, 

p sd = m jd , h (4-29) 

Here, h s<J/H denotes shielding gas enthalpy corresponding to the 
ambient temperature of shielding gas; m sd , mass flow rate of shielding 
gas 

2XrU p ad dr (4-30) 

In this calculation, total power of plasma jet P' 3/0 shown in 
Equation (4-28) , contains the power within the whole cross section of 
jet boundary layer where temperature varies radially from more than 
10 K at the inner layer of jet boundary layer to ambient temperature 
( sa Y 298 K) at the outer layer of jet boundary layer. The welding 
requires sufficiently high temperature difference to 
maintain heat conduction and to promote sufficient heat flow to melt 
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the workpiece metal. More important than promoting heat flow by 

convection is the maintenance of a conducting plasma path for current. 

Spot welding electrodes do not have to be hot at all to provide a molten 

weld nugget by resistance heating. Resistance heating due to plasma 

sheath voltage drop and work function drop are more important (Eckert 
and Pfender, 1967) . 

Let R c be a half size of the crown width. This is also equivalent 
to the radius of the inlet end of keyhole (See Figure 4b) . 

As R c < r 2 . 


p 3,o 2irru ( Ppl h pl + pSd h 8d )dr 

and as R c > r 2 , 

2 

o 


(4-31) 


- f 

a 


2*ru ( ,,, h fl t fti h ld )dr + *(R t * - r 2 ‘) u H fti ( 4 - 32 ) 


The difference between P' 3 , 0 , calculated from Equation (4-28), and 
P a ,„, calculated from Equations (4-31) or ( 4 - 32 ) , is the convective 
heat loss or power increment due to injecting shielding gas Q tc , 
W ^thin the standoff column, i.e., 

Qs 'C =P ' 3 '0 - p 3,0 (4 _ 33) 

In short, workpiece heating is mainly contributed by the plasma 
jet, comprising a mixture of plasma and shielding gases, with the 
power ofP, #0 . Throughout the keyhole all plasma gas properties used 
in heat transfer calculations are based upon the mixture of plasma and 
shielding gas established at the inlet of the keyhole. 

A detailed descritption of how to determine radius R c will be 
disscussed in Section VI of this report. 
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(rv-E) Total Workpiece Power Loss 

Total power loss in the workpiece is comprised of the following 
items : 


( XV-E_1 ) Convective Heat Loss in Workpiece 

Local power loss to the workpiece from heat conduction can be 
computed from Equation (4-34), namely, 

^ fkp Urn — h w ) 

D c. 


Q w ,c = Nu z • A v 
z 


(4-34) 


where Nu z is the Nusselt number along the axial direction; A w stands 
for heat transfer area per unit length along the axial direction of the 
workpiece (A v * irR, ) i.e. the plasma contact area in the keyhole taken 
to resemble a finger ring comprimising somewhat more than one half of 
the circumferencial area of the keyhole; R is the radius of plasma 
contact area in the keyhole which varies with respect to axial 
distance z and is determined by the calculation; and D is the diameter 
of the keyhole equal to twice the effective radius of plasma contact 
area (2 R) . The rest of symbols are the same as used in the orifice. 

The Nusselt number is computed as in Equation (4-3), used in 
orifice, which is 

Nu z = 0.2233(z + )-°- 7 «55 + 3.66 f 4 --^ 


where 


z 1 

z + = 

R (Re-Pr) m 

Here R denotes effective radius of plasma contact area in the keyhole 

at the axial location of z. z is the distance between the entrance of 

keyhole to the half thickness of the section of interest in the 

workpiece considered. The rest of the symbols stand for the 'v-up 
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definition used in the orifice. 

Total convective heat loss in the workpiece is 


CL a 


Qv,c 


dz = 


NUj A v 


k 

D 


( h plasm« " ^v) 


dz (4-36) 


(XV— E— 2) Direct Power Loss from the Surface of Keyhole 


Direct power loss to the keyhole surface includes P H / : , the power 
generated by the workpiece work function; and part of P Vy 2 , the power 
generated at the workpiece anode and cathode drops, namely, 

Qv,G “ Py,l + **v,2 (4~ 37 ) 


where 

Q v c = Direct power loss to the keyhole surface 
/ = Fraction of power generated by anode and cathode drops 

that is transmitted to the keyhole surface 
(l-^)P W/2 = Power transmitted to the plasma from that generated 
at anode and cathode drops 
(XV— E— 3) Electron Collision Power Loss in Keyhole 

Power carried by free electrons into the surface of the keyhole 
can be calculated from the following equation: 


Qv , a 


I 

e 


(3 ' 

- KT 
2 


(4-38) 


where Q W/ « denotes power carried by free electrons in to the keyhole 
wall; T, the mean electron temperature entering the wall taken as the 
average temperature between plasma arc and wall of keyhole; and I, the 
straight polarity component of the transmitted current, which can be 
computed from the following equation: 
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I 


( 4 - 39 ) 


t t + t_ 

Thus, total power loss to the workpiece, Q w , is the summation l. o ;sa 

computed by Equations ( 4 - 36 ), ( 4 - 37 ), and (4-38), namely 

Q V " Qv,c + Qw, 0 + Qv,* (4-40) 

Power Discharge Loss 

The power carried away by the plasma jet discharge from fir. At 

end of keyhole, Q d t 8 , can be computed from the plasma arc enfch U • f ,,i 
mass flow rate, i.e., 

Qdis = (h pl m pl +h sd m sd ) dis (-lAi.) 

where the values of h p L m p x and h sd m, d are the results of integr ;-. v- r c 
the cross-section of the outlet end of the keyhole from r - 0 to r = r . 
R r is the radius of the outlet end of the keyhole. A detailed 

description of how to determine R r will be discussed in Section Vi. cn 
this report. 

Obviously, summation of power loss throughout the entire 'Vi'iil. 
welding system is equivalent to the total power input of the system, 
p tot»i» namely 

p tot*i - Q e + Q c + Qm + Q s + Q h + Q dis (4-42) 

Calculation of power loss at each part of VPPA welding system 
will be illustrated in Section VII of this report. 
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V. CALCULATION OF PLASMA ARC ENTHALPY 

an example, the density of argon gas at a pressure of 1 

atmosphere and a temperature of 298 °K is 

PM, r 101320 (N/m 2 ) • 39.95 (kg/kmol) 

P ~ = 

p R T 8315 (N-m/kmol °K) • 298 (°K) 

= 1.6336 (kg/m 3 ) 

where 

P = Atmospheric pressure 
M Ar = Molecular weight of argon 
R = Universal gas constant 
T = Absolute Temperature 

The mass flow rate of argon gas in terms of volumetric flow rate can be 
computed as 

1 Hi ) 

m - (V ft 3 /hr) hr/s m 3 /ft 3 p 

.3600 J (35.3147 J Ar 

1.6336 

= (V ft 3 /hr) 

(3600) (35.3147) 

= 1.285 • 10' 5 (V ft 3 /hr) (kg/s) 

" 1-285 • 10' 2 (V ft 3 /hr) (g/s) 

(V-A) Plasma Arc Enthalpy Within the Gap Between Electrode and 
Orifice 

Let us assume that entrance and exit ends of plasma arc enthalpy 
in the gap between electrode and orifice are h. , and l\, 
respectively. 

The plasma arc power entering the gap, P 1/t , is the power 
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generated at the electrode, P E , minus the power loss by electrode heat 
conduction, Q E , namely 

P l,l - P E - Qe (5-1) 

Plasma arc enthalpy at the entrance end of the gap, thus, can be 
calculated as 



m 


The plasma arc power at the exit end of the gap, P 1/0 , has been 
augmented through Joule heating, Pj, lf and diminished through 
radiative heat loss, Q G . It is 

p i,o = p i/i + p jm “ Qg (5-3) 

Thus, plasma arc enthalpy at the exit end of the gap can be calculated 
from the expression 


h 


1 , 0 



(5-4) 


m 

(V— B) Plasma Arc Enthalpy Within Orifice 

Let us assume that the enthalpy of the plasma at the entrance and 
exit ends of the orifice are h 2/i and h 2 0 , respectively. 

The power, P 2/ i/ and enthalpy, h 2/1 , of the plasma arc at the 
entrance end of the orifice are identical to the power, P x 0 and 
enthalpy, h x / 0 , respectively, of the plasma arc at the exit end of the 
gap between electrode and orifice. Thus, 

P 2 / i = p i / o 
^2 , i = ^i,o 

The plasma arc power at the exit end of orifice, P 2/0 , is 
increased by part of the power generated at the orifice anode drop, (1- 


(5-5) 
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*)P„ /2 , and by all the power generated by Joule heating within 
orifice, P J/ 2 , and is reduced by convective heat loss to the orifice, 
Q* /C , and free electron power loss to the orifice, q . ^ v . 5 , 


2 / 0 


" P *' 1 + (W)P *^ + - Qm, c - Q„, e 


: • 6) 

The corresponding plasma arc enthalpy at the exit end of orifice is 
P. 


1 2 , 0 


2 / o 


m 


: c > -7) 


(V-C) Plasma Arc Enthalpy Within Standoff Column 

Let us assume that the plasma enthalpies sat at the eiitx and 
exit ends of the standoff column are h 3 # t and h 3/0 , rec;^ y . 

The plasma arc power, P 3 # j , and enthalpy, h 3/ lr at the ta*: . se 
end of standoff column are identical to the plasma arc power, P 2 # # , nnd 
enthalpy, h 2/0 , respectively, at the exit end of orifice. Thus, 


P 3 , i ~ P 2 , 0 
h 3 , i = h 2 , o 


(s-e) 


Plasma arc power at the exit end of standoff column, % ■;*, 
increased by Joule heating in standoff column, P J/3 , and indeed 
shielding gas mixing power, P, d , and reduced by radiative heat loss, 
Q S/r , and convective heat loss, Q S/C , from the same column. Thus, 


3 , 0 


= P 


3 , i 


+ P 


J / 3 


+ P, 


d “ ®s,r “ Qs, 


(5-9) 


The corresponding plasma arc enthalpy at the exit end of standoff 
column is 

*pi h P i + ®«d h sd 


*3,0 


fflp j + m 


* d 


( 5 - 10 ) 


where the values of i pl , ® 8e |, ^pi^px, and m s<J h sd are the results of 
integration at the cross-section of outlet end of standoff column from 
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r * 0 to r - R c where R c is the half crown width. 

(V-D) Plasma Arc Enthalpy Within Keyhole of Workpiece 

Let us assume that the plasma enthalpy at the entrance and exit 
ends of the keyhole of the workpiece are h 4/ i and h 4/0 , respectively. 

The plasma arc power, P 4 t j , and enthalpy, h 4 t j , at the entrance 
end of the keyhole of workpiece are identical to the plasma arc power, 
P 3 ,o' and enthalpy, h 3 t 0 , respectively, at the exit end of standoff 
column. Thus, 


p 4 , l p 3 , o 
^4 , i — ^3 , o 


( 5 - 11 ) 


Plasma arc power at the exit end of the keyhole of workpiece, 
P 4 , o • i s increased by a portion of the power generated at the workpiece 
anode and cathode drops, (1-/) P v t 2 , and is reduced by convective heat 
loss to the keyhole surface, Q V/ e and the loss of hot electrons to the 
keyhole , Q w / e . Thus , 

P 4,o * P 4,i + (W) p y ,2 ~ Q v ,c - Qw,. (5-12) 

The corresponding plasma arc enthalpy at the exit end of the 
keyhole of the workpiece is 
m p ih p i + m Jd h, d 


l 4 , o 


m 


pi 


m 


( 5 - 13 ) 


3 d 


where the values of m p 1 , m S(d , m pl h pl , and m s< jh s<J are the results of 
integration at the cross-section of outlet end of keyhole from r = 0 to 


r = R r . 

(V-E) Plasma Arc Enthalpy Discharged from the Workpiece 

Plasma arc power, Qjij, and enthalpy, h d j s discharged from the 

keyhole of workpiece is identical to the plasma arc power, P 4 0 , 

enthalpy, h 4/0 , at the exit end of keyhole, namely 

30 


and 



(5-14) 


Q d 1 « * P 4 ; o 
l « = ^4,0 

Calculation of plasma arc enthalpy at each part of VPPA welding 
system will be illustrated in Section VII of this report. 
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VI. CALCULATION OF CROWN AND ROOT BEAD WIDTHS 


Bead width is calculated from the heat absorbed in the workpiece, 
Q h . As shown in Equation (4-40) Q H is the summation of convection 
heating (the empirical relation actually includes radiation.) from 
the plasma, "direct 1 ’ heating of the surface by current passing through 
localized surface voltage drops, and heat carried into the surface by 
hot electrons from the plasma. 

The heat input source in the keyhole during the actual welding 
process is neither a point source nor a line source along the axial 
direction in which the input energy is not uniformly distributed. To 
make a calculation closer to the actual situation, we adopted a 
calculated result at various locations of z along the axial direction 
as a heat source distribution. As the magnitude of the heat source 
per unit length along the z axis is closely related to the area of heat 
trandfer (radius of keyhole at various locations of z) , both heat 
absorbed by the keyhole per unit lenth along the axial direction and 
the radius of the keyhole at various locations of the z-axis shall be 
solved by simultaneous iteration. 

The temperature field within the confined space between the two 
surfaces (assumed nonconducting) of the workpiece is constructed out 
of moving heat sources in an infinite continuum. "Image" sources are 
added outside the workpiece space so as to produce planes of symmetry 
across which no heat flows at the location of the workpiece surfaces. 
Given coordinate z zeroed at the top of the workpiece and extending 
down into the workpiece body and locating the heat source Q H/ z at 
location z, then an image heat source above the top of the plate of 
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strength Q H/ z located at -z creates a plane of zero heat flow at z =0, 
the location of the top surface of the workpiece. The latter plane 
will remain and a further plane of zero heat flow at the site of the 
bottom surface of the workpiece will be created if an infinite array of 
image heat source pairs is constructed along the z-axis, each heat 
source pair of strength Q H/Z spaced at ± z about an infinite array of 

points spaced at 2L 4 and including z = 0. See Figure 5 for 

illustration. 

The temperature T c at location (x c , y c , z c ) within the workpiece 
due to the welding power input to the workpiece is then: 

Q'w,z dz f u(Rj + x c ) 

exp - 

4ttKRj ( 2a 

where 

T 0 = Ambient temperature of workpiece 
K = Thermal conductivity of workpiece 
u = Weld speed (in the x-direction) 
a = Thermal diffusivity of workpiece 
Rj = Distance between the i th point heat source and the point 

of interest ! 

i 

R t is computed as follows: | 

i 

*i - -J x c 2 + Yc 2 + (2i “ 2 C ) 2 (6-2) I 

| 

1 + (- 1) 1+1 • 

Zi = [i + ] L 4 + (-1) 1 z (6-3) i 

2 , 

Weld bead width is determined by locating the largest value of y c 

for which T c = T m , the melting temperature of the workpiece, at z = 0 



( 6 - 1 ) 
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(crown) and z c - L 4 (root) . The weld widths are twice the values 
of these y c 's. A converging computer search routine is used to 
make the computation. 

One approximation is used in the computations. The s umm ation 
is truncated. 35 source points are used, 17 image sources on each 
side of the weld heat source. The values of z. are generated by a 
recurrence relation: 


Z \ ~ Z f*2 + 


where the starting values are: 


z 0 = z 


z, 2 z 


(6-3) 


(6-4) 

(6-5) 


where 0 < z < L, 


z o z ) the weld heat source position located in the workpiece. 
After truncation. Equation (6—1) becomes 


i- -17 


Q v , z dz R, + x 

— — r — exp [ — u ( — 

4xici? 1 2a 


( 6 - 6 ) 


As we indicated earlier for the computation of heat source 
distribution along the z-axis, the temperature distribution within 
the workpiece can be obtained by integrating Equation (6-6) which 


rp rp C 0«T, g^Z ( R± + 

T '- T ° = ( 


V' r r .. IXj + x c) 


E C wv, r 


(6-7) 


Adopting temperature T c to be the melting temperature T m , one can 
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compute the melting isotherm surface from Equation (6-7) . The 
geometric distribution of the keyhole can also be obtained for the 
maximum radius of the melting widths at various locations of z. 
Crown and root widths, 2R C and 2R r , at the upper and lower ends of 
the workpiece can be calculated by a similar method. 

Equation (6-7) is used to compute the temperature distribution 
within the workpiece in which Q' uz , denotes the power conducted 
away through the workpiece, or the strength of heat source per unit 
length at location z. This strength of the heat source in the 
workpiece is the heat absorption Qy z contributed from the plasma 
jet, minus a power loss due to the radiative heat transfer from the 
surfaces of the workpiece Q R . In the actual numerical integration 
of Equation (6-7), Q R is included in Q' uz at various locations of 
z by a method of weighing. 

Since the distributions of Qy G from Equation (4-37) and Qy e 
from Equation (4-38) inside the keyhole are hard to determine, we 
assume the parabolic type of square approximation (see Figure 6 for 
schematic expression of the workpiece power distribution) , 

Q z = a (1 - -y-) 2 Q (6-8) 

L t 

where Q denotes either Q w G or Qy e * a is a constant which will be 
determined from the integration of Equation (6-8) based on the 
value of Q in which the total value of either Qy G or Q w e shall be 
conserved. This is 
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£ 4 1*4 

0 = f Q z dz = fail - JL) 2 * Qdz 

0 0^* 

which leads to a = 3/L 4 , and Equation (6-8) becomes 



Radiation losses from both sides of the workpiece surfaces are 
computed. The radiative heat transfer loss can be computed as 


0 R = ZApaTf ( 6 - 10 ) 

where Aj denotes a small element i of surface area; €j is emissivity 
of the workpiece; T { stands for the average temperature (°K) of the 
surface area element A ( ; and a is Stefan-Boltzman constant (= 5.67 
x 10' 8 W/m 2, k 4 ) . 

Workpiece power absorption per unit length along the z-axis 
due to melting heat can be computed from the following relation: 

Qm.z * ( 6 - 11 ) 

where 

u = Traveling speed of welding 
p„ = Density of workpiece 
Q f = Melting heat of workpiece 
d z = Bead width at position z 

Based on this analysis, power transfer per unit length in 
workpiece Q' u z can be computed as follows: 


2 ~ Q*f.c, + T ^ ~ r ^ + - On, r “ ( ~!t' ~ ) i2) 

"^4 *“4 Uu 


36 



The total net power absorbed by workpiece is then 


Ql - / Ql. t dz ( 6 - 13 ) 

0 

Substituting Equations (6-12) in (6-13), we have 

OL = 0 V - 0 R - 0 a ( 6 - 14 ) 

where 

~ f @w, c,dz + Q v ,G + Qw. a 
0 

- ®V.C + Qff.G + Q*. a 

is similar to Equation (4-40) , and Q R is from Equation ( 6 — 10 ) . 
Here, 

*■« 

@a ~ j Qm. z^Z 
0 

is the total melting power. 
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VII. CALCULATION OF PARTIAL PENETRATION OF THE WORKPIECE 
In Section VI, we have discussed the weld widths and power 
loss in the workpiece. An estimation of full penetration or 
partial penetration shall be made before we can proceed the 
calculation of Section VI. Before the keyhole is formed for 
P ar "ti a l penetration, the plasma jet impinges perpendicularly on the 
surface of the weld pool, and is deflected across the weld surface 
radially away from the stagnation area. Thus, the heat transfer 
from the plasma jet to the surface of the weld pool and the -model 
of the heat source used to calculate the temperature field 
workpiece are different from those of the weld process in td---' '• 
keyhole has been formed in the workpiece. 

(VII-A) Weld Pool Convective Heat Transfer 
For the case of partial penetration, a concave shape melting 
puddle is formed while the plasma jet impinges on the surface of 
the workpiece. At this moment, there is a very complicated f'i.ox-/ 
pattern developed within the melting puddle of the molten liquid 
me ^ a l* It is anticipated that heat transfer process, which causer 
the convective heat flow from the plasma jet to the workpiece, i:: 
also very complicated. It is hard to adopt a very suitable 
mathematical formulation to describe convective heat loss from the 
plasma jet to the workpiece. In this study in the determination of: 
partial and full penetration of the workpiece, an axial symn dry 
circular plasma jet, which impinges perpendicularly on the surface 
of the weld pool of the workpiece, is assumed to compute the 
convective heat loss problem. 
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The power transferred to the workpiece by convection from the 


plasma arc jet can be written as: 


- hA v (HpXagna ~ H w ) / 


(7-1) 


where 

= Plasma arc jet impingement surface area 

= JTRj 2 

Rj ** Outer boundary radius produced by plasma arc jet on the 
surface of the workpiece (see Figure 4b) 

H piasm» = Average enthalpy of the plasma arc jet on the 
stagnation area A w 

H M = Enthalpy of the plasma arc jet at the temperature of 
the weld pool 

= Specific heat of the plasma corresponding to a 

temperature T a at which the enthalpy of plasma is equal 
to the mean enthalpy at the temperature of the arc jet 
and the weld pool 

h = Convective heat transfer coefficient 

h can be computed from the following relation (Metcalfe and 
Quigley, 1975) : 


N u = 0.78 R g °- s P°- 3:i 

where the dimensionless parameters are defined as follows: 


(7-2) 


= hi2PJ 
k 


0 . PW2 R,) 

K = 

8 p 
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Here, v is the average velocity of the plasma arc jet within the 
outer boundary on the surface of the workpiece. The rest of the 
parameters are the same as that defined earlier. In this 

calculation, the parameters of the plasma are determined based on 
the average enthalpy within the cross sectional area of radius Rj. 

Total power absorbed by the weld pool or the workpiece, Q w , 
shall be the summation of Equations (7-1), (4-37) and (4-38), i.e., 

Q" = Qv.c + Qw.G * Qv.a ( 7 - 3 ) 

The Model of Temperature Pield and Partial Penetr't r 
In the determination of the formation of keyhole from the 
plasma arc jet impingement on the surface of the workpiece, it is 
assumed that there is a point source with an intensity of Qy vht b 
moves with a speed of u acting on the surface of the workpiece. 
Point source model, proposed by Rosenthal, is adopted to compute 
the temperature field of the workpiece created by this moving point 
source. The depth of the weld pool, z m , created by the rnovii: 
point source can be computed by the constant temperature surface of 
the workpiece molten temperature which is the solid-molten pool 
interface of the workpiece weld pool. The impingement of the high 
speed plasma jet creates an intensive disturbance on the molten 
liquid of the workpiece and forms a concave hole on the surface of 
the workpiece. It may create a distributed heat source from the. 
surface of the workpiece to the depth of the weld pool if the 
impingement of the plasma jet can reach the bottom of the molten 
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pool. To simplify the calculation , we may assume that there is a 
uniformly distributed heat source moving with a speed of u, and 
with an intensity of Q u /Z m (See Figure 7) . The temperature field 
of the workpiece can be computed based on the effect of the moving 
uniformly distributed heat source. Determination of the depth of 
the partial penetration and the formation of the keyhole can be 
resulted from the distribution of the constant temperature surface 
of the workpiece molten temperature. 

The temperature field produced by the point source . t the 
surface of the workpiece moving with uniform velocity of u 'r.\ ie 
a moving coordinate with a speed u moved with a point 
Location of the point source is at the origin) is 


T c ~ T 0 


iR > * x ° )] 


- E 


(7—4) 


where 

k = Heat conductivity of the workpiece 
a = Thermal diffusivity of the workpiece 
u = Moving speed of the point source 
Rj = [* c 2 + y c 2 + (z c - z ,.) 2 ] 1/2 

x » Y/ z = Location of the point of interests 

z,- = (i + 1 ± f-ll w 1’L 

2 

L 4 = Thickness of the workpiece 

Here, 

Qu‘ - Q« - Q. 


Qy can be computed from Equation (7-3 ) , and 


41 


ORIGINAL PAGE fS 
OF POOR QUALITY 



Qh, = Latent power of the workpiece 


3 f Q e upd z dz 

O 

where 

Q f = Latent heat of the workpiece 
p = Density of the workpiece 

d t = Width of the molten liquid of the workpiece at the 
corresponding location of z 

Put T c = T m (molten temperature of the workpiece) in Equation (7-4) , 
one can solve the constant molten temperature surface of the 
workpiece. The depth of the workpiece molten pool under the effect 
of the heat source, thus, can be determined. 

The temperature field created by the moving uniformly 
distributed heat source with intensity of Qy'/Z,,, located from z = 
0 to Z m is 

m 

2 

/-r^ lr exp 1 [ --& ,dz (7 ' 5> 

i--« o l 

where 

R x = [ x c 2 + y c 2 + (z i - z c ) 2 ] l/2 

z t = (i + - 1 - 1 - ] -L 4 + (-l) i z 

o*z*z m 

Qj = Q w - Q m 
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*1 

O m = f Q e upd z dz 
o 

z t = Molten metal depth of the workpiece under the effect of 
uniformly distributed heat source 

The determination of full penetration or the formation of the 
keyhole is dependent upon whether z t is equal to the workpiece 
thickness L 4 . In the realistic calculation, adoption of 35 terms 
of series summation in Equation (7-5) is good enough to assure* the 
series convergence accuracy of the temperature field determirrv. n. 
In other words, the adoption of -17 < i < 17 together with solving 
simultaneous equations of Equation (7-5) and the computatio \ 

with iteration of related variables are required to solve the 
temperature filed distribution. 

It is due to the fact that there is no experimental 
available for the variation of the depth of partial penetration 
with respect to control parameters of the welding processes. Th 
is no way to compare the model results with the all of the stages 
of penetration processes during the plasma arc jet welding. In 
this study, model comparison with the experimental results of 
partial penetration and the possibility of the formation of the 
keyhole are made. 

Equation (7-5) is used to calculate the probability of full 
penetration as the lower surface of the workpiece reaches the 
molten temperature. For the case of partial penetration, it may be 
due to the fact that the power supply is too small and/or the weld 
speed is too high. So long as the keyhole is formed, Section VI 
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Will be used to compute the power absorbed by the workpiece in the 
keyhole, the weld widths, etc. 
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VIII. CALCULAITON OF CROWN AND ROOT HEIGHTS 


(VIII— A) Volume Conservation Equation of Materials 

Assuming that reinforcement profiles of crown and root beads are 
paraboloid, as shown in Figure 8. These paraboloid profile can be 
expressed in the following equation: 

(2y\ n 


Z = h. 


1 - 


<v 


( 8 - 1 ) 


where 


h c = Crown height 
d c = Crown bead width 
n = Numerical value greater than 1 
The cross-sectional area of the crown reinforcement is 


A e 



'2y' 

n “ 

UJ 



dy 


n 

h c d c 

n + 1 

Similarly, cross-sectional area of the root reinforcement is 

n 

A- = h r d r 

n + 1 


where 

h r = Root height 
d r = Root bead width 

Volume conservation requires the equation of the material volume 

provided by welding wire to the reinforcement volume crown and rocn 

beads plus any gaps of holes in the weld site prior to welding: 

45 



n 


n 


(8-2) 


where 


n + l 


d c h e 


n + l 


d r h r 


TV V y 

- dv 2 

4 U 


W 


1 


d v = Diameter of wire 
V v = Feeding speed of wire 
u = Welding speed 
Wg = Initial gap at the weld seam 
L 4 = Thickness of workpiece 
(VIII— B) Force Equilibrium Equation 

It should be noted that for the sake of tractability the force- 
equilibrium analysis has been greatly simplified. Longitudinal 
curvature of the weld puddle surface has been neglected as have the 
effects of Marangoni circulations and thermal variation of surface 
tension. There is a great deal of room for improvement of this first, 
very rough analysis. The cross-sectional areas a, b, and c, defined 
in Figure 9 which are given by 


d r /2 



= h c d r - 


1 k „ ^ n+l 

h c d c 


n + 1 




(8-3) 


A b = d r • L 4 


(8-4) 
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The slopes of the reinforcement surfaces at points B (crown) 
and D (root) , shown in Figure 9, are given by 


(8-5) 


dz 

tan 9 c = — 
dy 


2nh e 




y=d r /2 




n - 1 


and 


tan 0 r = 


dz 


dy 


2nh r 


|y=d r /2 d r 

It is due to the fact that both 0 C and 0 r are relatively small, and 
we can assume sin0 c = tan0 c and sin0 r = tan0 r . 

With reference to the free-body diagram shown in Figure 9, the 
force equilibrium equation for the body comprising areas a, b, and c is 
given by 


h c d r - 


n + l 


h c d c 




n + 1 




+ d. 


n 


L 4 + 


n + 1 


d r h r 


= 4n 

h r h c 


n - 1 - 


.<*r 

idj 



<r - pd r 


p g 


(8-6) 


where the first term represents a gravitational force, g being the 
gravitational acceleration component perpendicular to the workpiece 
surface and directed from crown to root. The second term represents a 
surface tension force, a being the molten weld puddle surface tension. 

"third term represents a pressure force, p being the mean pressure 
acting on the weld crown surface. 

(VXII-C) Calculation of Crown and Root Heights in Horizontal 

Welding (with Consideration of Gravity Force) 

Experience suggests that the effect of external pressure p is 
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much smaller than that of the surface tension. Therefore, the p term 
will be deleted from the force equilibrium equation (8-6) . Crown and 
root heights, thus, can be found by simultaneously solving Equations 
(8-2) and (8-6) which yield the following expressions: 

* V, 
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4 <r 

1 (n + 1) — 

d r P<3 J 


" d v 2 L 4 

4 U 
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- d r .L 4 
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d c + 


4n <r 

' 1 

r dr l 
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d r 2 J 
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IdJ 
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( 8 - 8 ) 
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n + 1 

(VTII-D) Calculation of Crown and Root Heights in Vertical Welding 
(Without Gravity Effect) 

The gravity force component g perpendicular to the workpiece 
vanishes for the case of vertical welding and also for welding carried 
out in space. Ignoring both gravity force and external pressure in 
Equation (8-6) one obtains 


hr = h c 


ug 


(8-9) 


and substituting this equation to Equation (8-2) , 


h c = 


7T 

- d v 2 
4 

V v 
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n 

r* i 
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[ dr j 

n * 

n + 1 

a c + 

. 

d r 

i d cJ 



( 8 - 10 ) 
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(VIII-E) Depression of Crown Height During Solidification 

Figure 10 shows the configuration of the cross-section of puddle 
during the process of solidification. The solid-liquid interface 
lies along the melting temperature isotherm, which generally expand 
toward the hotter crown and contract towad the colder root. They come 
together, thus, first at the root as in Figure 10(B) . This leaves a 
sector of molten metal which shrinks upon solidification. The height 
of the solidified crown must fall to accomodate the shrinkage volume 
reduction as shown in Figure 10(C) . The cross-sectional area of the 
molten sector in Figure 10(B) can be approximated as 
1 

A c = - L< • d c ' (8-11) 

2 


where L 4 stands for the thickness of the workpiece; and d c ' , the crown 
width in a state of molten metal. Given a volumetric shrinkage, av/v, 
a material property, the reduction in area aA c of the molten material 
during solidification is 



f AV' 

1 

'av’ 

aA c = A C 

— 

= ~ l 4 • d c ' 



iV J 

2 



( 8 - 12 ) 

From Equation (8-5) for the cross-sectional area of the crown 
reinforcement, the reduction in the area can be related to the change 
in the crown reinforcement height Ah c ' 


aA c = a 
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In + 1 


d c 'h c ' 


n + 1 


dc'*h c ' 


( 8 ** 13 } 


Equating Equations (8-12) and (8-13), leads to 


AVI 


(V ) 
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n + 1 


d c 'ah c ' 
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which gives upon simplification 


ah c ' = 


n + l 


2n 


( av 


v ) 


(8-14) 


Using King's (1987) , value for specific volumetric shrinkage for 
aluminum during the solidf ication, av/v, is 6.7%. Therefore, crown 
height shrinkage for solidification of aluminum Equation (8-14) 
becomes 


n + l 

ah e ' = 0.067 L 4 

2n 

Detailed calculation results of crown and 
illustrated in Section IX of this report. 


(8-15) - 

root heights are 
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IX, LEADING EDGE ANGLE COMPUTATION 
The leading edge angle of the keyhole of workpiece can be 
computed from thermal power equilibrium requirements as follows: 

tan8 - P^^CT.-T,) » c, (T.-TJ I 

pj^dl a 

where 

P * Density of workpiece ( k9 /m 3 ) 
v * weld speed (m/s) 

^ = Melting latent heat of workpiece ( kJ / k9 ) 
c * ~ Specific heat of solid workpiece (Kj/kg-K) 

~ Specific heat of liquid workpiece (Kj/kg*K) 

— temperature of workpiece (K) 

T 0 = Original temperature (room temperature) of workpiece (K) 

T* = Tem P era ^ ur e of superheated molten fluid at leading edge (K) 

P u = Thermal power absorbed by workpiece (Kw) 

D = Keyhole diameter (m) or orifice diameter (m) 

^4 ~ Thickness of workpiece (m) 

9 = Leading edge angle 

Theoretical computations and experimental measurements from 
rubber molds of keyholes (Nichols Research, 1992) are found to be 
in good agreement, at least with respect to velocity dependence 
except for four measurements at low plasma flow rates, which show 
experimental leading edge angles of 6 > or = 30». in other words. 
Equation (9-1) may hold only above some threshold value of plasma 
flow rate. Above the threshold value, for reasons not yet 
understood a "beam absorption" model may allow computation of the 
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leading edge angle of the VPPA keyhole; while another concept 
appears to be reguired to deal with this problem when the plasma 
flow rate is below the threshold value. 
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X. RESULTS OF COMPUTATION AND COMPARISON WITH 
MEASUREMENTS MADE BY THE GDI TEAM 

A computer implementation of the proposed mathematical model 
for VPPA modeling computation is included in the Appendix A; and 
for VPPA partial penetration to the formation of the keyhole is 
included in the Appendix B. Input data for the computer program is 
provided for in Appendix C for aluminum; and Appendix D for mild 
steel. Nomenclature are listed in Appendix E. 

For the purpose to explain VPPA modeling computation, and VPPA 
partial penetration computation listed in Appendices A and B, two 
block diagrams and corresponding explanations are included in 
Appendices F and G to explain the structures of VPPA partial 
penetration and VPPA modeling computation, respectively. 

Model computations have been made and compared with the 
experimental measurements accomplished by Nichols Research 
Corporation (1992) . These comparisons are included in Appendix H 
of this report. 
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Appendix A 

********* COMPUTER PROGRAM FOR VPF'A MODELING *********** 


c 


PROGRAM WELDBOTH. FOR, COMPILED ON FEB. 26,1993 

Line source , keyhole diameter is variable along the axis 
of keyhole 


D I MENS I ON ROW 1(11), P0W2 (11), P0W3 (11), P0W4 (11) 

DIMENSION ENT 1 (11) , ENT2 (11) , ENT3 (11) , ENT4 (11) , TABMPL (26) 
DIMENSION QLQS1 (10) , QL0S2 (10) , QL0S3 ( 10) , QL0S4 (10) , TAE'MSD (26) 
D I MENS I ON TEMP 1(11), TEMP2 (11), TEMF'3 (11), TEMF'4 (11), BHOLE (10) 
DIMENSION EL0S2 (10) ,EL0S4(10) ,QSEGM(10) , XSOUR(IO) , BHOLEO(IO) 
DIMENSION AREA (10) , DUDL1 (10) , DUDL2 ( 10) , DUDL3 ( 10) , QSOUR (10) 
DIMENSION R 1 (11) , R 2 ( 1 1 ) , FLWPL (11), FLWSD (11), ENERG (11) 

COMMON/ DAT 10 /NX, NY, NZ, X (25) , Y (20) , TM, TO, THK, EPSL , DENS , HFUS I 


COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
OPEN (5, 
OPEN ( 6 , 
OPEN (7, 
READ (5, 
READ (5, 
READ (5, 
READ (5, 
READ (5, 
READ (5, 
READ (5, 
READ (5, 
READ (5, 
READ (5, 
READ (6, 
READ (6, 
READ (6, 
READ (6, 
READ (6, 
READ (6, 


MMM 

TABTEM ( 26 ) , TAB I PL ( 26 ) 
ITMAX, DHMAX , DTMAX 
T AMB I , DNOZL , FMF'L 
TABCSD ( 26 ) , TABKSD ( 26 ) 
TW4, MREG4, NREG4 
F'LM, SDM 

FILE—' C: WELDTAB. USE’ ) 

F I LE= ’ C : WELD I N . USE ’ ) 
FILE—’ C: WELDOU. USE ' ) 

* ) MMM 

(TABTEM ( I ) , 1=1, MMM) 
(TABIPL(I) , 1=1, MMM) 
(TABCPL ( I ) , 1=1 , MMM) 
(TABKF'L ( I ) , 1 = 1 , MMM) 
(TABISD ( I ) , 1=1 , MMM) 
(TABCSD(I) , 1=1, MMM) 
(TABKSD (I) , 1 = 1, MMM) 
(TABMPL ( I ) , 1=1 , MMM) 
( TABMSD ( I ) , 1 = 1 , MMM) 


TABCPL (26) , TABKF'L (26) , TABISD 


? 6 ) 


*) 
*) 
*) 
*) 
*) 
*) 
*) 
* ) 
* ) 
*) 
*) 
*) 
*) 
*) 
*) 


NREG 1 , NREG2 , NREG3 , NREG4 
VA, VC 

FIA1, FIA2, FI A3 
TIME, TIMER 

D I SL 1 , D I SL2 , D I SL3, D I SL4 
VVE , WK 


READ (6, *) 
READ (6, *) 
READ (6, *) 
READ (6, *) 
READ (6, *) 
READ (6, *) 
READ (6, *) 
READ (6, *) 
READ (6, *) 
READ (6, *) 
READ (6, *) 


APIL, APILR, AMA I , AMAIF: 
DNOZL , DSHID, VF'L, VSD 
F'LMOL, SDMOL, UF) 

ITMAX, DHMAX, DTMAX 

TW1 , TW2 , TW4 , HSUR2, HSUR4 

DXW, TAMBI , DF'OLAR, AMDW 

F'S I A , F'AMB I , CF'SD 

ALFA, VWELD 

NX, NY 

X (1) , X (NX) 

Y ( 1 ) , Y (NY) 


F'age 


1 



READ (6,*> TM , TO 

READ (6,*) THK, EF'SL 

READ (6,*> DENS, HFUSI 

VQLF'L— VF'L/ 1 . 27133/1 . E5 

, v'QL3D = VSD/ 1.271 33/ 1 . E5 

FMPL=V0LF'L*PAMBI*PLM0L/8315. /TAMBI 

FMSD=U0LSD*PAMBI*SDM0L/8315. /TAMBI 

UH=V0LSD*4. /3. 14159/ (DSHID**2-DN0ZL**2) 

AF:EAW=3. 14159*DF'0LAR**2/4. 

WRITE (7, 1030) D I SL 1 , D I SL2 , D I SL3 , D I SL4 

WRITE (7, 1040) APIL, APILR, AMAI , AMAIR 

WRITE (7, 1060) FIA1 , FIA2, FIA3 

WRITE (7, 1070) VA, VC 

WRITE (7, 1130) F'LMOL , SDMOL , UR 

ELOS= 1 . 5*WK/WE/ 1 . E10 

GHW= ( (AMAI+AF'IL) * (VC-FIA1 ) *TIME+ ( AMAI R- APILR) * (VA+FIA1 ) * 

$ T I MER ) / ( T I ME+T I MER ) / 1 000 . 

CONDW=AMDW*AREAW* (TW1-TAMBI) /DXW 
P0W1 ( 1 ) =GHW— CONDW 
ENT 1 ( 1 ) =P0W1 ( 1 ) /FMF’L 

CALL INSEF' (TAB I PL, TABTEM, ENT 1 (1) , TEMPI (1) ) 

ELQS1=EL0S*TEMP1 (1 ) * ( AMAIR-APILR ) *TIMER/ (TIME+TIMER) 
AMF'CP1= ( (AMAI+AF'IL) *TIME+ (AMAIR-APILR) *TIMER) / (TIME+TIMER) 
DX1=DISL1/FLQAT (NREG1 ) 

GE0L1=3. 14159* ( ( DF'OLAR+DNOZL ) /2. ) **2/4. /DX1 
V0LUM1=3. 14159* ( (DNOZL+DPOLAR) /2. ) **2/4. *DX1 
MREG 1 =NF:EG 1 + 1 
SUMQ 1 =0 . 

DUBC=0. 

DO 40 M=2 , MREG 1 
N=M— 1 

DUDL1 (N) =2000. 

OLDS 1 ( N ) =0 . 

NN1 =0 

10 CONTINUE 
NN 1 =NN 1+1 

I F ( NN 1 . GT . 200 ) THEN 

WRITE (*,*)’ NN1.GT.200. Err a stop' 

STOP 
END IF 

DP1=DUDL1 (N) *DX 1 *AMF‘CF'l / 1000. — QL0S1 (N) 

POW 1 ( M ) =POW 1 ( N ) +DP 1 
ENT 1 (M) =P0W1 (M) /FMF'L 
ENTCF'1= (ENT 1 (N) +ENT1 ( M ) ) / 2 . 

CALL I NSEP ( TAB I PL , TABTEM , ENTCF' 1 , TEMCP 1 ) 

QL0S1 ( N ) = VOLUM 1*4. *1 . E7* (TEMCP 1 / 15000. ) **16/ 

$ ( 1 . + ( TEMCP 1/1 5000 . ) **16) 

CALL I NSEP (TABTEM, TABMPL, TEMCP 1 , F'LMHO) 

QHM1=1 . /PLMH0/GE0L1 
DUDL=AMF'CP1*0HM1/DX1 
DIFU=1 . — DUDL/DUDL1 (N) 

IF (AB3 (DIFU) . GE. 0.005) THEN 
DUDL1 (N> = (DUDL+DUDL1 (N) ) /2. 

GO TO 10 


F'age 



40 


C 

c 

c 

c 

c 

c 


END I F 

DUBC=DUBC+DUDL1 (N) *DX1 
SLIMQ 1 =3IJMQ 1 +QL0S 1 (N) 

ChLL I MSEF' ( TAB I PL , T ABTEM , ENT 1 < M ) 
CONTINUE 

( ROW 1(1), 

( ENT 1 ( I } , 

(0L0S1 ( I ) 

(TEMPI ( I ) 

(DUDL1 (I) 


TEMPI (M) ) 


1 170) 
1180) 
1 190) 
1200 ) 
1 220 ) 
1210 ) 


SUMQ1 , DUBC 


1=1 , MREG1 ) 
1=1 , MPEG 1) 

, 1 = 1 , NREG1 ) 
, 1 = 1, MREG1 ) 
NREG1 ) 


1 = 1 


WRITEC 
WRITEC 
WRITEC 
WR I TE ( 7 ’ 

WRITE (7, 

WRITE (7 
C 

P0W2 ( 1 ) =F'0W 1 (MPEG 1 > 

ENT2 ( 1 ) =ENT 1 (MREG1 ) 

TEMP2 ( 1 > = TEMP 1 ( MREG 1 ) 

DX2«d7sl2/"SIt ( WEG 2? IR " flPILR> * TI " ER ’ ' <TIME+TIMER > 

AREA2=3. 14159*DN0ZL*DX2 
GE0L2=3 . 14159*DN0ZL##2/4. /DX2 
MREG2=NF:EG2+ 1 
2 WRITE (7, 1230) 

GHN1=<AF'IL*FIA2*TIME+APILR*FIA2*TIMER) / (TIME+TIMER) /lOOO 
DUCG=0. 

DO 80 M=2 , MREG2 
N=M— 1 


50 


$ 


XM=FLOAT (N) *DX2-DX2/2. 

DUDL2 (N) =2000. 

QL0S2 (N) =0. 

NN2=0 
CONTINUE 
NN2=NN2+ 1 

I F ( NN2 . GT . 200 ) THEN 

WRITEC*,#)-’ Erro Stop ENT3=-’ ENT3 

STOP ' - 

END IF 


IF (NREG2. EQ. 1) THEN 


- N> - ELas * TEMp 2 (N) * <APIL*TIME+APILR*TIMER) / (TIME+TIMER) 
ELSE IF (NREG2. NE. 1 . AND. N. EQ. 1 ) THEN 
EL0S2 (N) =EL0S#TEMF'2 (N) #AF'IL*TIME/ (TIME+TIMER) 

ELSE IF (NREG2. NE. 1 . AND. N. EQ. NREG2) THEN 

= EL0S#TEMP2 (N) *AF'ILR#TIMER/ (TIME+TIMER) 

END I F 

DF'2=DUDL2^N^#DX2#AMPCP2/ 1000. + ( 1 . -PS I A) *GHN2/FL0AT (NREG2) 


P0W2 ( M ) =F'0W2 ( N ) +DP2 
ENT2 (M) =F'0W2 (M) /FMF'L 
ENTCF- (ENT2 (N) +ENT2 (M) ) /2. 

CALL I NSEP ( TAB I PL , TABTEM , ENTCP , TEMCP ) 
CALL I NSEP (TABTEM, TABKPL, TEMCP, AMD) 
CALL I NSEP (TABTEM, TABCPL, TEMCP, CF') 
CALL I NSEP ( T ABTEM , T ABMF'L , TEMCP F'LMHO ) 
C0EFC=2. # FMF'L #CP/3. 14159/AMD 


Page 


80 


C 

C 

C 

C 

C 

C 

C 

C 

c- 


90 


C 


XPLUS=C0EFC/XM 

VNUX=0 . 2233*XF’LUS**0. 7455+3.66 
ALFA X = VMU X * AMD/ DNOZL 

QL0S2 (N> =ALFAX*AREA2* (ENTCP-HSUR2) /CP 
0HM2=1. / PLMHQ / GE0L2 
DUDL=AMPCP2*0HM2/DX2 
DIFU=1 . -DUDL/DUDL2 (N) 

IF (ABS (DIFU) . GE.O. 005) THEN 
DUDL2 (N) = (DUDL+DUDL2 (N) )/2. 

GO TO 50 
END I F 


DUCG=DUCG+DUDL2 (N) *DX2 
SUMQ2=SUMQ2+QL0S2 (N> +EL0S2 (N) 

CALL INSEP < TAB I PL, TABTEM, ENT2 (M) , TEMP2 (M) ) 
CONTINUE 


WRITE (7, 1250) 
WRITE (7, 1260) 
WRITE <7, 1270) 


WRITE (7, 1280) 
WRITE (7, 12<k>) 


WRITE (7, 1300) 
WRITE (7, 1310) 
WRITE (7, 1320) 


(P0W2 ( I ) , 1=1 , MREG2) 

( ENT2 ( I ) , 1=1 , MREG2) 
(QLQS2 < I ) , 1 = 1 , NREG2) 
(TEMP2 Cl) , 1=1, MREG2) 
( EL0S2 ( I ) , 1=1 , NREG2) 
(DUDL2 ( I ) , 1=1 , NREG2) 
SUMQ2 , DUCG 


R0=DN0ZL/2. 

P0W3 ( 1 ) =P0W2 ( MREG2 ) 

MITER=0 

CONTINUE 

MITER=MITER+1 
IF (MITER. GT. 200) THEN 

WRITE'*,*)’ Erro Stop MITER. GT. 200, DPP= ! 
STOP 
END IF 

ENT3 ( 1 ) =P0W3 < 1 ) /FMPL 

CALL I NSEP ( TAB I PL , TABTEM , ENT3 ( 1 ) , TEMP3 ( 1 ) ) 
DENS0=PAMBI*PLM0L/S315. /TEMP3 ( 1 ) 

UO=FMPL / DENSO / 3 . 1 4 1 59 / RO * * 2 
ENERG ( 1 ) =FMPL*U0**2/2. /1000. 

PPP=P0W2 (MREG2) -ENERG ( 1 ) 

DPP=ABS ( 1 . — PPP/P0W3 ( 1 ) ) 

P0W3 ( 1 ) =PPP 

IF (DPP. GE.O. 005) GO TO 90 
COEFM=UH / UO 

WRITE (7, 1340) UH,UO,COEFM 
VMID1= ( 1 . — COEFM) / (l.+COEFM) 

0MID3=1 . +0. 8*C0EFM-0. 45*C0EFM**2 
DX3=DISL3/FL0AT (NREG3) 

MREG3=NREG3+ 1 
SUMQ3=0 . 

DUGE=0. 

R1 ( 1 ) =R0 
R2 ( 1 ) =R0 
FLWPL ( 1 ) =FMPL 
FLWSD ( 1 ) =0 . 


, DPP 
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T I ME + TIMER) 


120 


150 


AMPCP3= ( AM A I * T I ME + AM A I R# T I MEF: ) / 

DO 190 M=2 , MREG3 
N=M— 1 

X3=FL0AT (N) *DX3 
VMID2=0. 27* X3 

VM I D4=VM I D 1 * ( 0 . 4 1 6+0 . 1 34*C0EFM+0. 02 1 * VMI D2* VM ID1*VM ID3/RO ) 

R 1 ( M ) =F:0-VM I D2* VM I D4 

VMIDA=0. 021 *VMID3/R0 

VMIDB=0. 416+0. 134*C0EFM 

VM I DC=R 1 (M) -RO 

VM I DD= ( - VM I DB+SQRT (VMIDB**2-4. *VMIDA*VMIDC) ) /2 /VMIDA 

R2 (M) =R1 (M) +VMIDD 

N0M=0 

TBEG I =TEMP3 ( 1 ) 

DENBEG=DENSO 

DR= (R2 (M) — R1 (M) ) /10. 

DUDL2 < N ) =2000. 

GE0L3=3. 14159*R2 (M) **2/0X3 

CONTINUE 

NOM=NOM+ 1 

IF (NOM. GT. 200) THEN 

WRITE (*,*.)’ Erro Stop ITER. IN STANDOFF EXCEED 200’ 

i LJ i 

END IF 

FLWF'L (M) =3. 14159*R1 (M) **2*UO*DENBEG 
FLWSD (M) =0. 

ENERG (M) =FLWPL (M) *U0**2/2. / 1000. 

CALL INSEF' (TABTEM, TABIF'L, TBEG I , EBEGI ) 

POWER=FLWPL (M) *EBEGI 
DO 150 K=1 , 10 

RMID=FLOAT <K) *DR-DR/2. +R1 (M) 

ATA= (R2 (M) -RM ID) / (R2 (M) — R1 ( M ) ) 

UMID=UO- OJO-UH) * ( 1 . -ATA* * 1 . 5) **2 
TMID=TBEGI— (TBEGI-TAMBI ) * ( 1 . -ATA) 

F'GPL=PAMBI / ( 1 . + ( 1 . /ATA-1 . ) *PLMOL/SDMOL) 

PGSD=PAMB I -PGPL 
DENF'L=PGPL*PLM0L/8315. /TMID 
DENSD=PGSD*SDM0L/83 15. /TMID 
V0LUM=2. *3. 14159*RMID*DR*UMID 
FLWF'L (M) =FLWPL (M) +VOLUM*DENF'L 
FLWSD (M) =FLWSD (M) +VOLUM*DENSD 

ENERG (M) =ENERG (M) + ( DENF'L+DENSD ) *V0LUM*UMID**2/2. /lOOO 

IF (TMID. GE. 3000. ) THEN 

CALL INSEF' (TABTEM, TABIF'L, TMID, ENTPL) 

CALL INSEF' (TABTEM, TABISD, TMID, ENTSD) 

ELSE 

ENTF'L=T ABCF'L ( 1 ) *TMID 
ENTSD=T ABCSD ( 1 ) *TMID 
END I F 

F'OWER-POWER+VOLUM* ( DENF'L*ENTF'L+DENSD*ENTSD) 

CONTINUE 

POWF:CF'= ( P0WER+P0W3 ( N ) ) /2. 

PLCF - (FLWF'L (M) +FLWPL (N) ) /2. 

SDCF'= (FLWSD (M) +FLWSD (N) ) /2. 
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CALL AVERT (F'LCP, SDCP, POWRCP , TABTEM, TAB I PL, TAB ISO, TCP) 

VRADI-3. 14159* (R2(M) +R2 (N) ) **2/4. *DX3 

QL0S3 (N ) =VRAD 1*4. *1 . E7* (TCP/ 1 5000. ) ** 16 / ( 1 . + (TCP/ 1 5000. ) **16) 
CALL INSEP (TABTEM, TABMPL, TCP, PLMHO) 

CALL INSEP (TABTEM, TABMSD , TCP, SDMHO) 

CPMHO= ( (FLWPL (M) +FLWPL (N) ) *PLMHQ+ (FLWSD (M) +FLWSD (N) ) *SDMHO) / 
$ (FLWPL (M) +FLWPL (N) +FLWSD (M) +FLW3D (N) ) 

0HM3= 1 . /CPMH0/GE0L3 
DUBL3 (N) =QHM3* AMPCP3/DX3 
DP3=0HM3* AMPCP3* *2/ 1 000 . 

P0W3 (M) =F'QW3 (N) +DP3-QLDS3 (N) + (FLWSD (M) -FLWSD (N) ) *CPSD*TAMBI- 
$ ( ENERG ( M ) -ENERG ( N ) ) 

PMF‘= 1 . — POWER /PQW3 (M) 

IF (ABS (PMP) . GE. 0.010) THEN 
TBEGI=TBEGI * ( 1 . +PMP/10. ) 

DENBEG=PAMBI*PLM0L/8315. /TBEGI 
UO=UO*FMPL/FLWPL (M) 

GO TO 120 
END IF 

WRITE (#, 1390) N, NOM, PMP 
SUMQ3=SUMQ3+QL0S3 (N) 

CALL AVERT (FLWPL (M) , FLWSD (M) , POWER, TABTEM, TABIPL, TABI3D, 

* TEMP3 (M) ) 

ENT3 (M) =F'0W3 (M) / (FLWPL (M) +FLWSD (M) ) 

DUGE=DUGE+DUDL3 (N) *DX3 
190 CONTINUE 

WRITE (*, 1480) (FLWPL (I) , I=1,MREG3) 

WRITE (7, 1480) (FLWPL(I) , I=1,MREG3) 

WRITE (7, 1490) (FLWSD ( I ) , 1 = 1 , MREG3) 

PAID=FLWSD (MREG3) *CPSD*TAMBI 


RBOUND=RO 

ITWID=1 

DO 191 N= 1 , NREG4 
BHOLE (N) =RBOUND 
191 CONTINUE 
195 CONTINUE 

ATARB= (F:2 (MREG3) -RBQUND) / (R2 (MREG3) -R1 (MREG3) ) 
TBOUND=TBEG I - ( 1 . -ATARB) * ( TBEG I-TAMB I ) 

WRITE (7, 1540) TBOUND 
DR= ( R2 ( MREG3 ) -RBOUND ) / 10. 

F‘LL0S=0. 

SDL0S=0. 

QL0SJ=0. 

QGETJ=0. 0 
SDGET=0. 0 
DO 220 K= 1,10 

RMID=FLOAT (K) *DR-DR/2. +RBOUND 
ATA= (R2 (MREG3) -RMID) / (R2 (MREG3) -R1 (MREG3) ) 
UMID=U0- (UO-UH) * (1. — ATA** 1.5) **2 
TMID=TBEGI- (TBEG I-TAMB I ) * ( 1 . -ATA) 

PGPL=F'AMBI / ( 1 . + ( 1 . /ATA-1 . ) *PLMOL/SDMOL) 
F'GSD=PAMB I — F'GPL 
DENPL=PGPL*PLM0L/83 1 5 . /TMID 
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DENSD=PGSD SSDMOL/83 15. / TM I D 
V0LUM=2 . * 3 . 14159*RMID#DR*(JMID 
P L L 0 B = P L L 0 S + V 0 L U M * D E N F' L 
SDLuS=3DL0S+V0LUM>KDENSD 
I F ( TM I D . GE . 3000 . ) THEN 

CALL INSEP (TABTEM, TAB I PL, TMID, ENTPL ) 
CALL INSEP (TABTEM, TABISD, TMID, ENTSD) 
ELSE 

ENTPL=T ABCPL ( 1 ) *TMID 
ENTSD=TABCSD (1 ) *TMID 
END IF 


QLOSJ=QLOSJ+VOLUM# ( DENF'L*ENTF’L+DENSD*ENTSD) 
220 CONTINUE 


C 


c — 


GHK 1 = 
GHK2= 


( AM A I *FI A3*TIME+AMAIR* (UR-FIA3) *TIMER) / (TIME+TIMER 
(AMAI *VA*TIME+AMAIF:*VC#TIMER) / (TIME+TIMER) /1000. 


) / 1000 . 


225 CONTINUE 


P0WER=P0W3 (MREG3) -QLOSJ+OGETJ 
PLM=FMFL— PLLOS 

SDM=FLWSD (MREG3) -SDLOS+SDGET 
P0W4 ( 1 ) =F'OWER+ ( 1 . -PS I A) *GHK2 
ENT 4(1) =P0W4 ( 1 ) / ( PLM+SDM ) 

CALL AVERT (PLM, SDM, P0W4 ( 1 ) , TABTEM, TAB I PL, TABISD, TEMF'4 < 1 ) ) 
SUMQ4=0 . 

DP4=0. 


EL0SS=EL0S*TEMF'4 ( 1 ) Xc AMAI *TI ME/ (TIME+TIMER) 
GHK12=GHK1+F'SIA*GHK2 
C WRITE • 7 , 1555) 

DL=DISL4/FL0AT (NREG4) 

MREG4=NREG4+1 
DO 230 1=2, MREG4 
J=I-1 

XSOUR ( J ) =FLOAT ( J ) *DL-DL/2. 

AREA (J) =3. 141 59*BH0LE ( J ) *DL 
230 CONTINUE 

C WRITE (7, 1560) ( XSOUR <K> , K»1 , NREG4) 

C WF;ITE (7, 1570) RBQUND 

C WRITE (7, 1580) (AREA (K) , K=1 , NREG4) 

WRITE (7, 1590) PLM, SDM, TEMF'4 ( 1 ) 

DO 240 M=2,MREG4 
N=M— 1 


% 


OSEGM (N) — JKGHK12# ( 1 . -XSOUR (N) /DISL4) XcXc2XcDL/DISL4 
EL0S4 (N) =3. *ELOSS* ( 1 . -XSOUR (N) /DISL4) **2*DL/DISL4 
CmLL C0NV4 ( XSOUR (N) , AREA (N) , DF'4, F'0W4 (N) , P0W4 (M) , ENT4 (M) , 

QL0S4 (N) , HSUR4, EL0S4 (N) , BHOLE (N) ) 

CALL AVERT (PLM , SDM, P0W4 (M) , TABTEM, T AB I F'L, TABISD, TEMF'4 (M) ) 

QSOUR (N) -0L0S4 (N) +EL0S4 (N) +QSEGM (N) 

SUMQ4=SUMQ4+QS0UF.' < N ) 

240 CONTINUE 

ITWID=ITWID+1 MQ4 '' Q3aUR '’ XS0UR ’ NREG4, DISL4, ALFA, VWELD , BHOLEO , ITWID) 
" R0=? ? RBQUND*' , R'BQUND, ’ R2(MREG3)=’ ,F:2(MF:EG3) 

DKo-U . U 
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c 

c 

c 

c 

c 

c 

c 

c 


DO 250 J= 1 , NREG4 

DRB=AMAX1 ( DRB , ABS ( BHOLEO ( J ) -BHOLE ( J ) ) /BHOLE < J ) ) 
250 CONTINUE 

I F ( DRB . GE . 0 . 002 ) THEN 
DO 260 J = 1 , NREG4 

BHOLE ( J ) =BHOLE ( J) +. 5* (BHOLEO (J) -BHOLE ( J) ) 

260 CONTINUE 

RBOUND= BHOLE ( 1 ) 

IF (RBOUND. GE. R2 (MREG3) ) THEN 
DENSD=F'AMB I *SDM0L/83 15./ TAMB I 
QLQSJ =0 . 0 


PLLQS=0. 0 
SDL0S=0. 0 

SDGET=3. 141 59* (RBOUND* *2-R2 (MREG3) **2) *UH*DENSD 
QGETJ=SDGET*TAIiBI*CPSD 
GO TO 225 
ELSE 

GO TO 195 
END I F 
END IF 

WRITE (7, 1600) 


WRITE (7, 1610) 
WRITE (7, 1620) 
WRITE (7 , 1630) 
WRITE (7, 1640) 
WRITE (7, 1650) 
WRITE (7, 1660) 
WRITE (7, 1670) 
WRITE (7, 1680) 


□IMP IN 

( P0W4 ( I ) , 1=1 , 
( ENT 4(1) , 1=1, 
(TEMF’4 ( I ) , 1 = 1 
( QL0S4 ( I ) , 1=1 
( EL0S4 ( I ) , 1=1 
(QSEGM(I) ,1=1 
(QSOUR(I) , 1=1 
SUMQ4 


MREG4) 
MREG4) 
MREG4) 
NREG4) 
NREG4) 
NREG4) 
NREG4 ) 


PT0T0U=CQNDW+SUMQ1+SUMQ2+SUMQ3+SUMQ4+P0W4 (MREG4) 

$ +QLOSJ+ENERG (MREG3) 

PTOT I N=GHW+ ( AMF'CP 1 *DUBC+AMF'CP2*DUCG+AMF'CF'3*DUGE > / 1 000 . 
* +GHN 1 +GHN2+GHK 1 +GHK2+PA I D 

WRITE (7, 1690) PTOT I N , PTQTOU 
F'ERCQW=CONDW/PTOT IN 
PERQG=SUMQ 1 /F'TOT I N 
F'ERQN=SUMQ2/PT0TIN 
F'EF:QSR=SUMQ3 /F'TOT IN 
F‘ERQSC=QLQSJ / PTOT I N 
PERQW=SUMQ4/PT0TIN 
PERQD=F'0W4 (MREG4) / F'TOT IN 

F'ERSUM=PERCQW+PERQG+PERQN+PERQSR+PERQSC+F'ERQW+PERQD 
WR ITE (7, 1 700 ) CONDW , F'ERCQW 


SUMQ1 

SUMQ2 

SUMQ3 

QLOSJ 

SUMQ4 


WRITE (7, 1710) 

WRITE (7, 1720) 

WRITE (7, 1730) 

WRITE (7, 1740) 

WRITE (7, 1750) 

WRITE (7, 1760) 

WRITE (7, 1770) 

WRITE (7, 1780) 
DUBE=DUBC+DUCG+DUGE 
DUAF=DUBE+VA+VC+FIA3— FIA1 
DUAFR=DUBE+VA+VC— F I A3+F I A 1 +UR 


F'ERQG 
PERQN 
PERQSR 
PERQSC 
PE ROW 
P0W4 (MREG4) , F'ERQD 
F'ERSUM 

AMF'CP 1 , AMF'CF'2, AMF'CF'3 


Page 


8 



WR i TE •; 7 , 17 90 ) DUBE , DUAF , DU AFF; 

DUF'F=DUBC+VA+VC+F I A2-F I A 1 
DUF'F.-F I A2 -F I A 1 -DUBC-DUCG 

F'TOT AL= ( ( AM A I * DU AF+ AF' I L # DUF'F ) * T I ME+ < AM A I RKDUAFR+AP I LR * DURR ) # 
$ TIMER) / (TIME+TIMER) / 1000. 

C F'TOT AL = ( AM A I KDUAF*. T IME + AMA I RfcDUAFRstcT IMER) / (TIME+TIMER) /lOOO. 

ATAQ4=SUMQ4/F'T0TAL 
WRITE (7, 1800) DUF'F , DUF'R 
WRITE (7, 1810) F'TOT AL , SUMQ4 , AT AQ4 
WR I TE ( 7 , # ) ’ CROWN WIDTH®’ , 2. #BHGLE ( 1 ) , ’ (mm) ’ 

WRITE (7, *>’ ROOT WIDTH®’ , 2. *BHOLE (NREG4) , ’ (mm) ’ 

STOP 

C 

1000 FORMAT <///5X, ’ INPUT DATA:’) 

1010 FORMAT (/2X, ’ F'LSMA (kg/s) ’, E12. 5, 2X, ’ SHIELD (kg/s) ’ , E12. 5) 
1020 FORMAT ( /5X, ’ TW OF W, NOZZLE, HOLE, INITIAL SOLID’ , 3F8. 1 , / 

* 5X , ’ ENT . OF ARC IN TW OF NOZZLE, HOLE’ , 2E12. 5) 

1030 FORMAT </5X, ’Ll, L2,L3,L4 (m> ’ , 2X, 4F10.8) 

1040 FORMAT (/5X, ’PILOT I (+)=’, F8. 2, 5X, ’ PILOT I (-) =’ , FS. 2, / 

* 5X, ’ MAIN I ( + ) =’ , F8. 2, 5X, ’ MAIN I (-)=’, F8. 2) 

1050 FORMAT (/5X, ’t ( + )=’ ,E11.4,5X, ’t (-)=’ , El 1.4) 

1060 FORMAT </5X, ’WORK FUNCTION 1,2,3 ®’,3F5.2) 

1070 FORMAT (/5X, ’ ANODE DR0F'=’,F5.2^ 

* 2X , ’ CATHOD DROP®’ ,F5.2) 

1080 FORMAT (/5X, ’ e= ’ , E 1 2. 5, 5X , ’ k = ’ , E12. 5) 

1090 FORMAT (/5X, ’ NREG 1 = ’ , 13, 2X, ’ NREG2=’ , 13, 2X, ’ NREG3® ’ , 13, 

* 2X, ’NF:EG4=’ ,13) 

1100 FORMAT (/5X, ’NOZZLE DI AM. (m> ’ , F8. 6, 2X, ’ SHIELD DI AM. (m) ’ , F8. 6) 
1110 FORMAT (/5X, ’LENGTH OF POLARITY W =’,E12.5, 

* /5X, ’ ENVIROMENT TEMP. =’ , F6. 1 , 2X , ’ AREA OF W POLARITY’, 

$ E12.5, /5X, ’AMD OF POLARITY W =’,E12.5) 

1 120 FORMAT ( /5X , ’ ITMAX , DHMAX , DTMAX=’ , 14, El 2. 5, F6. 4) 

1125 FORMAT ( /5X , ’F'SIA=’ ,F6.4,2X, ’ F'AMB I (N/m**2) =’ , F10. 1, / 

$ 2X,’CF'. OF HE IN AMBIENT ( kJ/kg ) =’ , E12. 5, 

$ 2X , ’ (m/s) ’ , E12. 5) 

1130 FORMAT (/5X, ’ PLASMA MOL. WEIGHT ( kg/kmol ) =’ , F8. 4, 2X, 

$ /5X, ’SHIELD MOL. WEIGHT ( kg/kmol )=’, FS. 4, 

$ /5X, ’UR®’ ,F5. 1) 

1140 FORMAT (/2X, ’ GHW , CONDW, ELOS 1 (kw> ’ ,3E12.5) 

1150 FORMAT ( //5X, ’ BEGINING OF SEGMENT ONE’) 

1170 FORMAT (/2X, ’F'OWl®’ , 4E12.5, /7X,4E12.5) 

1180 FORMAT (/2X, ’ENT1®’ , 4E12. 5, /7X , 4E12. 5) 

1190 FORMAT (/IX, ’QL0S1 = ’ ,4E12.5, /7X, 4E12.5) 

1200 FORMAT (/2X, ’ TEMPI®’ ,4E12.5, /8X,4E12.5) 

1210 F0RMAT(/2X, ’ TOTAL LOSS IN GAP I ’ , E12. 5, 2X, ’ DUBC (v) ’ , F6. 2) 

1220 FORMAT ( /2X , ’ DUDL 1 ’ ,4E12.5, /7X.4E12.5) 

1230 FORMAT (//5X, ’NOZZLE COMPUTATION’) 

1240 FORMAT </5X, ’NN2.GT. 100’ ) 

1250 FORMAT (/5X, ’F'0W2=’ , 4E 1 2 . 5/ 10X , 4E 1 2 . 5) 

1260 F0RMAT(/5X, ’ ENT2=’ , 4E 1 2. 5/ 1 OX , 4E1 2. 5) 

1270 FORMAT (/4X, ’QL0S2=’ , 4E12. 5/ 10X . 4E12. 5) 

1280 FORMAT (/2X, ’ TEMF'2 ’ ,4E12.5, /7X,4E12.5) 

1290 F0RMAT(/2X, ’ EL0S2 ’ , 4E12. 5, /7X , 4E12. 5) 

1300 FORMAT </2X, ’ DUDL2 ’ , 4E 1 2 . 5 , /7X , 4E1 2. 5) 
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1310 FORMAT 
1320 FORMAT 
1340 FORMAT 
1350 FORMAT 
1370 FORMAT 
Cl 360 FORMAT 
1380 FORMAT 
1390 FORMAT 
$ 

1450 FORMAT 
1410 FORMAT 
1430 FORMAT 
1440 FORMAT 
1460 FORMAT 
1470 FORMAT 
1480 FORMAT 
1490 FORMAT 
1500 FORMAT 
1510 FORMAT 
1520 FORMAT 
1530 FORMAT 
1540 FORMAT 
1550 FORMAT 
1555 FORMAT 
1560 FORMAT 
1570 FORMAT 
1580 FORMAT 
1590 FORMAT 
1600 FORMAT 
1610 FORMAT 
1620 FORMAT 
1630 FORMAT 
1640 FORMAT 
1650 FORMAT 
1660 FORMAT 
1670 FORMAT 
1680 FORMAT 
1690 FORMAT 
$ 

1700 FORMAT 
1710 FORMAT 
1720 FORMAT 
1730 FORMAT 
1740 FORMAT 
1750 FORMAT 
1760 FORMAT 
1770 FORMAT 
1780 FORMAT 
1790 FORMAT 
1800 FORMAT 
1810 FORMAT 
END 

C 

C 


TOTAL LOSS IN NOZZLE = ’ , E 1 2. 5, 2X , ’ DUCG ( 
.•’STANDOFF COMPUTATION’) 

' UH, UO , COEFM ’ , 3E 1 2 . 5 > 

’ LINE 285 TBEGI ’ , El 2. 4) 

’ AT A , TBEG I , TAMB I , TM I D ’ , 4E 1 1 . 4 ) 

M , R 1 , R 2 , D R , R M I D , A T A ’ , 1 2 , 2 E 1 2 . 5 , / 5 X , 3 E 1 2 . 5 ) 
’IN MAIN, NOM , MAR , MHE , POWER ’ , 14, 2X , 3E10. 3) 

’ ITER. IN STANDOFF, SEG. =’ , 12, 2X, ’NOM=’ ,14, 

’ PMF'=’ , F10. 6) 

’TOTAL LOSS IN STANDOFF = ’ , E 1 2 . 5 , 2 X , ’ DUGE ( v ) 
’ P0W3=’ , 4E12.5, /10X, 4E12. 5) 

’ QL0S3= ’ , 4E 12.5, /10X, 4E12. 5) 

4E 1 2 . 5 , / 1 OX , 4E 1 2 . 5 ) 

4F10.8, / 1 1 X , 4F 10.8) 

4F10. 8, /I IX, 4F10.S) 

4E12.5, /11X,4E12.5) 

4E12.5, / 1 1 X , 4E1 2. 5) 

4E12.5, / 1 1 X , 4E 12. 5) 

4E12. 5, /7X, 4E12.5) 


( / 5 X , 

( / / 5 X 
( / 5 X , 

( / 2 X , 

( /IX, 

ax,-’ 

( / 2 X , 

(/2X, 

2X , 

(/2X, 

( /5X , 

( / 4 X , 

( / 4 X , ’ TEMF'3= ■’ 
</5X, ’ R 1 < m ) = ’ 
(/5X, ’ R2 <m) =’ 
C/5X, •’ FLWF'L= ’ 
( /5X , ’FLWSD=’ 
( / 5 X , ’ ENERG= ’ 
(/2X, •’ DUDL3’ , 


) ’ ,F6. 


, F6. 2) 


</5X, ’POWER AID FROM SHIELD, PAID <kW) = 7 , E12. 5) 
( /2X , ’ VELOC. , TEMP. IN CORE AT LAST’ , 2E12. 5) 
(/5X,’T OF BOUNDARY, TBOUND=’ , E12. 5) 

(/5X, ’CONVECTIVE LOSS OF JET <kW)’,E12.5> 
<//10X, ’WORKPIECE COMPUTATION’) 

(/2X, J XSOUR ’ , 2X , 5F8 . 6 ) 

RADIUS OF KEYHOLE (m) =’ , E12. 5) 


(/5X, 
(/2X, 
(/4X, 
( /5X , 
(/2X, 
(/2X, 
( / 2 X , 
( /2X , 


’ AREA (m**2) ’ , 4E12. 
’ F'LM , SDM , TEMF'4 ( 1 ) •’ 
’IMPINGE LOSS (lew) 
’ F'0W4= ’ , 4E12. 5, /7X 
4E12. 5, /7X 
, 4E12. 5 
, 4E1 


ENT 4= ’ , 
TEMF'4= ’ 
QL0S4=’ 
( / 2 X , ’ EL0S4= ’ 
( / 2 X , ’ QSEGM= ’ 
( / 2 X ’ QSOUF:=’ 


. U 

, 4E12. 5 
i 4E12. 5 
, 4E12. 5 


5X, ’TOTAL LOS 


i, /12X, 4E12.5) 
3E14.5) 

, F8. 5) 

4E12. 5) 

. 4E12.5) 

/3X, 4E12. 5) 

/8X, 4E12.5) 

/ 8 X , 4E12.5) 

/8X, 4E12.5) 

/8X , 4E 12.5) 


S IN WORKPIECE =’,E12.5) 


, E12. 


i / 

(5X, ’TOTAL POWER INCLUDING HELIUM, PTOTIN=’ 

5X , ’ TOTAL POWER FROM LOSS, F'TOTOU=’ , E12. 5) 

(/2X, ’POLARITY CONDUC. LOSS AND PERCENT . ” , 2E 12.5 
(/2X,’GAP LOSS AND F'EF:CENTAGE’ , 2E12. 5) 

(/2X, ’NOZZLE LOSS AND PERCENTAGE’ , 2E12. 5) 

(/2X, ’RAD I AT. IN STANDOFF AND PERCENTAGE’ , 2E 12. 5) 

( /2X , ’ CONVEC. IN STANDOFF AND PERCENTAGE’ ^ 2E1^ 5) 
(/2X, ’WORKPIECE LOSS AND PERCENTAGE’ , 2E 12. 5) *" ’ 

< /2X ,’ DISCHARGE LOSS AND PERCENTAGE’ , 2E12. 5) 

( / 5 X , ’ SUM . PERCENTAGE’ ,F8. 5) 

C/5X, ’ AMF'CP 1,2,3’ ,3E12. 5) 

(/5X, ’ DUBE, DUAF, DUAFR’ , 3E12.5) 

(/4X, ’ DUPF= ’ , E12. 4, 4X , ’DUPR=’ ,E12.4) 

</2X, ’ POWER, Q4,ATA04=’ , 3E14.5) 
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n n n 


c 


200 


SUBROUT I NE C0NV4 ( XM I D , AM I D , DP , F' A , F'B , HB , F'CONV , HW , OELOS , RHOLE ' 
COMMON MMM ' , n,nui_c. . 

COMMON TABTEM (26) , TAB IF'L (26) , TABCPL (26) , TABKF'L (26) , TAB I SD <^6' 
COMMON I TMAX , DHMAX , DTMAX ’ 

COMMON T AMBI , DNQZL , FMF'L 
COMMON T ABCSD ( 26 ) , TABKSD ( 26 ) 

COMMON TW4, MREG4, NREG4 
COMMON PLM , 3DM 
PB=F'A 
NQM=0 
CONTINUE 
NOM=NOM+ 1 


IF ( NOM . GT . I TMAX ) THEN 
WRITE (#,#) ’ Erro Stop 
STOP 
END IF 


ITERATIONS EXCEED I TMAX IN C0NV4 : ' 


PCP= (F'A+F'B) /2 . 

HCF'=PCP/ (F'LM+SDM) 

CALL AVERT (PLM, SDM, F'CP, TABTEM, TAB I PL, TABISD, TCP) 
CALL I NSEF' (TABTEM, TABKF’L, TCP, F'LAMD) 

CALL I NSEF' (TABTEM, TABKSD, TCP, SDAMD) 

AMD= (F'LAMD*F'LM+SDAMD*SDM) / (F'LM+SDM) 

CALL I NSEF' (TABTEM, TABCF'L, TCP, F'LCF') 

CALL I NSEF' (TABTEM, TABCSD, TCP, SDCP) 

CF'= (PLCF'*PLM+SDCF'*SDM) / (F'LM+SDM) 

C0EFC=2. * (F'LM+SDM) *CF'/3. 14159/AMD 
XPLUS=CQEFC/XMID 
VNUX=0. 2233*XPLUS**0. 7455+3.66 
ALFA X=VNUX # AMD/ RHOLE/2. 

PCONV=ALFAX!KAMID* (HCF'-HW) /CP 
F'OU=F'A-F'CONV+DP-QELOS 
HB=F'OU/ (F'LM+SDM) 

DIFF'=ABS (F'OU - F'B) /F'B 
F'B=POU 

IF (D IFF'. GE. DHMAX) GO TO 200 
END 


C 


10 


SUBROUTINE AVERT (FLOF'L, FLOSD, F'SUM, T ABT , T ABF'L , T ABSD , TVER A ) 
DIMENSION TABT (26) , TABF'L (26) , TABSD (26) 

COMMON MMM 

WRITE (7, 3000) FLOF'L , FLOSD , F'SUM 
NN=0 


TVERA= 10000. 

CONTINUE 
NN=NN+ 1 

IF (NN. GT. 200) THEN 

WRITEOK,*)’ Erro Stop ITER. EXCEED 200 IN AVERT, THEN STOP’ 
END IF 

CALL I NSEF' ( T ABT , T ABF'L , TVERA , EF'L ) 

CALL I NSEF' (TABT, TABSD, TVERA, ESD) 
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PCOMP=FLOPL*EPL+FLOSD*ESD 
DPOWF-1 . -PCOMP/PSUM 
I F ( AB3 f DPOWR ) . GE . 0 . 0 1 0 ) THEN 
TVERA=TVERA* < 1 . +DP0WR/10. ) 

GO TO 10 
END IF 
RETURN 
END 


SUBROUTINE INSEP ( A, B, C, D) 

DIMENSION A (26), B (26) 

COMMON MMM 

IF ( C. LT . A ( 1 ) . OR.C.GT. A (MMM) ) THEN 

WRITE (*,*)’ Erro Stop INDEP. IS BEYOND SCOPE 7 

WRITE (*, *) 7 C= 7 ,C, 7 A(1)= 7 ,A(1), 7 A (MMM) = 7 , A (MMM ) 

STOP 
END IF 

DO 30 K=2 , MMM 

IF ( (C.LE.A(K) . AND.C.GE. A(K-l) ) .OR. 

$ (C.GE.A(K) . AND.C.LE. A(K-l) ) ) THEN 

D=B (K-l) + (B(K> -B CK-1 ) ) / (A (K) -A (K-l ) ) * (C-A(K-1 ) ) 

GO TO 40 
END IF 

50 CONTINUE 
10 CONTINUE 
RETURN 

*010 FORMAT (/5X, 7 INDEP. C= 7 , E12. 4, 5X , 7 DEPEND. D= 7 ,E12.4) 
END 


SUBROUT I NE W I DTH ( OOR I G , QSOUR , HH , NHH , H , ALPHA , VWELD , RMAX , I TW I D ) 
DIMENSION ZZ<34> , TEMP <25, 15, 11) , TEMPO ( 10, 25, 15, 11) 

DIMENSION BEAD (30, 10) , QSOUR ( 10) , QMOLT (10) , HH ( 10) , Q4K ( 10) , QS ( 10) , 
* Z ( 1 1 ) , RMAXO ( 1 i ) , RMAX (11) 

COMMON/DAT 10/NX, NY,NZ, X (25) , Y (20) , TM, TO, THK, EPSL, DENS, HFUSI 
IF ( ITWID. GT. 1 ) GO TO 66 

IF (NX . GT. 1 ) THEN 
DO 1 1=1, NX 

X ( I ) =X ( 1 ) +FLOAT ( I — 1 ) # (X (NX)— X (1) ) /FLOAT (NX-1 ) 

1 CONTINUE 

END IF 

IF (NY. GT. 1 ) THEN 
DO 5 J = 1 , NY 

Y ( J > = Y ( 1 ) +FLOAT ( J-l ) # ( Y (NY) -Y ( 1 ) ) /FLOAT (NY- 1 ) 

5 CONTINUE 
END IF 

NZ=NHH 

IF (NZ. GT. 1 ) THEN 
DO 6 K=1 , NZ 

Z (K) =FLOAT (K-l ) *H/FLOAT (NZ-1 ) 

6 CONTINUE 
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END IF 
NOM= 1 

DO 8 L= 1 , NHH 

HH (L) =FLOAT (L-l > *H/FLOAT (NHH-1) 

8 CONTINUE 



TMOLT=TM-TO 
V=VWELD#0. 0254/60. 

V2A=-V/ (2. *ALPHA) 

N I TER=0 
NZ=NHH+1 
QTOT AL=QQR I G 
DO 10 1 = 1 , NZ 
RNAX ( I ) =0. 0 

Z ( I ) =FLOAT ( 1-1 ) *H/FLOAT (NZ-1 ) 

10 CONTINUE 

DO 40 N= 1 , NHH 
QS < N ) =QSQUR ( N ) 

DO 30 L=1 , NZ 
DO 25 K=1 , NY 

DO 20 J=1 , NX 

TEMPO (N,J,K,L) =0 . 

20 CONTINUE 

25 CONTINUE 

30 CONTINUE 
40 CONTINUE 

DO 65 N= 1 , NHH 
HHH=HH (N) 

ZZ ( 1 ) =-16. *H-HHH 
ZZ (2) =-16. *H+HHH 
DO 50 *>3,34,2 

ZZ (K) =ZZ <K-2) +2. *H 
ZZ (K+l ) =ZZ ( K— 1 ) +2. *H 
50 CONTINUE 

DO 64 K= 1 , NX 
DO 63 L=1 , NY 
DO 62 J = 1 , N Z 
DO 60 M=1 , 34 
ZZZ=Z ( J) -ZZ (M) 

X X X = X <K> 

VYY=Y ( L) 

R=5QRT (XXX**2+YYY>|c*2+ZZZ**2) 

RX=R+XXX 

RX=V2A*RX 

TEMPO < N , K , L , J ) = TEMPO (N, K, L , J ) +EXP (RX ) /R 
60 CONTINUE 

62 CONTINUE 

63 CONTINUE 

64 CONTINUE 

65 CONTINUE 



66 CONTINUE 
NITER-NITER+1 

IF (NITER. GT. 200) THEN 
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Erro Stop 


NITER. GT. 200 IN WIDTH, THEN STOP 


WR I TE ( * , * ) 7 
STOP 
END IF 

DO 63 N= 1 , NHH 
C D4K (N) =QS (N) / (4. *3. 14159 *THK) 

QNOLT <N> = (RMAX (N) +RMAX (N+l) ) *H/FLOAT (NHH) *V*DENS*HFUSI 
Q4K (N) = (QS (N) -QNOLT (N) ) / (4. *3. 14159*THK> 

68 CONTINUE 

DO 75 K=1 , NX 
DO 74 L= 1 , NY 
DO 73 J = 1 , N Z 

TEMP (K,L, J) =0.0 
DO 72 N= 1 , NHH 

TEMP (K,l_, J)=TEMP(K, L, J) +Q4K (N) * TEMPO ( N , K , L , J ) 

72 CONTINUE 

73 CONTINUE 

74 CONTINUE 

75 CONTINUE 

DO 95 K= 1 , NZ 
RMAXO (K) =0. 0 
DO 90 1=1, NX 
DO 80 J=2 , NY 
L=J— 1 

IF ( ( TEMP ( I , L , K ) . LT . TMOLT . AND. TEMP ( I , J, K) .GE.TMOLT) .OR. 

$ (TEMP ( I , L, K) . GT . TMOLT. AND. TEMP ( I , J, K) . LE. TMOLT) ) THEN 

YMOLT=Y(L)+(Y(J)-Y(L) ) / (TEMP C I , J, K) -TEMP ( I , L, K) ) * 

* (TMOLT-TEMP ( I , L, K) ) 

GO TO 85 

END IF 

80 CONTINUE 

85 RMAXO (K) =AMAX 1 (RMAXO (K) , YMOLT ) 

90 CONTINUE 
95 CONTINUE 
DRMAX=0 

DO 100 N=1 , NZ 

DRMAX=AMAX1 ( DRMAX , ABS (RMAX (N) —RMAXO (N) ) ) 

WRITE (*, *> 7 N= 7 , N , 7 2 ♦ RMAX (N)= 7 ,2. *RMAX (N) , 

* 7 2*RMAX0(N)= 7 ,2. *RMAXO(N> 

100 CONTINUE 

WIDTHC=2. JftRMAXO ( 1 ) 

WIDTHR=2. * RMAXO (NZ) 

write'*,*) 7 ROOT WIDTH= 7 , WIDTHR, 7 CROWN WIDTH= 7 , WIDTHC 
wr i te ( 7 , * ) 7 ROOT W I DTH= 7 , W I DTHR , 7 CROWN W I DTH= 7 , W I DTHC 
QRCR0=0. 

QRR00=0. 

DO 120 1=2, NX 
11=1-1 

DO 110 J=2, NY 
J1=J-1 

TRM I D= ( TEMP ( 1 1 , J , 1 ) +TEMP ( 1 1 , J 1 , 1 ) +TEMP ( I , J 1 , 1 ) 

* +TEMP ( I , J , 1 ) ) /4. 

TCMID= (TEMP ( 1 1 , J , NZ) +TEMP ( 1 1 , J 1 , NZ) +TEMF' ( I , J1 , NZ ) 

* +TEMF' ( I , J , NZ ) ) /4 . 

DX = X ( I 1 ) -X ( I ) 
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DY=Y ( J ) — Y ( J 1 > 

DAREA=2 . *ABS <DX*DY> 

QRCRO=QRCRO+ABS (DAREA) *EFSL*0. 00567* ( (TCMID+TO) / 100. ) **4 
QF:ROQ=QRROO+ABS ( DAREA ) *EPSL*0. 00567* ( (TRMID+TO) /100. ) **4 
110 CONTINUE 
120 CONTINUE 

QSQS=QQF: I G-QRCRO-QRROO 
C QSQS=QOR I G 

DQS=1 . -QSQS/QTOTAL 
DB=DRMA X / W I DTHC 

IF (ABS (DOS) . GE. O. 005. OR. DB. GE. 0. 002) THEN 
WRITE (*,*)’ DQS=’ , DOS, ' DB=’,DB 
QT0TAL=0. 0 
DO 130 N= 1 , NHH 

QS <N) =QS (N) * ( 1 . -DQS/QORIG) 

QTOTAL=QTOT AL+QS (N) 

130 CONTINUE 

DO 140 N=1 , NZ 

RMAX (N)=RMAX (N) +. 5* (RMAXO <N) -RMAX (N) ) 

140 CONTINUE 

GO TO 66 
END IF 
RETURN 
END 
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c 

c 

c 

c 

c 

c 

c 
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PROGRAM WELDO. FOR, COMPILED ON MAF:CH 16,1993 

Point and Line Source Models are Combined to Estimate Whether 
a Keyhole Can be Formed. 

Impinging Convective Heat Transfer on the Surface of Workpiece 
Is Taken as Convection before Keyhole Formed. 

D I MENS ION POW 1(5), P0W2 ( 5 ) , P0W3 ( 5 ) , P0W4 (10) 

DIMENSION ENT 1(5) , ENTS (5) , ENT3 (5) , ENT4 ( 10) , TABMF'L (20) 

DIMENSION QL0S1 (10) ,QL0S2(10) ,QL0S3(10) ,QL0S4(10) ,TABMSD(20) 
DIMENSION TEMPI (5) , TEMP2 (5) , TEMP3(5) , TEMP4 ( 10) 

DIMENSION EL0S2 (10) ,EL0S4(10> ,QSEGM(10) ,X SOUR (10) 

D I MENS I ON AREA ( 1 0 ) , TEMF'G (11), DUDL 1(10), DUDL2 (10) , DUDL3 ( 1 A ) 
DIMENSION QSOUR (5) , QNET (5) 

DIMENSION F: 1 ( 1 1 ) , R2 ( 1 1 ) , FLWF’L (11), FLWSD (11), ENERG (11) 
common/dat 10/NX, NY, X (25) , Y(15) , TM, TO, THK, EPSL, DENS, HFUSI 
COMMON MMM ... 

COMMON TABTEIi (20) , TABIF'L (20) , TABCPL (20) , TABKF'L (20) , TAB I . >0) 
COMMON ITMAX, DHMAX, DTMAX 
COMMON TAMBI, DNOZL, FMPL 

COMMON TABCSD (20) , TABKSD (20) , TABNF'L (20) , TABNSD (20) 

COMMON TW4, MREG4, NREG4 

COMMON F'LM.SDM 

OPEN ( 5, F I LE= 7 C : WELDTAB . DAT - ) 

OPEN ( 6 , F I LE= 7 C : WELD I N . USE 7 ) 

OPEN (7, FILE= 7 C: WELDOU. USE 7 ) 

READ (5,*) MMM 

READ (5, #) (TABTEM (I), 1=1, MMM) 

READ (5, *) (TABIF'L (I) , 1 = 1, MMM) 

READ (5,*) (TABCPL ( I ), 1=1 , MMM) 

READ (5,*) (TABKF'L ( I ) , 1 = 1, MMM) 

READ ( 5 , * > ( TAB I SD ( I ) , I = 1 , MMM ) 

READ (5, * ) (TABCSD ( I ), 1=1 , MMM) 

READ (5, #> (TABKSD (I) , 1=1, MMM) 

READ ( 5 , * ) (TABMF'L ( I ), 1 = 1 , MMM) 

READ (5,*) (TABMSD ( I ) , 1=1 , MMM) 

READ (5, *) (TABNF'L ( I ) , 1 = 1 ’ MMM) 

READ (5,*) (TABNSD ( I ), 1=1 , MMM) 

READ (6, * ) NREG1 , NREG2, NREG3, NREG4 

READ (6 , %) VA, VC 

READ ( 6 , * ) FIA1 , FIA2, FIA3 

READ (6,*) TIME, TIMER 

READ (6, * ) DISL1 , D I SL2 , D I SL3, D I SL4 

READ (6,*) WE, WK 

READ (6, * ) AF'IL, AF'ILR, AMAI , AMAIR 

READ (6, *) DNOZL, DSHID, VPL, VSD 

READ (6,*) F'LMOL, SDMOU UR 

READ (6,*) ITMAX, DHMAX, DTMAX 

READ ( 6 , * ) TW 1 , TW2 , TW4 , HSUR2 , HSUR4 
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READ ( 6 , t ) 
READ (6, *) 
READ (6, *> 
READ ( o , t ) 
READ (6, * ) 
READ ( 6 , # ) 
READ (6, * ) 
READ < 6 . * ) 
READ (6, *) 


DXW, TAMBI , DPOLAR, AMDW 
PS I A , F'AMB I , CPSD 
ALFA, VWELD 
NX, NY 
X ( 1 ) , X (NX) 

Y ( 1 ) , Y (NY) 

TM, TO 
THK, EP3L 
DENS, HFUSI 


WRITE (7, 1000) 

WR I TE ( 7 , 1 00 1 ) MMM 
WRITE (7, 1002) 

WRITE (7, 1003) 

WRITE (7, 1004) 

WRITE (7, 1005) 

WRITE (7, 1006) 

WRITE (7, 1007) 

WRITE (7, 1008) 

WRITE (7, 1009) 

WRITE (7, 1010) 

WRITE (7, 1011) 

WRITE (7, 1012) 

WRITE (7, 1013) NREG1 
WRITE (7, 1014) NREG2 
WR I TE ( 7 , 1 0 1 5 ) NREG3 
WRITE (7, 1016) NREG4 
WR I TE ( 7 , 1017) VA 
WRITE (7, 1018) VC 
WRITE (7, 1019) FIA1 
WR I TE ( 7 , 1 020 ) FIA2 
WRITE (7, 1021 ) FI A3 
WRITE (7, 1022) TIME 
WR I TE ( 7 , 1 023) T I MER 
WRITE (7, 1 024 ) D I SL 1 
WRITE (7, 1025) DISL2 
WRITE (7, 1026) DISL3 
WRITE (7, 1 027 ) D I SL4 
WRITE (7, 1028) WE 
WRITE (7, 1029) WK 
WRITE (7, 1030) AF'IL 
WRITE (7, 1031) AF'ILR 
WR I TE ( 7 , 1 032 ) AM A I 
WR I TE ( 7 , 1 033 ) AM A I R 
WRITE (7, 1034) DNQZL 
WR I TE ( 7 , 1 035 ) DSH I D 
WRITE (7, 1036) VF'L 
WRITE (7, 1037) VSD 
WRITE (7, 1038) F'LMOL 
WRITE (7, 1039) SDMGL 
WRITE (7, 1040) UR 
WRITE (7, 1041 ) ITMAX 
WRITE (7, 1042) DHMAX 
WRITE <7, 1043) DTMAX 
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WRITE (7, 1044) TW 1 
WRITE <7, 1045) TW2 
WRITE '.7, 1046) TW4 
WRITE (7, 1 047 ) HSUR2 
WRITE (7, 1048) HSUR4 
WRITE (7, 1049) DXW 
WRITE (7, 1050) TAMBI 
WRITE (7, 1051 ) DF'OLAR 
WRITE (7, 1052) AMDW 
WRITE (7, 1053) PS I A 
WRITE (7, 1054) F'AMBI 
WRITE (7, 1055) CPSD 
WRITE (7, 1056) ALFA 
WRITE (7,1 057 ) VWELD 
WRITE (7, 1058) NX 
WR I TE ( 7 , 1059 ) NY 
WRITE (7, 1060) X ( 1 ) 

WRITE (7, 1061) X (NX) 

WRITE (7, 1062) Y(l) 

WRITE (7, 1063) Y (NY) 

WRITE (7, 1064) TM 
WRITE (7, 1065) TO 
WRITE (7, 1066) THK 
WRITE (7, 1067) EPSL 
WRITE (7, 106B) DENS 
WRITE (7, 1069) HFUSI 



VOLF'L=VF'L/ 1 . 27133 / 1 . E5 
V0LSD=VSD/1 . 27133/ 1 . E5 

FMPL=V0LPL*F'AMBI*PLM0L/8315. /TAMBI 
FMSD=V0LSD*PAMBI*SDM0L/B315. /TAMBI 
UH=V0LSD*4. /3 . 14159/ (DSHID**2-DN0ZL**2> 

FIRST SEGMENT CALCULATION 

AREAW=3. 14159*DP0LAF:**2/4. 

ELOS= 1 . 5IWK/WE/1.E10 

GHW= ( (AMAI+AF'IL) * (VC-FI A1 ) *TIME+ ( AMAIR-AF'ILR) * ( VA+FIA1 ) * 

% TIMER) / (TIME+TIMER) / 1000. 

CQNDW=AMDW*AREAW* (TW1 -TAMBI ) /DXW 
POW 1 ( 1 ) =GHW— CQNDW 
ENT 1(1) =F'0W1 ( 1 ) /FMF’L 

CALL INSEP (TAB I PL, TABTEM, ENT1 (1) , TEMPI (1) ) 

EL0S1=EL0S#TEMF'1 ( 1 ) * (AMAIR-AF'ILR) *TIMER/ (TIME+TIMER) 

AMF'CP 1 = ( (AMAI+AF'IL) *TIME+ (AMAIR-AF'ILR) *TIMER) / (TIME+TIMER) 
DX 1 =DI SL1 /FLOAT (NREG1) 

GE0L1=3. 14159* ( ( DF'OLAR+DNOZL ) /2. ) **2/4. /DX1 
V0LUM1=3. 14159* ( ( DNOZL+DF'QLAR ) /2. ) **2/4. *DX1 
MREG 1 =NREG 1+1 
SUMQ 1 =0 . 

DUBC=0. 

DO 40 M=2, MREG1 
N=M— 1 

DUDL1 (N) =2000. 
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QLQS 1 ( N ) =0. 

N N 1=0 

10 CONTINUE 
NN1 =NN 1+1 

IFlNNl.GT.200) THEN 

WR I TE <*,*) ’ NN 1 . GT . 200 . Erro stop ' 

STOP 
END IF 

DF'1=DUDL1 (N) *DX1 *AMF'CP1 / 1000. -QL0S1 (N) 

P0W1 (M) =P0W1 (N) +DP 1 

ENT 1 (M)=P0W1 (M) /FMPL 

ENTCP 1 = ( ENT 1 (N) +ENT1 CM) ) /2. 

CALL I NSEP ( TAB I PL , TABTEM , ENTCP 1 , TEMCP 1 ) 

DL0S1 (N) =V0LUM1*4. * 1 . E7* ( TEMCP 1 / 15000. ) **16/ 
$ ( 1 . -M TEMCP.l / 1 5000 . ) **16) 

CALL I NSEP ( TABTEM , T ABMPL , TEMCP 1 , PLMHO ) 

OHM 1=1 . /PLMHQ/GEOL 1 
DUDL=AMPCP 1 *OHM 1 /DX 1 
DIFU= 1 . -DUDL/DUDL1 (N) 

IF CABS (DIFU) . GE. 0.005) THEN 
DUDL1 ( N ) = ( DUDL+DUDL 1 ( N ) ) / 2 . 

GO TO 10 
ENDIF 

DUBC=DUBC+DUDL 1 ( N ) * D X 1 
SUMQ 1 =SUMQ1 +QL0S 1 (N) 

CALL I NSEP ( TAB I PL , TABTEM , ENT 1 CM) , TEMP 1 CM) ) 

40 CONTINUE 


C WRITE (7, 1 170) ( ROW 1 ( I ) , I = 1 , MREG 1 ) 

C WRITE (7, 1 180) ( ENT 1 ( I ) , I = 1 , MREG 1 ) 

C WRITE <7, 1 190) (QL0S1 (I) , 1=1 ,NREG1) 

C WRITE (7, 1200) (TEMPI C I ) , 1=1 , MREG1 ) 

C WRITE (7, 1220) ( DUDL 1 C I ) , I = 1 , NREG 1 ) 

C WRITE (7, 1210) SUMQl.DUBC 

C 

C SECOND SEGMENT CLCULATION - 

C 


PQW2 C 1 ) =PQW 1 C MREG 1 ) 

ENT2 ( 1 ) =ENT 1 ( MREG 1 ) 

TEMF'2 ( 1 ) =TEMP 1 (MREG1 ) 

AMPCP2= (AM A I *TIME+ ( AMAIR-APILR) *TIMER) / (TIME+TIMER) 

DX2=D ISL2/FL0AT (NREG2) 

AF:EA2=3. 14159*DN0ZL*DX2 
GE0L2=3 . 1 4 1 59*DN0ZL* *2/ 4. / D X 2 
MREG2=NREG2+1 
C WRITE (7, 1230) 

GHN1=(APIL*FIA2*TIME+APILR*FIA2*TIMER) / (TIME+TIMER) / 1000. 
GHN2=(APIL*VA*TIME+APILR*VA*TIMER) / (TIME+TIMER) / 1000. 
SUMQ2=GHN 1 +PS I A*GHN2 
DUCG=0 . 

DO 80 M=2 , MREG2 
N=M- 1 

XM=FLOAT (N) *DX2-DX2/2. 

DUDL2 ( N ) =2000 . 

QL0S2 (N) =0. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 


NN2=0 


50 CONTINUE 
NN2=NN2+ 1 

I F ( NN2 . GT . 200 ) THEN 

WRITE (*, *) ’ Erro Stop ENT3=-’ , ENTT 

STOP 
END IF 

IF (NREG2. EQ. 1 ) then 

EL0S2 (N> =EL0S*TEMP2 (N) * ( AP I L * T I ME+AP I LR * T I MER ) / (TIME+TIMEP' 
else IF (NREG2. NE. 1 . AND. N. ED. 1 ) then 

EL0S2 (N> =ELQS*TEMF'2 (N) KAPIL* TIME/ (TIME+TIMER) 
else IF ( NREG2 . NE. 1 . AND. N. EQ. NREG2) then 

EL0S2 (N) =EL0S*TEMP2 (N) *AF'ILR*TIMER/ (TIME+TIMER) 
endi-f 

DP2=DUDL2 (N) *DX2* AMPCP2/ 1 000 . + ( 1 . -PSIA) *GHN2/FL0AT (NREG^) 

* -ELQS2 ( N ) -QL0S2 ( N ) 

P0W2 ( M ) =P0W2 ( N ) +DP2 
ENT2 ( M ) =P0W2 (M) /FMF'L 
ENTCF'= (ENT2 (N) +ENT2 (M) ) /2. 

CALL INSEP (TABIF'L, TABTEM, ENTCF', TEMCP) 

CALL INSEP (TABTEM, TABKPL, TEMCP, AMD) 

CALL INSEP (TABTEM, TABCPL, TEMCP, CP) 

CALL INSEP (TABTEM, TABMF'L, TEMCP, F'LMHO) 

C0EFC=2. *FMPL*CP/3. 14159/AMD 
XPLUS=C0EFC/XM 

VNUX=0. 2233*XPLUS**0. 7455+3. 66 
ALFAX=VNUX # AMD/DNOZL 

QL0S2 (N) =ALFAX*AREA2* (ENTCP-HSUR2) /CP 
0HM2= 1 . / F'LMHO /GE0L2 
DUDL=AMF'CF'2*0HM2/DX2 
DIFU=1 . -DUDL/DUDL2 (N) 

IF (ABS (DIFU) . GE. 0. 005) THEN 
DUDL2 (N) = (DUDL+DUDL2 (N) ) /2. 

GO TO 50 
END IF 


DUCG=DUCG+DUDL2 (N) *DX2 
SUMQ2=SUMQ2+QL0S2 (N) +EL.0S2 (N) 
CALL INSEP (TABIF'L, TABTEM, ENT2 (M) 
SO CONTINUE 

(P0W2 ( I ) , 

( ENT2 ( I ) , 

( QLQS2 ( I ) 

(TEMP2 ( I ) 

( EL0S2 ( I ) 

( DUDL2 ( I ) 


TEMF'2 (M> ) 


WRITE (7, 1250) 
WRITE (7, 1260) 
WRITE (7, 1270) 
WRITE (7, 1280) 
WRITE (7, 1290) 
WRITE (7, 1300) 
WRITE (7, 1310) 
WRITE (7, 1320) 


SUMQ2, DUCG 


1=1 , MREG2) 

1 = 1, MREG2) 

, 1 = 1 , NREG2) 
, 1 = 1 , MREG2) 
, 1 = 1, NREG2) 
1 = 1, NREG2) 


THIRD SEGMENT CLCULATION 


R0=DN0ZL/2. 

P0W3 ( 1 ) =F'0W2 ( MREG2 ) 
MITEF:=0 
90 CONTINUE 
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M I TER=M I TER+ 1 
IF (MITER. GT. 200) THEN 

WR I TE ( * , * ) 7 Erro Stop MITER. GT . 200, DPP 35 ' , DPP 

STOP 
END IF 

ENT 3 ( 1 ) =P0W3 ( 1 ) /FMF'L 

CALL I NSEP ( TAB I PL , TABTEM , ENT3 ( 1 > , TEMP3 ( 1 ) ) 

DENSO=F'AMB I #PLMQL/S3 1 5 . /TEMP3 ( 1 ) 

U0=FMF'L/ DENSO /3. 14159/R0**2 
ENERG ( 1 )=FMPL*U0#*2/2. /1000. 

PPP-P0W2 (MREG2) -ENERG ( 1 ) 

DF'P=ABS ( 1 . -PPP/P0W3 ( 1 ) ) 

P0W3 ( 1 ) =PPP 

I F ( DPP . GE . 0 . 005 ) GO TO 90 
COEFM=UH/UO 

VMID1= ( 1 . — COEFM) / (l.+COEFM) 

VMID3=1 . +0. 8*C0EFM-0. 45*C0EFM**2 
DX3=D I SL3/ FLOAT (NREG3) 

MREG3=NREG3+1 

SUMQ3=0. 

DUGE=0. 

R1 ( 1 ) =R0 
R2 ( 1 ) -R0 
FLWPL ( 1 ) =FMPL 
FLWSD ( 1 ) =0. 

AMPCP3=(AMAI*TIME+AMAIR*TIMER> / (TIME+TIMER) 

DO 190 M=2 , MREG3 
N=M-1 

X3=FL0AT (N) *DX3 
VMID2=0.27*X3 

VM I D4=VM I D 1 * (0.41 6+0 . 1 34*C0EFM+0 . 02 1 # VM I D2* VM I D 1 * VM I D3 /RO ) 

R1 (M)=R0-VMID2*VMID4 

VM I DA=0 . 02 1 * VM I D3/R0 

VMIDB=0. 416+0. 1 34#C0EFM 

VMIDC=R1 (M) -RO 

VM I DD= ( - VM I DB+SQRT (VM 100**2-4. *VMIDA*VMIDC> ) /2. /VMIDA 

F:2 ( M ) =F: 1 (M) +VMIDD 

N0M=0 

TBEG I =TEMP3 ( 1 ) 

DENBEG=DENSO 

DR= ( R2 ( M ) -R 1 (M) ) /10. 

DUDL2 (N) =2000. 

GE0L3=3. 141 59KR2 ( M ) **2/DX3 
120 CONTINUE 
NOM=NOM+ 1 

IF (NOM. GT. 200) THEN 

WRITE (*,*>’ Erro Stop ITER. IN STANDOFF EXCEED 200’ 
STOP 
END IF 

FLWPL (M) =3. 14159*R1 (M) * *2*U0*DENBEG 
FLWSD (M) =0. 

ENERG (M) =FLWFL (M) *U0**2/2. /1000. 

CALL INSEP (TABTEM, TAB I PL , TBEG I , EBEGI ) 

POWER=FLWPL (M) *EBEGI 
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DO 150 K= 1 , 10 

F:M ID-FLOAT (K) * DR-DR/2. +R1 (M) 

AT A = ( R2 < M ) -ft II ID) / (R2(M) -R 1 (M) ) 

UMID=UO- (UO-UH) * ( 1 . -ATA** 1 . 5) **2 
TMID=TBEGI- (TBEGI-TAMBI ) * ( 1 . -ATA) 

F'GPL=PAMB I / ( 1. + (1. / ATA— 1 . ) *F'LMOL/SDMOL ) 

F'GSD=PAMB I -F'GF'L 
DENPL=PGF'L*PLM0L/8315. /TMID 
DENSD=PGSD*SDM0L/8315. /TMID 
V0LUM=2. *3. 14159*RMID*DR*UMID 
FLWPL (M) =FLWPL ( M ) +VOLUM*DENPL 
FLWSD < M ) =FLWSD ( M ) +VOLUM*DENSD 

ENERG ( M ) =ENERG (M) + (DENPL+DENSD) *V0LUM*UMID**2/2. /1000. 

IF ( TM I D . GE . 3000 . ) THEN 

CALL I NSEP < T ABTEM , TAB I PL , TM I D , ENTF'L ) 

CALL I NSEP ( T ABTEM , TABI SD , TMID, ENTSD) 

ELSE 

ENTF'L = T ABCF'L ( 1 ) *TMID 
ENTSD=T ABCSD < 1 ) *TMID 
ENDIF 

POWER=POWER+VOLUM* (DENPL*ENTPL+DENSD*ENTSD) 

150 CONTINUE 

F'OWRCP= (P0WER+R0W3 (N) ) /2. 

PLCF'= (FLWPL (M> +FLWPL ( N ) ) / 2 . 

SDCP= < FLWSD ( M ) +FLWSD ( N ) ) /2. 

CALL AVERT ( F'LCP , SDCP , PQWRCP , T ABTEM , TAB I F'L , TAB I SD , TCP ) 

VRAD I =3. 14159* (R2(M)+R2«N) ) **2/4. *DX3 

QLDS3 (N) =VF;ADI *4. * 1 . E7* (TCP/ 15000. ) **16/ ( 1 . + (TCP/ 15000. ) **16) 
CALL INSEP (TABTEM, TABMF'L, TCP, PLMHO) 

CALL INSEP (TABTEM, TABMSD, TCP, SDMHO) 

CPMHO= ( (FLWPL (M) +FLWPL (N) ) *PLMHO+ (FLWSD (M) +FLWSD (N) ) *SDMHQ) / 
$ ( FLWPL ( M ) +FLWPL ( N ) +FLWSD ( M ) +FLWSD ( N ) ) 

0HM3= 1 . /CPMHQ/GEQL3 
DUDL3 (N) =0HM3* AMPCP3/DX3 
DP3=0HM3 * AMPCP3 * *2/ 1 000 . 

P0W3 (M) =F'0W3 (N) +DP3-QL0S3 (N) + (FLWSD (M) -FLWSD (N) ) *CPSD*TAMBI- 
% ( ENERG ( M ) -ENERG (N) ) 

F’MF'= 1 . — P0WER/P0W3 (M) 

IF ( ABS (F'MP) .GE. 0.010) THEN 
TBEG I =TBEG I * ( 1 . +PMP/10. ) 

DENBEG=PAMBI*PLM0L/B315. /TBEG I 
U0=U0*FMPL/ FLWPL (M) 

GO TO 120 
END I F 

WRITE (*, 1390) N , NOM , F’MP 
SUMQ3=SUMQ3+QL0S3 (N) 

CALL AVERT (FLWPL (M) , FLWSD (M) , POWER, TABTEM, TABI F'L, TABI SD, 

$ TEMF'3 (M) ) 

. ENT3 ( M > =F'0W3 ( M ) / (FLWPL (M) +FLWSD (M) ) 

DUGE=DUGE+DUDL3 (N) *DX3 
190 CONTINUE 

PAID=FLWSD (MREG3) *CF’SD*TAMBI 
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GHK1=(AMAI*FIA3*TIME+AMAIR* (UR-FIA3) *TIMER> / (TIME+TIMER) /1000. 
GHK2= (AMAI *VA*TIME+AMAIR*VC*TIMER) / (TIME+TIMER) / 1000. 

GHK 1 : = ( AM A I tF I A3* TIME- AM A I R#FI A3*TIMER ) / (TIME+TIMER) /1000. 

GHK2= (AMAI *VA*TIME+AMAIR* (UR+VC) *TIMER> / (TIME+TIMER) /1000. 


P0WER=PDW3 (MREG3) 

PLM=FMPL 
SDM=FLWSD (MREG3) 

PQW4 ( 1 > =POWER+ ( 1 . -PS I A ) *GHK2 
ENT4 ( 1 ) =P0W4 ( 1 ) / ( PLM+SDM ) 

CALL AVERT (PLM, SDM, PQW4 < 1 ) , TABTEM, TABIF'L, TABISD, TEMP4 ( 1 ) ) 
EL0SS=EL0S*TEMP4 ( 1 ) * AMAI* TIME/ (TIME+TIMER) 

GHK 1 2=GHK 1 +PS I A*GHK2 

CALL CON VO (R2 (MREG3) , ENT 4(1) , TEMP4 ( 1 ) , PLM, SDM, HSUR4, QCONV ) 
SUMQ4=C!C0NV+EL0SS+GHK12 
WRITE (*,*> ’ SUMQ4=’ , SUMQ4, ’ VWELD=’ , VWELD 
CALL DEPTH CSUMQ4, DISL4, NREG4, ALFA, VWELD) 

STOP 


-++++++++ 


1000 FORMAT (/3X, ’ SYMBOL UNIT 

$’ DATA’/) 

1001 FORMAT ( /3X , ’ MMM 
$15) 

1002 FORMAT (/3X, ’ TABTEM K 
$’ Table attached’) 

1003 FORMAT </3X, ’TABIF'L kJ/kg 
$’ Table attached’) 

1004 FORMAT ( /3X TABCF'L kJ/kg. K Specific heat of plasma gas 
$’ Table attached’) 

1005 FORMAT ( /3X , ’ TABKF'L kW/M. K 
$’ Table attached’) 

1006 FORMAT (/3X, ’ TABISD kJ/kg 
$’ Table attached’) 

1007 FOF:MAT ( /3X , ’ TABCSD kJ/kg. K Spec i -fie heat o-f shield gas 
$’ Table attached’) 

1008 FORMAT ( /3X ’ TABKSD kW/M. K 

$ ’ Table attached’) 

1009 FORMAT (/3X, ’ TABMPL MHO/M 
$’ Table attached’) 

1010 FORMAT (/3X, ’ TABMSD MHO/M 
$’ Table attached’) 

101 1 FORMAT ( / 3X , ’ TABNF'L kg/M. S 
$’ Table attached’) 

1012 FORMAT (/3X, ’ TABNSD kg/M.S 
$’ Table attached’) 


DESCRIPTION 

Number o-f point in table 
Temperaeture table 
Enthalpy o-f plasma gas 


Thermal conductivity o-f plasma gas ’ 
Enthalpy o-f shield gas ’ 


Thermal conductivity o-f shield gas ’ 
Elictric conductivity o-f plasma gas’ 
Electric conductivity o-f shield gas’ 
Viscosity o-f plasma gas ’ 
Viscosity o-f shield gas ’ 


1013 FORMAT (/3X, ’NREG1 

Segment 

number 

between electrode 

$ /3X , ’ 

and ori-fice 



? , 15) 

1014 FORMAT ( /3X , ’ NREG2 
$15) 

Segment 

number 

o-f 

orifice 

1015 FORMAT </3X, ’ NREG3 
$15) 

Segment 

number 

o-f 

standof -f 

1016 FORMAT ( /3 X „ ’ NREG4 
$15) 

Segment 

number 

o-f 

wor kpi ece 
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1017 FORMAT ( / 3X , 7 Vfi 
Sel2. 5) 

1 0 1 3 FORMAT C / 3 X , VC 
*e 12.5) 

1019 FORMAT ( / 3 X , 7 F I A 1 
*el2. 5) 

1020 FORMAT ( /3X , 7 FIA2 
$el2.5> 

1021 FORMAT ( /3X , 'FI A3 
*el2.5> 

1022 FORMAT (/3X, 7 TIME 
*el2. 5) 

1023 FORMAT (/3X, 7 TIMER 
*e 12.5) 

1024 FORMAT ( /3X, 7 DISL1 
*/3X , 7 

1025 FORMAT (/3X, 7 DISL2 
*el2.5> 

1026 FORMAT ( /3X , 7 DISL3 
Sel2.5> 

1027 FORMAT ( /3X , 'DISL4 
*el2. 5) 

1028 FORMAT (/3X, 7 WE 
*el2.5> 

1029 FORMAT ( /3X , 7 VVK 
*el2.5> 

1030 FORMAT (/3X, ’ARIL 
*/3x , 7 

1031 FORMAT ( /3X, 7 API LR 
*/3x, 7 

1032 FORMAT </3X, 7 AMAI 
*/3x, 7 

1033 FORMAT </3X, 7 AMAIR 
*/3x, 7 

1034 FORMAT ( /3X , 7 DNOZL 
*el2.5> 

1035 FORMAT (/3X, 7 DSHID 
$el2. 5) 

1036 FORMAT ( /3X , 7 VPL 
*e 12.5) 

1037 FORMAT </3X, 7 VSD 
*el2. 5) 

1038 FORMAT </3X, 7 PLMOL 
$e 12.5) 

1039 FORMAT ( /3X , 7 SDMOL 
$el2. 5) 

1040 FORMAT ( /3X , 7 UR 
$e 12.5) 

1041 FORMAT </3X, 7 ITMAX 
*15) 

1042 FORMAT (/3X, 7 DHMAX 
Se 12. 5 ) 

1043 FORMAT (/3X, 7 DTMAX 
*el2. 5) 


Anode potential drop 
Cathode potential drop 


Electrode work -function 


Ori-fice work function 


Workpiece work function 


mS Time duration for reverse polarity 7 , 

M Distance betwteen electrode and 

orifice 7 , E12. 5) 

M Axial length of orifice 

M Distance of standoff 7 , 

M Thickness of workpiece 7 , 

Coulomb Electron charge 7 , 

erg/K Boltzman constant 7 , 

A Pilot electric current in straight 7 , 

polarity 7 ,el2.5) 

A Pilot electric current in reverse 7 , 

polarity 7 ,,el2.5) 

A Main electric current in straight 7 , 

p o 1 ari ty 7 , el2. 5) 

A Main electric current in reverse 7 , 

pol ar i ty 7 , el2. 5) 

M Orifice diameter 

M External diameter of shield 7 , 


scf h 


Plasma gas flow rate 


scf h 


Shield gas flow rate 


kg/kmol Molecular weight of plasma gas 7 , 

kg/kmol Molecular weight of shield gas 7 , 

V Reverse polarity voltage rise 7 , 

Maximum number of iteration 7 , 

Maximum relative error of enthalpy 7 , 


Maximum error of temperature 
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1 044 FORMAT ( / 3 X , 7 T W 1 
*el2. 5) 

K 

1045 FORMAT < /3X , 7 TW2 
$e!2. 5) 

K 

1046 FORMAT ( /3X , 7 TW4 
*el2.5) 

K 

1047 FORMAT </3X, 7 HSUR2 
$el2. 5) 

kJ/kg 

1048 FORMAT ( /3X , 7 HSUR4 

kJ/kg 


104? FORMAT C/3X, 7 DXW M 

*el2.5) 

1050 FORMAT ( /3X , 7 T AMD I K 

$el2. 5) 

1051 FORMAT (/3X, 7 DPOLAR M 
$e 12. 5 ) 

1052 FORMAT (/3X, 7 AMDW kW/M.K 

$el2. 5) 

1053 FORMAT </3X, 'PSIA 
$/3x, 7 

1054 FORMAT (/3X, ’FAME I 
*el2.5> 

1055 FORMAT ( /3X , 7 CPSD 
*/3>:, 7 

1056 FORMAT C/3X, 7 ALFA 
$e 12.5) 

1057 FORMAT C/3X, 7 VWELD ipm 
*el2.5) 

1058 FORMAT </3X, ’NX 
$/3 x , 7 

1059 FORMAT ( /3X , 7 NY 

1060 FORMAT ( /3X , 7 X ( 1 > M 
$/3:-: , 

1061 FORMAT (/3X, 7 X (NX) M 
$/3x, 7 

1062 FORMAT ( /3X , 7 Y ( 1 ) M 
$/3;-c , 7 

1063 FORMAT ( /3X , 7 Y (NY) M 

$/3x , 7 

1064 FORMAT ( /3X , 7 TM K 

*el2.5) 

1 065 FORMAT ( / 3 X , 7 TO K 
$el 2. 5) 

1 066 F ORM AT ( / 3 X , 7 THK k W / M . K 

*el2.5> 

1067 FORMAT </3X, 7 EPSL 
$el2. 5) 

1068 FORMAT (/3X, 7 DENS kg/M3 

$e 12.5) 

1069 FORMAT </3X, 7 HFUSI kJ/kg 

$el2. 5) 


El ec trodt? wall temperature 7 , 

Orifice wall temperature 7 ? 

Workpiece wall temperature 7 

Plasma enthalpy correponding to TW2 7 , 
Jet enthalpy correponding to TW4 7 , 

Electrode length 

Enviromental temperature 7 , 

Electrode diameter 7 

Thermal conductivity of electrode 

Fraction power loss due to anode 7 , 

and cathode drop 7 ,el2.5) 

N/M2 Envi ronmental pressure 7 

kJ/kg. K Specific heat of shield gas with 
environmental temperature 9 ,612-5) 

kW/M.K Thermal dittusivity ot workpiece % 


Welding speed 


Number at computing nodal in X 
coordi nate * 15 ) 

Number o-f computing nodal in V 
coordi nate ", 15) 

Beginning position o-f computed 
nodal in X coordinate 


Ending position of computed 


nodal in X coordinate 


Beginning position of computed 


7 ,012.5) 
7 , el2. 5) 


nodal in Y coordinate 


Ending position of computed 


nodal in Y coordinate 


7 , e!2. 5) 
7 , 012 . 5 } 


Melting temperature at workpiece 
Initial temperature o-f workpiece 
Thermal conductivity at workpiece 
Emissivity ot workpiece 
Density ot workpiece 
Melting latent heat ot workpiece 


1170 FORMAT (/2X, ’F0W1*’ , 4E12. 5, /7X , 4E12. 5) 
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11 30 FORMAT ( / 2 X , ENT 1 = ’ , 4E 12 . 5 , /7X , 4E 1 2 . 5) 

1190 FORMAT </2X, ’QL0S1=’ ,4E12.5, /7X,4E12.5) 

1200 FORMAT </2X, ’TEMPI-’ , 4E1 2. 5, /7X , 4E12. 5) 

1210 FURMAT (/2X, ’ TOTAL LOSS IN GAP I , E 1 2 . 5 , 2 X , ’ DUBC ( v ) ’ , F6 . 
1220 FORMAT ( /2X , ’ DUDL 1 ’ , 4E 12. 5, /7X , 4E1 2 . 5 ) 

1300 FORMAT < /2X , DUDL2 , 4E12. 5, /7X , 4E12. 5) 

1310 FORMAT ( /5X , ’ TOTAL LOSS IN NOZZLE = ’ , E12. 5, 2X , ’ DUCG ( v) ’ , 
1 250 FORMAT ( /5X , ’ F'0W2= ’ , 4E 12.5/1 OX , 4E 12. 5 ) 

1260 F0RMATC/5X, ’ENT2=’ , 4E 12. 5/ 10X , 4E1 2. 5) 

1270 FORMAT ( / 4 X , ’QL0S2-’ , 4E12. 5/ 10X , 4E12. 5) 

1 280 FORMAT ( / 2X , ’ TEMF'2 ’ , 4E1 2 . 5 , / 7X , 4E 1 2 . 5 ) 

1290 FORMAT ( / 2X , ’EL0S2’ , 4E12. 5, /7X, 4E12.5) 

1320 FORMAT <//5X, ’ STANDOFF COMPUTATION’ > 

1230 FORMAT (//5X, ’ NOZZLE COMPUTATION’) 


1240 FORMAT 
1340 FORMAT 
1350 FORMAT 
1370 FORMAT 
1360 FORMAT 
1380 FORMAT 
1390 FORMAT 
$ 

1510 FORMAT 
1450 FORMAT 
1410 FORMAT 
1430 FORMAT 
1440 FORMAT 
1460 FORMAT 
1470 FORMAT 
1480 FORMAT 
1490 FORMAT 
1500 FORMAT 
1520 FORMAT 
1530 FORMAT 
1540 FORMAT 
1550 FORMAT 


/5X 
/5X 
/2X 
/IX 
IX, 
/2X 
/2X 
2X 
/2X 
/2X 
/5X 
/4X 
/ 4X 
/5X 
/5X 
/5X 
/5X 
/5X 
/5X 
/2X, 
/5X , 
/5X 
//I 


, ’NN2.GT. 100 ’ ) 

, ’ UH, UO, COEFM ’ , 3E12. 5) 

, ’LINE 285 TBEGI ’ , E12. 4) 

, ’ ATA, TBEGI , TAMBI , TMID’ , 4E11.4) 

’ M , R 1 , R2 , DR , RM I D , AT A ’ , 12, 2E12.5, /5X,3E12.5> 

, ’ IN MAIN, NOM, MAR, MHE, POWER’ , 14, 2X, 3E10. 3) 

, ’ ITER. IN STANDOFF, SEG.=’ , 12, 2X, ’NOM=’ , 14, 

, ’PMF'=’ ,F10.6> 

, ’ DUDL3’ , 4E12. 5, /7X, 4E12. 5) 

, ’TOTAL LOSS IN STANDOFF , E12. 5, 2X, ’ DUGE (v) 
, ’P0W3=’ , 4E12. 5, /10X, 4E12. 5) 


’ QL0S3= 
’ TEMP 3= 
’ R 1 ( m ) = 
’ R2 (m) = 
’ FLWF'L = 
’ FLWSD= 
’ ENERG= 
’POWER 


, 4E 12.5, /10X, 4E12.5) 
, 4E12. 

, 4F 1 0 . 

, 4F10. 

4E12. 5. 

4E12. 5. 


,5, /10X 
8 , / 1 1 X 
8 , / 1 1 X 
5, /I IX 
/I IX 


VELOC. , TEMP. IN 
T OF BOUNDARY, 

, ’CONVECTIVE LOS 
OX, ’ WORKPIECE 


4E12.5) 

4F10. 8) 

4F10.8) 

4E12.5) 

4E12.5) 

, 4E12.5, /I IX, 4E12. 5) 

ID FROM SHIELD, PAID (kW> =’ , E12. 5) 


CORE AT LAST’ , 2E1: 
TBOUND=’ , El 2.5) 
OSS OF JET ( k W ) ’ , E 1 2 
COMPUTATION’) 


1560 FORMAT (/2X, ’ XSOUR’ , 2X, 5F3. 6) 

1570 FORMAT </5X, ’RADIUS OF KEYHOLE (m)=’,E12.5) 

1580 FORMAT ( /2X , ’ AREA (m**2) ’ , 4E12.5, /12X,4E12.5) 

1590 FORMAT ( / 4X , ’ PLM, SDM , TEMF'4 ( 1 ) ’ ,3E14.5) 

1600 FORMAT C/5X, ’ IMPINGE LOSS (kw)’,F8.5) 

1660 F0RMAT(/2X, ’ QSEGM= ’ , 4E12. 5, /8X , 4E12. 5) 

1680 FORMAT (/5X, ’TOTAL LOSS IN WORKPIECE =’,E12.5) 

1610 FORMAT ( / 2X , ’ F'0W4= ’ , 4E 12.5, /7X , 4E 12 . 5) 

1620 FORMAT </2X, ’ ENT4=’ ,4E12.5, /7x’ 4E12.5) 

1630 FORMAT (/2X, ’ TEMF'4= ’ ,4E12.5, /8X,4E12.5) 

1640 FORMAT </2X, ’QL0S4=’ , 4E 1 2 . 5, /8X , 4E 1 2. 5) 

1650 FORMAT </2X, ’EL0S4=’ ,4E12.5, /8X,4E12.5) 

1670 FORMAT (/2X, ’ QS0UR=’ , 4E12. 5, /8X, 4E12. 5) 

1690 FORMAT (5 X, ’ TOTAL POWER INCLUDING HELIUM, PT0TIN=’ , E12. 5 
* 5X,’ TOTAL POWER FROM LOSS, F'TOTOU=’ , E12. 5) 

1700 FORMAT (/2X, ’POLARITY CONDUC. LOSS AND PERCENT. ’ , 2E 1 2 . 5 ) 
1710 FORMAT </2X, ’GAP LOSS AND PERCENTAGE’ , 2E 12. 5) 


1 1 


2 ) 

F6. 2) 


’ , F6.2) 
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1 720 
1730 
1 740 
1750 
1 760 
1 770 
1 780 
1 790 
1 800 


FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
END 


/ 2 X 

/ •— i 


/2X 

/5X 


/5X, 
/ 4 X , 


NOZZLE 
RAD I AT. IN 
CONVEC. IN 
WORKPIECE 
DISCHARGE 


SS AND PERCENTAGE’ ,2E 12. 5) 


STANDOFF 
STANDOFF 
LOSS AND 
LOSS AND 
SUM. PERCENTAGE’ , F8. 5) 

AMF'CF' 1,2,3’ , 3E 12.5) 

DUBE, DUAF, DUAFR’ , 3E12. 5) 
DUF'F= ’ , E12. 4, 4X , ’ DUF'R=’ .E12 


AND PERCENTAGE 
AND PERCENTAGE 
PERCENTAGE’ , 2E 
PERCENTAGE’ , 2E 


4) 


’ , 2Ei: 
’ , 2E 1 : 
12 . 5 ) 
12.5) 


SUBROUTINE AVERT (FLOF'L, FLOSD, F'SUM, TABT, TABPL, TABSD, TVER A) 
DIMENSION TABT (20) , TABPL (20) , TABSD (20) 

COMMON MMM 

WR I TE ( 7 , 3000 ) FLOPL , FLOSD , F'SUM 
NN=0 

TVERA= 10000. 

10 CONTINUE 
NN=NN+ 1 

IF(NN. GT. 200) THEN 

WRITE (*,*) ’ Erro Stop ITER. EXCEED 200 IN AVERT, THEN STOP’ 
STOP 

END IF 

CALL I NSEF' ( TABT , TABPL , TVER A , EF'L ) 

CALL I NSEF' (TABT, TABSD, TVER A, ESD) 

PCOMP=FLOPL*EPL+FLOSD*ESD 
DP0WF:=1 . -F'COMF'/ F'SUM 
I F ( ABS ( DF'OWR ) . GE. 0 . 0 1 0 ) THEN 
TVERA=TVERA* ( 1 . +DP0WR/10. ) 

WR ITE (* , 3020 ) NN , DF'OWR 
GO TO 10 
END IF 

WRITE (7, 3030) NN, DF'OWR, T VERA 
RETURN 

3000 FORMAT </2X, ’ IN AVERT, MF'L, MSD, F'SUM’ , 3E1 1 . 4) 

3020 FORMAT (2X, ’ NN=’ , 13, 2X, ’ DF'OWR=’ ,E12.5) 

3030 FORMAT (/2X, ’SUB. AVERT, NN, DF'OWR, TVERA’ , I4,2E14.5> 

END 


SUBROUTINE I NSEF' (A, B, C, D) 

DIMENSION A (26) , B (26) 

COMMON MMM 

IF(C. LT. A ( 1 ) . OR. C. GT. A (MMM) ) THEN 

WRITE (*,*) ’ Erro Stop INDEF' . IS BEYOND SCOPE’ 

WRITE OK, *) ’ C=’ , C, ’ A ( 1 ) =’ , A ( 1 ) , ’ A (MMM) =’ , A (MMM ) 

STOP 
END IF 

DO 30 K=2, MMM 

IF ( (C.LE. A(K) . AND.C.GE. A(K— 1) ) .OR. 

$ (C.GE.A(K) . AND.C.LE.A(K-l) ) ) THEN 


F'age 


12 



n n n 


c 


A ( K— 1 ) ) * (C-A (K-l ) ) 


D=B ( K— 1 > + ( B ( K > — B ( K— 1 ) ) / ( A ( K ) - 
WRITE (*, 4010) C,D 
GO TO 40 
END IF 
CONTINUE 

C WRITE (*,4020) 

40 CONTINUE 
RETURN 

4020 FORMAT ( / 5X , ’ WRITE FROM FORMAT 31 IN SUB INSERT 
4010 FORMAT (/5X, ’ INDEP. C=’ , E12. 4, 5X, ’ DEPEND. D=\,E12.4) 

END 


1 


SUBROUTINE DEPTH (QORIG, H, NHH, ALPHA, VWELD) 

DIMENSION Z (20) , ZZ (34) , HH(15) , TEMP (25, 10, 10) 

DIMENSION XMOLT (50) , YMOLT (50) , ZMOLT (50) , ZMY (50) D7(15) 
* DHZ (15) , XMOLTO (50) , ZM0LT0(50) , TEMPO (10, 25, 10, 10) *" 
common/dat 10/NX, NY, X (25) , Y ( 15) , TM, TO, THK, EPSL, DENS, HFUSI 
OPEN (8, FILE= : ’ C: WEL. 0U1 ' ) ' 

OPEN ( 9 , F I LE= •’ C : WEL . 0U2 •’ ) 

OPEN ( 1 0 , F I LE= ’ C : WEL . 0U3 ’ ) 

IF (NX. GT. 1 ) THEN 
DO 1 1=1, NX 

X ( I ) = X ( 1 ) +FLOAT ( 1-1 ) * (X (NX) -X ( 1 ) ) /FLOAT (NX-1 ) 
CONTINUE 
END IF 


IF (NY. GT. 1 ) THEN 
DO 5 J = 1 , NY 

Y ( J) =Y ( 1 ) +FLOAT ( J-l ) * (Y (NY) -Y ( 1 ) > /FLOAT (NY-1 ) 

5 CONTINUE 
END IF 

NZ=NHH 

IF (NZ.GT. 1 ) THEN 
DO 6 K= 1 , NZ 

Z (K) =FLOAT (K-l ) *H/FLOAT (NZ-1 ) 

6 CONTINUE 
END IF 


NOM= 1 

DO 10 L=1 , NHH 

HH (L) =FLOAT (L-l ) *H/FLOAT (NHH-1 ) 
10 CONTINUE 

C 

TMOLT =TM-T 0 
V=VWELD*0. 0254/60. 

V2A=-V/ (2. * ALPHA) 

QS=QORIG 

Q4K= 1 . / (4. *3. 141 59* THK ) 

C 

DO 70 N= 1 , NHH 
HHH=HH(N> 

ZZ (1 ) =-16. *H-HHH 
ZZ (2) =-16. *H+HHH 
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DO 20 K=3, 34, 2 

ZZ ( K > — ZZ (K-2) +2. *H 
ZZ (K+l ) = ZZ (K-l) +2. *H 
20 CONTINUE 

DO 60 K= 1 , NX 
DO 50 L=1 , NY 
DO 40 J= 1 , NZ 
DO 30 M= 1 , 34 
ZZZ=Z ( J ) -ZZ (M) 

X X X = X (K) 

YYY=Y (L) 

R=SQRT <XXX*#2+YYY**2+ZZZ**2) 

RX=R+X X X 

C IF (RX . GT. 10. ) GO TO 60 

RX=V2A*RX 
EXPAN=EXP(RX) 

TEMPO (N, K, L, J) =TEMFO <N, K, L, J ) +Q4k#EXPAN/R 
30 CONTINUE 

40 CONTINUE 

50 CONTINUE 

60 CONTINUE 
70 CONTINUE 
QM=0 . 0 
C 

110 CONTINUE 

DO 150 K=1,NZ 
DO 140 J=1 , NY 
DO 130 1=1, NX 

TEMP ( I , J, K) =0. 0 
130 CONTINUE 

140 CONTINUE 
150 CONTINUE 

DO 190 1=1, NX 
DO 180 J=1 , NY 
DO 170 K= 1 , NZ 

IF (NOM. GT. 1 ) THEN 
DO 160 N= 1 , NOH 

TEMP ( I , J, K>=TEMP ( I , J, K) +QS/ZMAXO#DHZ (N) * TEMPO <N, I , J, K ) 
160 CONTINUE 

ELSE 

TEMP ( I , J , K ) =QS* TEMPO ( 1 , I , J , K) 

END I F 

170 CONTINUE 

180 CONTINUE 
190 CONTINUE 
ZMAX=0. 0 
M 1 =0 
M2=0 
M3=0 

DO 220 K=1 , NZ 
YMOLT <K) =0. 0 
DO 210 1=1, NX-1 

I F ( TMOLT . GT . TEMP ( 1 , 1 , K ) . AND. TMOLT. LE. TEMP ( 1 + 1 , 1 , K> ) THEN 
XM1 =X ( I)+(X (1+1) -X (I) ) * (TMOLT-TEMP(I, 1,K> ) / 
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ro m 


( TEMP (1+1,1, K ) -TEMP ( I , 1 , p> ) 
M1=M1+1 

XMQLT (Ml ) =-1000#XMl 
ZMOLT ( M 1 ) = 1 000 . tZ ( K ) 

ZMAX=AMAX1 ( ZMAX , Z <K) ) 

11 = 1+1 
END IF 


IF (TMOLT. LE. TEMP ( I , 1 , K ) . AND. TMOLT. GT. TEMP ( 1 + 1 , 1 K ) 
XM2=X ( I ) + (X ( 1+1 ) -X ( I ) ) * (TMOLT-TEMP ( I , 1 , K ) ) /’ 

( TEMP ( 1 + 1 , 1 , K) -TEMP ( 1 , 1 , K) ) 

M2=M2+ 1 


$ 

$ 


XMOLTO (M2) =-1000. *XM2 
Z MOLT 0 ( M2 ) = 1 000 . # Z ( K ) 


ZMAX=AMAX 1 (ZMAX, Z (K) ) 
12=1 
END IF 

DO 200 J=1 , NY-1 


IF ( (TMOLT. GE. TEMP < I , J , K ) . AND. TMOLT. LT. TEMP ( I , J + l 
( TMOLT . LE . TEMP ( I , J , K ) . AND . TMOLT . GT . TEMP ( I ! J + 1 
YM=Y ( J) + (Y (J+l ) -Y (J) ) / ( TEMP ( I , J+l , K) -TEMP ( I J 

(TMOLT-TEMP ( I, J,K) ) ' ? ’ 

IF ( YMQLT ( K > .LT. 1000. #YM) JY=J 
YMOLT ( K ) =AMAX 1 (YMOLT(K) , 1000. 0*YM) 

Z M Y ( K > = 1 000 . 0* Z ( K ) 


M3=K 
ENDIF 

00 CONTINUE 

10 CONTINUE 
220 CONTINUE 

C 

WRITE (*,*) ’ M 1 = ’ , M 1 , ’ 11=’, n » 

IF (Ml . LT. NZ > THEN 


12=’ , 12 


$ 


DO 240 1 = 11, 12 
DO 230 I M=2 , NZ 

IF (TMOLT. LE. TEMP (1,1, IM-1 > . AND. TMOLT. GT. TEMP ( I 1 
ZM1-Z ( IM-1 >+(Z ( IM) -Z (IM-1 >)* (TMOLT-TEMP ( I 1 IM- 
(TEMP ( 1 , 1 , IM) -TEMP (1,1, IM-1 ' ^ 

M1=M1+1 


XMOLT (Ml ) =-1000. *X ( I ) 
ZMOLT (Ml ) =1000. 0*ZM1 
ZMAX=AMAX 1 (ZM1 , ZMAX) 
GO TO 240 
ENDIF 

230 CONTINUE 

240 CONTINUE 
ENDIF 


C 

c 

c 


248 


DZ ( 1 ) =0. 5* ( HH ( 2 ) — HH ( 1 ) ) 

DZ (M2) =ZMAX-0. 5* (HH (M2) +HH (M2-1 ) ) 

^ ( * ’ * ) DZ ( 1 ) = ? , DZ ( 1 ) , ’ DZ (M2) =’ DZ (M ? ) 

DO 248 L=2, M2-1 

DZ (L) =0. 5* (HH (L) +HH (L+l ) ) -0. 5* (HH (L-l ) +HH (L) ) 
WRITE ( * , * ) DZ (L)=’ ,DZ (L) 

CONTINUE 


1 cr 
1 -J 


) THEN 


, K ) > . DR . 

!> K) ) ) T> ..M 
<) > 


IM) ) THEN 
1 ) ) / 
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QM=0 . 0 

DO 250 K= 1 , M2 

QM=QM+2. *V*DZ <K) *YMOLT (K) *HFUSI*DENS*0. 001 
250 CONTINUE 

QSQS=QQF: I G-QM 

WRITE <*,*> 7 QSQS= 7 ,QSQS, 7 QOR I G= 7 , QOR I G , QM= 7 ,QM 

DQS=1 . -QSQS/QS 
IF (ABS (DOS) . GT. 0. 05) THEN 
QS=QS* ( 1 . -0. 2*DQS> 

GO TO 110 
END I F 

IF (NON. EQ. 1 ) THEN 
N0M=N0M+1 
NQH=M2 
ZMAXO=ZNAX 

DO 252 1=1, NOH 
DHZ ( I ) =DZ ( I ) 

252 CONTINUE 
GO TO 110 
END IF 


WRITE(*,*) 7 ZMAX= 7 , ZMAX, 7 ZMAX0= 7 , ZMAXO 

WRITE (*, *> 7 ZMAX/H= 7 , ZMAX/H 

IF ( (l.O-ZMAX/H) . LT. l.E-6)THEN 

WRITE (*,#)’ A keyhole has been -formed ! ! ! 7 

ELSE 

WRITE < * , # ) 7 No keyhole ! ! ! 7 

END IF 

300 CONTINUE 
MM=M1+M2 

WRITE (8, # > 7 NN= 7 , MN 
DO 310 1=1, M2 

XMOLT ( I+Ml ) =XMOLTO (M2-I+1 ) 

ZMOLT ( I+Ml ) =ZM0LT0 (M2-I+1 ) 

310 CONTINUE 

DO 320 J=1,MM 

WRITE (B,*) XMOLT (J) , ZMOLT (J) 

320 CONTINUE 

C 

WRITE ( * , * ) 7 M3= 7 , M3, 7 JY= 7 ,JY 

IF (M3. LT. NZ) THEN 
DO 333 J=JY, 1,-1 
ZM0=0. 0 
DO 332 1=1, NX 
DO 331 K=2, NZ 

IF ( (TMOLT. LE. TEMP ( I , J , K-l ) . AND. TMOLT. GT. TEMP ( I , J , K) ) .OR. 
% (TMOLT. GE. TEMP ( I, J, K-l) . AND. TMOLT. LT. TEMP ( I , J , K) ) ) THEN 

ZM3=Z (K-l )+(Z(K)-Z (K-l) ) * ( TMOLT— TEMP ( I , J , K- 1 ) ) / 

$ (TEMP ( I , J , K) -TEMP ( I , J , K-l ) ) 

ZMO=AMAX 1 ( ZMO, ZM3) 

GO TO 332 
END IF 

331 CONTINUE 

332 CONTINUE 


Page 16 



n n n 


M3=M3+ 1 

YMOLT ( M3 ) = 1 OOO . 0 #Y<J> 

Z M Y ( M3 ) = 1 000 . 0 # Z MO 

C WRITE'#, *> 7 J= 7 , J, 7 YM= 7 , YMOLT (M3) , 7 ZM= 7 , ZMY (M3) 

333 CONT I NUE 

END IF 

DO 335 J=1 , M3 

WR I TE ( ? , * ) YMOLT < J ) , ZMY ( J ) 

335 CONTINUE 

C 

IF ( ( 1 . O-ZMAX/H) . LT. 1 . E-6) THEN 
WR I TE ( 1 0 , 3000 ) QSQS, QOR I G , QM 
ELSE 


WRITE (10, 3005) QSQS, QQRIG, QM, 100*ZMAX/H 
END IF 
RETURN 


3000 FORMAT (/3X, ’ 

QSQS 

QORIG 

QM 

Key ho 1 e ( Y / N ) ? ’ 

* 

//3X, 3F10. 

*— 1 

Y’ ) 



3005 FORMAT ( /3X , “ 

QSQS 

QORIG 

QM 

Keyhole (Y/N)? 7 

$ 

//3X , 3F10. 

•— * * 

N 7 , F6. 1 , ’ 

7.’ ) 



END 


SUBROUT I NE CON VO ( R J , H J , T J , PLM , SDM , HW , PCONV ) 

COMMON MMM 

COMMON TABTEM (20) , TAB I PL (20) ,TABCPL(20> , TABKF'L (20) ,TABISD(20) 
COMMON ITMAX, DHMAX, DTMAX 
COMMON TAMBI , DNOZL, FMF'L 

COMMON TABC3D (20) , TABKSD (20) , TABNF'L (20) , TABNSD (20) 

COMMON TW4 , MREG4 , NREG4 
N0M=0 

200 CONTINUE 
NOM=NOM+ 1 

IF(NOM.GT. ITMAX) THEN 

WRITE'#,*) 7 Erro Stop ITERATIONS EXCEED ITMAX IN C0NV4 7 
STOP 
END IF 

VJM=PLM+SDM 
write ( * , * ) 7 TJ= 7 ,TJ 
HM=.5* (HJ+HW) 

POW=VJM*HM 

CALL AVERT ( PLM , SDM , F'OW , TABTEM , TAB I PL , TAB I SD , TM ) 

CALL I NSEP (TABTEM, TABCF'L, T J , PLCF') 

CALL I NSEP ( TABTEM , TABCSD , T J , SDCF' ) 

CF'= (PLCP*PLM+SDCP#SDM> /VJM 
CALL I NSEP ( TABTEM , TABCF'L , Tm , PLCF'M ) 

CALL INSEP (TABTEM, TABCSD, Tm, SDCF'M) 

CF'M= <PLCPM#PLM+SDCPM#SDM) / VJM 
CALL INSEP (TABTEM, TABKF'L, TJ , F'LAMD) 

CALL INSEP (TABTEM, TABKSD, TJ, SDAMD) 

AMD= (F’LAMD*PLM+SDAMD*SDM) /VJM 
CALL I NSEP ( TABTEM , TABNF'L , T J , PLN ) 

CALL INSEP (TABTEM, TABNSD, TJ, SDN) 
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VISC= (PLN*PLM+SDN*SDM> /VJM 
REN=2 . # VJM/RJ /VI SC/3 . 14159 
PRA=CP*VISC/AMD 
VNU=0. 78*REN**0. 5*PRA**0. 33 
ALFAJ = . 5* VNU* AMD/RJ 

PC0NV=3. 14159*RJ*RJ*ALFAJ* (HJ-Hw) /CPM 

RETURN 

END 
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*** 

INPUT DATA FOR ALUMINUM *** 


SYMBOL 

UNIT 

DESCRIPTION 

DATA 

MMn 


Number ot point in table 

IB 

TABTEM 

K 

T emperaetur e table 

Table attached 

TABIF'L 

kJ/kg 

Enthalpy of plasma gas 

Table attached 

TABCF'L 

k J / k g . 1 

K Speci-fic heat o-f plasma gas 

Table attached 

TABKF'L 

kW/M.K 

Thermal conductivity of plasma gas 

Table attached 

TABISD 

kJ/kg 

Enthalpy o-f shield gas 

Table attached 

TABCSD 

kJ/kg. I 

< Specific heat of shield gas 

Table attached 

TABKSD 

kW/M.K 

Thermal conductivity of shield gas 

Table attached 

TABMF'L 

MHO/M 

Electric conductivity of plasma gas 

Tabl e attached 

TABMSD 

MHO/M 

Electric conductivity of shield gas 

Table attached 

TABNPL 

kg/M. S 

Viscosity o-f plasma gas 

Table attached 

TABNSD 

kg/M. S 

Viscosity of shield gas 

Table attached 

NREG1 


Segment number between electrode 
and orifice 

1 

NREG2 


Segment number of orifice 

1 

NREG3 


Segment number of standoff 

1 

NREG4 


Segment number of workpiece 

7 

VA 

V 

Anode potential drop 

. 30000E+01 

VC 

V 

Cathode potential drop 

. 45000E+01 

FIA1 

V 

Electrode work function 

. 30000E+01 

FIA2 

V 

Orifice work function 

. 45000E+0 1 

FI A3 

V 

Workpiece work function 

. 38000E+01 

TIME 

mS 

Time duration for straight polarity 

Vari able 

TIMER 

mS 

Time duration for reverse polarity 

Var i able 

DISL1 

M 

Distance betwteen electrode and 
orifice 

Vari able 

DISL2 

M 

Axial length of orifice 

Variable 


1 


DISL3 

hi 

D 1 s t a r c; e a f s t a n d a f f 

Var i ab 1 e 

DI3L4 

M 

Thickness of workpiece 

. 6 3500 E -02 

WE 

Caul omb 

Electron charge 

. 1 &020E— 1 3 

VVK 

erg/ K 

Baltcman constant 

. 1 3800E- 1 5 

ARIL 

A 

Pilot electric current in straight 
polarity 

Var i ab 1 e 

APILR 

A 

Pilot electric current in reverse 
pol ar i ty 

Var i ab 1 e 

AMAI 

A 

Main electric current in straight 
pol ar i ty 

Var i abl e 

AMAIR 

A 

Main electric current in reverse 
pol ar i ty 

Var i able 

DNOZL 

hi 

Orifice diameter 

■ olQOOE — 02 

DSHID 

M 

External diameter o-f shield 

. 18000E-01 

VPL 

scf h 

Plasma gas -flow rate 

Vari abl e 

VSD 

sc-f h 

Shield gas flow rate 

Var i able 

F'LMOL 

kg/kmol 

Molecular weight of plasma gas 

. 39950E+02 

SDMOL 

kg/kmol 

Molecular weight of shield gas 

. 40030E+0 1 

UR 

V 

Reverse polarity voltage rise 

. 1 14S0E+02 

ITMAX 


Maximum number of iteration 

100 

DHMAX 


Maximum relative error of enthalpy 

. 50000E— 02 

DTMAX 

K 

Maximum error of temperature 

. 50 0 00 E +00 

TW1 

K 

Electrode wall temperature 

. 25000E+04 

TW2 

K 

Orifice wall temperature 

. 40000E+03 

TW4 

K 

Workpiece wall temperature 

. 91600E+03 

HSUR2 

kJ/kg 

Plasma enthalpy correponding to TW2 

. 20800E+03 

HSUR4 

kJ/kg 

Jet enthalpy correponding to TW4 

. 47600E+03 

DXW 

M 

Electrode length 

Vari able 

TAMBI 

K 

En vi r omental temperature 

. 29800E+03 

DPOLAR 

M 

Electrode di ameter 

. 39600E-02 

AMDW 

kW/M.K 

Thermal conductivity of electrode 

. 23500E— 01 


z 


PS I A 


Fraction power .Loss due to anode 
and cathode drop 


PAMEs I 

N/M2 

En vi r onmen t al pressure 

. 10132E+06 

CF’SD 

kJ/kg.K 

Specific heat of shield gas with 
environmental temperature 

. 51930E+01 

ALFA 

1'12/S 

Thermal diffusivity o-f workpiece 

. 52410E-04 

VWELD 

i pm 

Welding speed 

- 9 1 000E+0 1 

NX 


Number o-f computing nodal in X 
coor d i na t e 

15 

NY 


Number o-f computing nodal in Y 
coordi nate 

7 

X (1) 

M 

Beginning position of computed 
nodal in X coordinate 

70000E— 02 

X (NX) 

M 

Ending position of computed 
nodal in X coordinate 

. 70000E-02 

Y ( 1 ) 

M 

Beginning position of computed 
nodal in Y coordinate 

. 00000E+00 

Y (NY) 

M 

Ending position of computed 
nodal in Y coordinate 

. 70000E— 02 

TM 

K 

Melting temperature of workpiece 

. 91600E+03 

TO 

K 

Initial temperature of workpiece 

. 29S00E+03 

THK 

kW/M.K 

Thermal conductivity of workpiece 

. 12000E+00 

EPSL 


Emissivity of workpiece 

. 1 0000E+00 

DENS 

kg/ M3 

Density of workpiece 

. 26500E+04 

HFUSI 

kJ/kg 

Melting latent heat of workpiece 

. 3S700E+03 
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Attached Table for Aluminum 


T emper aeture 

(K) 




. 3000E+04 
. 1 00CE+05 
. 1500E+05 
. 2000E+05 


. 6 000 E +04 
. 1 100E+05 
. 1600E+05 
. 2200E+05 

. 7000E+04 
. 1 200E+05 
. 1700E+05 
. 240GE+05 

. 80 00E +04 
- 1 300E+05 
. 1800E+05 

. 9000E+04 
. 1 400E+05 
- 1 900E+05 

Enthalpy 

a-f PI asma Gas 

(kJ/kg) (ARGON) 



. 1562E+04 
- 6125E+04 
. 3553E+05 
. 5945E+05 


- 3138E+04 
. 8164E+04 
. 4384E+05 
. 6862E+U5 

* C-652E+04 
. 1 183E+05 
. 4957E+05 
. 8667E+05 

. 4229E+04 
. 1 778E+05 
. 5372E+05 

. 4966E+04 
. 2613E+05 
. 5661E+05 

Sped f i c 

Heat 

of F'l asma 

Gas ( k J / k g . K ) ( ARGON ) 


. 5 1 92E+00 
. 1511E+01 
. 9253E+01 
. 3136E+01 


. 5192E+00 
. 2721E+01 
. 7159E+01 
■ 6573E+01 

. 5401E+00 
. 4677E+01 
. 4731E+0J 
. 1 122E+02 

. 6280E+00 
. 7243E+01 
.321 1E+01 

. 8918E+00 
. 9253E+0 1 
. 2734E+01 

Thermal Conduct i vi t v o-f 

Plasma Gas (kW/M.K) 

(ARGON) 


. 1 440E-03 
- 6250E-03 
. 24 1 7E-02 
. 2644E-02 


• 1 660E— 03 
. 9610E-03 
. 2315E-02 
. 3094E-02 

. 2030E-03 
- 1403E— 02 
. 2252E-02 
. 3542E-02 

. 2720E-03 
. 190 IE— 02 
. 2294E-02 

. 4020E-03 
. 2297E-02 
. 2436E-02 

Enthal py 

of Shield Gas (kJ/kg) (HELIUM) 



- 2597E+05 
-5194E+05 
- S773E+05 
. 2624E+06 


• 3 1 1 6E+05 
- 5728E+05 
. 1C27E+06 
. 4322E+06 

- o<S35E+u5 
. 6295E+05 
. 1245E+06 
.61 9SE+06 

. 4154E+05 
. 6935E+05 
. 1562E+06 

. 4674E+05 
. 771SE+05 
. 201 1E+06 

Specific 1 

Heat 

o-f Sh i e 1 d i 

Gas ( kJ/kg. K) (HELIUM) 


- 5193E+01 
. 5266E+01 
. 1232E+02 
. 6880E+02 


.5193E+01 
- 5453E+01 
- 1 777E+02 
. 9526E+02 

. 5193E+01 
. 5924E+01 
. 2588E+02 
. 8384E+02 

. 5193E+01 
. 6952E+01 
. 3720E+02 

. 5229E+0 1 
. 8913E+01 
. 5183E+02 

Thermal Conductivity of Shield Gas (kW/M.K) 

(HELIUM) 


. 9100E-03 
. 1300E-02 
. 6900E-02 
- 3700E-0 1 


. 9300E-03 
. 1620E-02 
. 1174E-01 
• 3860E— 0 1 

. 9500E-03 
. 1940E-02 
. 1658E-01 
. 3420E-0 1 

. 1028E-02 
. 3060E-02 
. 2260E-01 

. 1 164E-02 
. 4980E-02 
. 2980E-0 1 

Electric Conductivity of 

Plasma Gas (MHO/M) 

(ARGON) 


. 1030E+02 
. 2730E+04 
. 7490E+04 
. 1 040E+05 


. 1010E+03 
. 3730E+04 
. B200E+04 
- 1 080E+05 

. 36 1 0E+03 
. 4740E+04 
. 88 1 0E+04 
. 1050E+05 

• 9230E+03 
. 5740E+04 
. 9730E+04 

. 1 770E+04 
. 6670E+04 
. 9S80E+04 
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El ectr 1 c 

. 4285E-01 
. 3793E+02 
- 2002E+04 
. 5720E+04 

Vi scosi ty 

. 1 290E-03 
. 2900E-03 
- 71 10E-04 
. 2300E-04 

Viscosi ty 

- 1000E-03 
. 1840E-03 
.21 8OE-03 
- 1950E-03 


Conductivity o+ Shield Gas (MHO/M) (HELIUM) 


. 7852E+00 
. 1000E+03 
. 2887E+04 
. 7093E+04 

ot Plasma Gas 


• 2070E+01 
- 2512E+03 
. 3648E+04 
* 8794E+04 

(kg/M.S) (ARGON) 


. 5458E+01 
. 6310E+03 
. 436SE+04 


. 2090E-03 
. 2840E-03 
. 4460E— 04 
. 2100E-04 


. 2330E-03 
. 2450E-03 
. 3150E-04 
- 1560E-04 


. 2560E-03 
. 1810E-03 
. 2570E-04 


o-f Shield Gas (kg/M.S) (HELIUM) 


. 1 400E-03 
. 1920E-03 
. 2200E-03 
. 1 600E-03 


. 1550E-03 
- 2000E-03 
. 2180E-03 
. 1 000E-03 


. 1620E-03 
. 2060E— 03 
.21 10E-03 


. 1439E+02 
. 1261E+04 
. 5137E+04 


. 2770E-03 
. 1 1 70E-03 
. 2380E-04 


. 1700E-03 
. 2120E-03 
. 2050E-03 
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Appendix Input Data in Program (I) (Mild Steel) 


SYMBOL 

DESCRIPTION 

MMM 

Number of point in table 

TABTEM 

Temperature table 

TAB I PL 

Enthalpy of plasma gas (Argon) 

TABCPL 

Specific heat of plasma gas 
(Argon) 

TABKPL 

Thermal conductivity of plasma 
gas (Argon) 

TABISD 

Enthalpy of shield gas (Argon) 

TABCSD 

Specific heat of sheild gas 
(Argon) 

TABKSD 

Thermal conductivity of sheild 
gas (Argon) 

TABMPL 

Electric conductivity of plasma 
gas (Argon) 

TABMSD 

Electric conductivity of sheild 
gas (Argon) 

NREG1 

Segment number between 
electrode and orifice 

NREG2 

Segment number of orifice 

NREG3 

Segment number of standoff 

NREG4 

Segment number of workpiece 

VA 

Anode potential drop 

VC 

Cathode potential drop 

FIA1 

Electrode work function 

FIA2 

Orifice work function 

FIA3 

Workpiece work function 

TIME 

Time duration for straight 
polarity 

TIMER 

Time duration for reverse 
polarity 

DISL1 

Distance between electrode and 
orifice 


DATA 


26 


Table Attached 


Table Attached 


Table Attached 


Table Attached 


Table Attached 


Table Attached 


Table Attached 


Table Attached 


Table Attached 



3V 


4.5V 


2.6V 


4.5V 


4.2 V 


19ms 



1 . 1938mx<t 























































Input Data in Program (II) 



AMAIR 


DNOZL 



PLMOL 


SDMOL 



ITMAX 


DHMAX 



DESCRIPTION 


Vertical length of orifice 


Distance of standoff 


Thichness of workpiece 


Electron charge 


Boltzman constant 


Pilot electric current in 
straight polarity 


Pilot electric current in 
reverse polarity 


Main electric current in 
straight polarity 


Main electric current in 
reverse polarity 


Orifice diameter 


External diameter of sheild 
gas 


Plasma gas flow rate 


Shield gas flow rate 


Molecular weight of plasma 
gas (Argon) 


Molecular weight of shield 
gas (Argon) 


Constant in computation of 
reverse polarity rise 


Constant in computation of 
reverse polarity rise 


Maximum number of iteration 


Maximum relative error of 
enthalpy 


Maximum error of temperature 


Electrode wall temperature 


Orifice wall temperature 


Workpiece wall temperature 


3.175mm 


Variable 


6 . 35mm 


1 ♦ 602xl0~ 19 coulomb 


1.38X10’ 16 erg/k 


Variable 


Variable 


Variable 


Variable 



Variable 


39.5 kg/kmol 


39.5 kg/kmol 



0.005 


2500K 


1727K 


a 



















D : Input Data in Program (III] 



DESCRIPTION 

Plasma enthalpy 
corresponding to TW2 (Argon) 

Jet enthalpy corresponding 
to TW4 (Argon) 


Electrode length 


Environmental temperature 


Electrode diameter 

Thermal conductivity of 
electrode 


Fraction of power loss due 
to anode and cathode drops 


Environmental press ure 

Specific heat of shield gas 
with environmental 
temperature (Argon) 

Number of computing nodal in 
X coordinate 

Number of computing nodal in 
X coordinate 


Number of computing nodal in 
Z coordinate 

Number of heat source 

Position of computed nodal 
in X coordinate 


Position of computed nodal 
in Y coordinate 


Position of computed nodal 
in Z coordinate 


Position of heat source 


Workpiece thickness 


Welding speed 


Thermal diffusivity of 
workpiece 


DATA 

190 kJ/kg 


900 kJ/kg 


68 . 072mm 


298K 


3 . 9624mm 


23.5 W/m.K 



101320 N/m 2 


0.5205 kJ/kg. K 



Variable 


Variable 


5 .45x10** m 2 /s 
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D : Input Data in Program (IV) 


SYMBOL 

DESCRIPTION 

DATA 

THK 

Thermal conductivity of 
workpiece 

21.63 W/m.K 

DENS 

Density of workplace 

7900 kg/m 3 

HFUSI 

Melting latent heat 

272 kJ/kg 

TMOLT 

Temperature rise 

1429K 

















Appendix D : Attached Table for Mild Steel (I) 


51 , 99L 
51: 


52: THERMAL CONDUCTIVITY 0 

53 : 

F PLASMA GA 

S ' ( kW/m . K ) 

c ARGoti ; 

54: .1013E-03 

55 : 

. 1057S-03 

. 1174E-03 

. 1432E-03 

. 1S33E-03 

56 

57 

. 2320E-03 

. 4279E-03 

. 6550E-03 

. 9916E-03 

. 1463E-02 

58 

59 

. 2050E-02 

. 2563E-02 

. 2S75E-02 

. 2480E-02 

. 2348E-02 

60 

61 

.2374E-02 

. 2507E-02 

-2697H-02 

. 3148E-02 

. 3632E-02 

62 

63 

. 4126S-02 

. 462 3 E -02 

. 5120E-02 

. 5618E-02 

. 6116E-02 

64 

65 

66 
67 

. 63S5E-02 

• 

• 


- 

68 

69 

ENTHALPY OF 

SHIELD GAS 

( kJ /kg ) 



70 

71 

. 15S2E+04 

.2082E+04 

. 2603E+04 

. 3125E+04 

. 3654E+04 

72 

73 

-4215E+04 

-4894E+04 

. 5931E+04 

. 7892E+04 

M173E+05 

74 

75 

. 1823E-i-05 

. 268CE+05 

. 3549E+05 

. 42S7E+05 

. 4874E+05 

76 

77 

. 5351E+05 

. 57S2E+-05 

.6143E+05 

. 638SE+05 

. 6590E+C5 

78 

79 

. 1112E+06 

. 1298S+06 

- 1441E+06 

. 1609E+06 

. 1831E+06 

80 

81 

82 

83 

- 1958E+06 





84 

85 

SPECIFIC HEAT OF SHIELD 

GAS (kJ/kg 

.:<) cargoes) 

86 

87 

. 5205E+00 

. 5205E+00 

. 5209E+00 

. 5238E+00 

. 5382ji+00 

88 

89 

. 5972E+C0 

. 737 1E+0G 

137QE+01 • 

• • 727 iSEf-01 

" . 51I5E+01 

90 

91 

-7791E-01 

-b37o2. - r\.)i 

. 3 163E+0 1 

. 65ol£L+0i 

. 523SS+01 

92 

93 

.4370E+C1 

. 3901E+01 

.377 1E+01 

-5252S+01 

...1139E+02 

94 

95 

. . 1I73E+02 

. 7515Et-0i 

. 7 329E+0 1 

.9700E+01 

. 1229E+-02 

96 

. 1255E-02 






97: 

98: 

99: 
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D- 

Attached 

Table for Mild Steel (II) 

- 

c i cc r 






3 1 1 






52: 
53 : 

THERMAL CON 

’vCTIV'TY 0 

F PLASMA GA 

S ( kW./tn . K ) C AZ&o/J ~) 

54 : 
55 : 

. iO 13E-03 

. 10572-03 

.11742-03 

. 14322-03 

. 1 9 oo w Go 

56: 

57: 

.23202-03 

.42792-03 

. 65502-03 

.99162-03 

• 1463E-02 

56: 

59: 

.20502-02 

.25632-02 

.26752-02 

.24802-02 

.23482-02 

60: 

61: 

. 237 4E-02 

.25072-02 

.26972-02 

.31482-02 

. 3632E-02 

62: 

63: 

. 41262-02 

. 43232-02 

.51202-02 

. 56132-02 

.61162-02 

64: 

. 63652-02 





65 : 






66: 






67: 






65: 

ENTHALPY OF 

SHIELD GAS 

/ ?. 7 \ 

v ) 

( AZ&ori ) 


69: 




70: 

71: 

. 15522+04 

. 20822+04 

. 26032+04 

.31252+04 

. 36542+04 

72 : 
73: 

.42152+04 

.43942+04 

.59312+04 

. 78922+04 

.11732+05 

. 74: 
75: 

. 18232-05 

. 26302-05 

• oo49i— Oo 

. 42872+05 

. 48742+05 

76: 

t i l 

. 52512-05 

• - > o w + 0 r. 

.61432+05 

. 6y3oi+0c 

. 85902+05 

73: 
73 : 

. 11122-06 

. 12982+06 

. 1441E+0S 

.1609E-06 

.18312+06 

80: 

. 19562+06 





81: 






82: 






83: 






84: 
85 : 

SPECIFIC HEAT OF SHIELD 

GAS (kJ/kg 

-K) CAZGW) 


86: 
87 : 

.52052+00 

. o20w^4-00 

. 5209E+00 

. 5236E+00 

.53322+00 

83: 

89: 

.53722+00 

. 797 1E+QQ 

. 1370E+01 

•2719E-01 

.51152+01. 

90: 

91: 

-779l2+0i 

. 8973E+0 1 

. 8 163E+0 i 

. 658 iE^O i 

. 52362+0 1 

92: 

S3: 

. 4.j70c.+0 1 

390IE+0I 

. 377 12+01 

. 5252E+01 

.11332+02 

94: 

95: 

.11732+02 

7515E+-01 

. 73292+01 

.9700Et-01 

.12232-02 

96: 

. 12652-02 






37 : 
96: 

cq . 
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D: Attached Table for Mild Steel (III) 


100 . 147L 






IOC : 
10 1 : 

- n n ;\M \ ... 

activity of 

•“•r: x i . jAc 

' kW/li . V 


102 : 
103: 

. 10 I3E-03 

. 1057E-03 

- 1174E-03 

. 1432E-03 

. 1933E-03 

104: 

105: 

o c o n. zr o -*5 

- — *-* w ^ — — ».♦ 

. 4273E-03 

. 6550E-03 

-991SE-03 

. 1463E-02 

106: 

107: 

. 2050E-02 

. 2563E-02 

. 2675E-02 

. 2480E-02 

. 2348E-02 

108: 
" 109: 

.2374E-02 

. 2507E-02 

.2S97E-02 

. 3 143E-02 

. 3632E-02 

110: 

111: 

. 4I26E-02 

- 4623E-02 

. 5 120E-02 

. 5618E-02 

.6116E-02 

112: 
113: 
114: 
115 : 

. 63S5E-02 





. 116: 
- 117 : 

ELECTRIC CONDUCTIVITY OF 

PLASMA GAS 

(MK0/m) CAAGO/V ) 

118: 

119: 

. 5 9 60S— 04 

. 1270E+00 

. 1030E+02 

. 1010E^03 

. . 36 10 tv' 03 

120: 

121: 

.923OS+03 

. 1770E+04 

. 2730E-04 - 

. 3730E+04 

.47405-I-04 

po . 

123: 

. 5740 E -04 

. 367CE+04 

. 7 490 S-*-04 

. 8200an-04 

.86l0E-r04 

124: 

1 o ^ . 

-L _ v . 

. 9 i 3 0 E 0 4 

. 9880E-1-04 

. 1040E+05 

.1080E+05 

. 10 5 0.7;- 05 

126: 

127: 

. 1020S+05 

. 1 Z.*r0 0 

. 1090 £+05 

. 1 130 E r-0 5 

. 1160E+05 

170. 

^ m 

. 11 70 £-05 





129: 
130 : 
131: 






1 OO . 
iuc . 

133 : 

ELECTRIC CC 

N DUCT I VI TY OF 

SHIELD GAS 

(MK0/mj (/ 

134: 
135 : 

. 5960E -04 

. 1270Ei-0c' 

1030E+02 

• 10 iv a+0o 

. 36 10 E, 03 

136: 

137: 

- 9^30e.f03 

. 1770E+04 

2730E-1-04 

. 37 30 E t-04 

. 4740E+04 

138: 

139: 

. 5740E+04 

.S670E4-04 

7430E-r04 

.820CE+-04 !! 

. S810E+04 

140: 
141 : 

-9730E+04 

-9330E-t-04 

1040E+05 

- 1080E4-05 

. 1050K ;0 5 

142: 

143: 

. 1020E-05 

• I040i-*-05 . 

I0S0S+05 

. 1130E+05 

.11607,-05 

144: 

. 1170E-1-05 
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< E * Nomenclature 


A c “ Crown reinforcement area (m 2 ) 

A e = Cross-sectional area of electrode (m 2 ) 
A x = Orifice circumferential area (m 2 ) 

A r ~ Root reinforcement area (m 2 ) 

C p = Constant pressure specific heat (J/kg*k) 
D = Orifice diameter (m) 
d e = Electrode diameter (m) 

D s / « = External diameter of shielding gas (m) 
d c = Crown bead width (m) 
d r = Root bead width (m) 
dv 

— = Voltage drop per unit length (V/m) 


— Feed wire diameter (m) 
g = Gravitational acceleration (m/s 2 ) 
h c = Crown height (m) 

h exit = Plasma jet enthalpy at the exit of workpiece (J/kg) 

h * = P (j/kg^ enthal Py corresponding to the orifice wall temperature 

^plasma = Plasma enthalpy (J/kg) 
h r = Root height (m) 


h * ' 1 ' h i , o “ Plasma enthalpy at entrance and exit ends, rest : - v 

°f gap between electrode and orifice (J/kg) 

h ‘''' ***'•■ P of“in« al (J/kg) entranCa andexit ends ' respectively, 
h3 ' ‘ ' hj ' * " P or:L“ P c y oi a uVn n wk C g ) and eXit endS ' --P^ively, 


1 


n «,i 


V' 


k = 




L e = 


L x = 


L 2 = 
L 3 = 


l 4 = 


m p j 

m sd 

M m 

M,d 
n = 
P = 

p n : 

P E = 


P H = 


P «d 
P 


J/ 1 


J / 2 


J/3 


M/ 1 


' h * ' 0 P ^ sma enthalpy at entrance and exit ends, respect - 
of workpiece keyhole (J/kg) p 

Im = !? a f n ?J? ctric current in straight and reverse 
polarities, respectively (A) 

Ip " = nibn> eleCtric current in straight and reverse 
polarities, respectively (A) 

Thermal conductivity of plasma jet (W/m*k) 

Thermal conductivity of electrode (W/m*k) 

Length of electrode (m) 

Distance between electrode and orifice (m) 

Axial length of orifice (m) 

Column distance of standoff (m) 

Thickness of workpiece (m) 

= Plasma gas mass flow rate (kg/s) 

= Shielding gas mass flow rate (kg/s) 

= Molecular weight of plasma gas (kg/kmol) 

= Molecular weight of shielding gas (kg/kmol) 

Parabolic exponent of reinforcement shape 
External pressure (N/m 2 ) 

= Partial pressure of plasma gas (N/m 2 ) 

Power input at electrode (W) 

Environment pressure (N/m 2 ) 

= Partial pressure of shielding gas (N/m 2 ) 

= orific^twf 119 P ° Wer inpUt in thS gap between electrode 

Joule heating power input in the orifice (W) 

= Joule heating power input in the standoff column (W) 

— work function induced power input (W) 


,vcly, 


and 


2 



», 2 


V, 1 


V, 2 


' 1 , 1 ' *1 


Orifice anode drop induced power input (W) 

= Workpiece work function induced power input (W) 

- workpiece anode and cathode drops induced power input. (w) 
P 


2 , 1 ' 


3 , 1 / 


4,1' 


' ° P of S aan a hPt P u° Wer a ? entrance and exit ends respect: j >>ely 
of gap between electrode and orifice (W) Y 

• " P o! S o?mce P °( W W) r “ e " trance and exit respectively, 


P o? S sLndofTo e o\^ e w anCe a " d SXit 6ndS reSy ‘ 

4 ' ° wo?kpiecTkeyhole r (W^ ^ 6Xlt SndS/ respec * 


'.vc.ly, 


i s 

Q c 

Qx 

Qx , c 

Q»,e 

®X , 0 


Q s , r 
Q « , c 

Qs 

Qh,c 

Q»,. 

Qw 

R 

r e 

R o 
R i 


Of 


(W) 


°y 

t .("> 


Plasma arc discharged power loss (at exit end of ke 
= Power loss in the gap between electrode and orifi< 

Total power loss in the orifice (W) 

= Convective heat transfer power loss in the orifice (w> 

= Free electron power loss in the orifice (W) 

= °tSS? ^ ower 1 . oss from the surface of the orifice contr -r 

anod^ d^ Wer work function and part of power V 

anode drop at orifice (W) p tx x 

= Standoff column radiative power loss (W) 

Standoff columln convective power loss (W) 

~ Standoff column total power loss (W) 

= Convective heat transfer workpiece power loss (W) 

Free electron power loss in workpiece (W) 

— Total power loss in workpiece (W) 

Universal gas constant (J/mol‘K) 

= Effective radius of plasma jet (m) 

= Orifice radius (m) 

= Distant between source and point of computation at workpiece 
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respectively" (s° r Straight and reverse polaritic 
T E - Effective temperature of plasma jet; Electrode temper ^ re 
T« - Environment temperature ( K ) C 

T» = Melting temperature of workpiece (K ) 

temperature; Initial temperature of workpiece (K) 

U " Weldi - — -U- ^t velocity (m/s) " 

U “ = -“^nt flow velocity (m/s) 

U ° = Jet core velocity (m/s) 

V A/ v b/ V c , V„, V £ , V 


F ' V ° * straight polarity electric 
shown in Figure 1, at point 
' D ' E ' F / and G, respect -> 


v »'/ V,', V c ' t V,', Vp', V r ' Tr / jj 

c P° larit y elect 

VOltclCTS . shnwn -i ** in * 


je, 

(V) 


voitage shown in Figure I, 
at points A', B ', C', t>' E' 
V , P1 F / and G', respective < v) 

g ^-Lasma volumetric x. 

electrode and ortfSr^vi, in the gap between 

V. = Plasma volume in the standoff column (m’/s, 
v w Wire feeding speed (m/s) 

«• - initial gap distance between two workpieces (m) 

X " “ MStanCe ° f Jet initial Region (m, 
y — Coordinate axis 

Yc = Coordinate to 110 «o * 

computation (m) f ° r Spatial temperature distribution 

2 = Coordinate axis 

computation^)^ US6d f ° r spatial temperature distribution 

z . = Heat source s P a tial point of distribution (m) 

Angle (rad, ; workpiece thermal diffusivity („»/., 

OF rvi ^ 


4 


6 = Angle (rad) 


aV *' iVc A ^° ) de and cathode electric potential drops, respectively 
"■*' "■' ^ Strai9ht “* 


aV/ ' * V " polarities^respectively^vf *" Stral9ht *" d ™ 

( = Emissivity of workpiece 

c R = Reverse polarity voltage rise (V) 

<r s — Surface tension coefficient (N/m) 

* - Fraction of heat source distributed in spatial coordinate 
= Electrode work function (V) 

= Orifice work function (V) 

= Workpiece work function (V) 

* ~ e”de 1 and°fa?ho!e C drops Workplece ' P ower lo3s contributed by 
tj = Dimensionless coordinate 

P P1 = Mass of plasma gas in unit volume of plasma jet 

P ad - Mass of shielding gas in unit volume of plasma jet (kg/m 3 ) 


5 



Basic Input Parameters (I) 
(Aluminum #2219-T87) 


Symbol 

Description 

Data Caput 

*■1 

Distance between electrode and nozzle 

1.1176 mm 


Vertical length of nozzle 

3.27 6 0 mm 

^4 

Thickness of workpiece 

6.35 mm 

l e 

Length of electrode 

69.6 am 

D 

Nozzle diameter 

3.17'. m 

d e 

Electrode diameter 

3.9624 mm 

k e 

Thermal conductivity of electrode 

23. G :W/)tt • 0 K 

m 

Argon gas mass flow rate 

Varicwai a 

♦ e 

Electrode work function 

3.0 V 

4>k 

Nozzle work function 

4.5 V 


Workpiece work function 

3.8 V 

*V a 

Anode potential drop 

3.0 V 

AV C 

Cathode potential drop 

4.5 V 

V 

Main electric current in straight polarity 

Variable 


electric current in reverse polarity 

Variable 

V 

Pilot electric current in straight polarity 

Variable 

V 

Pilot electric current in reverse polarity 

Variable 

t + 

Time duration for straight polarity 

19 ms 

t_ 

Time duration for reverse polarity 

4 ms 

t e 

Electrode temperature 

2500 °K 

T * 

Nozzle temperature 

400 °K 

T 0 

Environment temperature; Initial 

298 °K ! 


temperature of workpiece 
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(II) 


Symbol 


♦ 


Basic Input Parmeters 

(Aluminum # 2219 -T 87 ) 
Description 

wall tempera tur e C0rreSP0n< ^^ n9 to °rifice 

temperature alPy corres P°nding to keyhole 
cathode n drops° Wer l0SS due to ano <*e and 

“r*£“e in L h r t ; our « 

thickness P fraction of plate 

Welding traveling speed 


Thermal 


conductivity of workpiece 


a 


Pw 

Qf 

€ 

T_ 

m 


Thermal diffusivity of workpiece 

Density of workpiece 

heat of workpiece 
Emissivity of workpiece 
Melting temperature of workpiece 


Data Input 
kJ 

0.20U - — 

7 

kJ 

0.476 — 

'3 

0.0 

°*5 - i.o 

Vari %bXe 

kv 

0.13 

•«i • °ic 

5.670 lo-5 

2650 kg/m 3 

387 k v j,., 

0.1 

916 °K 


^GE IS 

of poor eoaury 
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Appendix F Block Diagram of Penetration Depth Calculation 

Structure Diagram of WELDO Computer Program 
(Computation of penetration depth) 



Explanation of WELDO Program Structure 
1. Main Program 


Input physical parameter tables of plasma and shield gas and initial parameters of <*(■; 

notenri^ mbUQ ° n of related Physical parameters of plasma and shield gas , such as e !- ‘ 
p nnal , power , power loss , enthalpy, etc., between torch and workpiece. 

2. Subroutine AVERT 


Compute mean temperature of plasma and shield gas mixture with given 
plasma and shield gas flow rate and power of the mixture. 


3. Subroutine INSEP 


Compute mterpalanvely temperature or enthalpy and their corresponding values of relate, 
physical parameters of plasma and shield gas with input values of plasma and shield ga 
material constants and given temperature or enthalpy of plasma and shield gas. 

4. Subroutine CONVO 

Compute power transmission from plasma jet to workpiece through heat process. 

5. Subroutine DEPTH 


Compute temperature distribution , geometry of molten pool of workpiece and 
depth of partial penetration of welding process. 


*.ss . 











Flow Chan of WELDO Computation Program 



Open data documents : 
WELDTAB.DAT 
WELD IN. USE 
WELDOU.USE 


i 




Input physical parameters of plasma and shield gas 

X==3 1 

Input initial parameters of welding process 

1 L_ 

Print input datas to WELDOU.USE 

[ J 

| Initial mass flow rale and velocity plasma and shield gas 

1 j 

Poxomcrer mstnoution of plasma at location between 
electrode and orific* 

r , 

Related parameters of plasma in the orifice 

X r 1 

Parameter distribution within stanoff column 

■ — 1 

v-au DEPTH to compute penetration depth 
and print output 


i 


I 


t 


I 

1 


I 










Appendix G Block Diagram of VPPA Modeling Structure 


Structure Diagram of Diagram of WELDBOTH 
(VPPA Modeling) 



Explanation of WELDBOTH Program Structure 

1. Main Program 

Input physical parameter tables of plasma and shield gas and initial parameters of welding process. 
Compute distribution of related physical parameters of plasma and shield gas .such as electric 
potential , power , power loss , enthalpy , etc., between torch and workpiece . 

2. Subroutine AVERT 

Compute mean temperature of plasma and shield gas mixture with given 
plasma and shield gas flow rate and power of the mixture . 

3. Subroutine INSEP 

Compute interpalatively temperature or enthalpy and their corresponding values of related 
physical parameters of plasma and shield gas with input values of plasma and shield gas 
material constants and given temperature or enthalpy of plasma and shield gas. 

4. Subroutine CONV4 

Compute power transmission from plasma jet to the wall of keyhole in workpiece through 
heat convection process. 

5. Subroutine DEPTH 

Compute crown and root widths based on temperature distribution of workpiece and geometry 
of keyhole. 












Flow Chan of WELDBOTH Computation Program 


C Stan ") 

I -I I 

Open data documents: WELDI^SE^ 

WELDO U. USE 

| — - i _ 

Input physical parameters of plasma and shield gas 

— I 

Input initial parameters of welding process 

— ' ‘ i 1 — 

Initial mass flow rate and velocity plasma and shield gas 

Parameter distribution of plasma at location between 
electrode and orifice 

— ♦ 

Related parameters of plasma in the orifice 

— " I 

Parameter distribution within stanoff column 


3 

Call CONV4 to compute distribution of convective 
heat loss in kevhote 

■ f 

Compute power loss distribution in keyhole 

y__ 

Call WIDTH to compute geometry of keyhole 








Appendix H 

Comparison Between Results of Model 
Computations and Experimental Measurements 


Figures H-l to H-63 show the comparison between the results of 
model computations and experimental measurements. „ odel 
computations are based on the mathematical formulations illustrated 

“ thlS r6POrt ' WhUe measurements are adopted from 

^^experiments accomplished b y Hichols Hesearch Corporation 

in Figures H-l to H-31, comparisons are made between the model 
computed values of forward electric voltage and reversed electric 
voltage versus that of the values of experimental measurements. 

e results show that the tendency and magnitudes of the model 

computations are verv mm-h -t-K 

veiry much the same as thaf a.. 

as that of the experimental 

measurements . 

In Figures H-32 to H-38, comparisons are made for leading edge 
ang e between model computated values and that of the values of 
experimental measurements. There are so many parameters which 
could be changed in the experimental measurements. It is very hard 

O Choose a series of experimental measurement with one parameter 

as a variable and with the rest of fho 

of the parameters unchanged. 

ompansons made in Figures H-32 to H-38 follow the order of the 
sequence of experimental measurement of leading edge angle between 
»odel computation and experiment. The results also show that the 
tendency and magnitudes of the model computations are very much the 
same as that of the experimental measurements. 


/ 



Figures H-39 to H-62 show the comparison for crown width and 
root width between model computed values and that of the values of 
experimental measurements. The results show that the model 
computation and the experimental measurement are relatively in good 
manners of agreement. It is noted that some parameters provide a 
great influence on the model computation of the crown and root 
widths. However, within the range of proper control of parameters, 
it can anticipate to predict a fairly reasonably good value", of 
model computation in crown and root widths agreeable in compe d 
with the experimental measurements. 

As to the problems of abnormal welding, model formiilr ! l r. 
illustrated in Section VII of this report have been adopts i to 
calculate the probability of the keyhole full and partdi. 
penetrations of the welding workpiece. Table H-A shows the results 
of model computations for keyhole full and partial penetration of 
the workpiece in comparison to experimental measurements of I* 
Numbers 7A to 7P listed in the report accomplished by Niohol? 
Research Corporation (1992) , while Table H-B illustrates the 
results of the similar model computations in comparison to 
experimental measurements of ID Numbers 8A to 8N listed in the s.-.amt 
report. It shows that the results of model computation are 

basically all agreeable with the experimental measurements based on 
the points of the possibility of keyhole full and pao! . 
penetrations of the welding workpiece. 

For the purpose of making easier comparison for the 
possibility of keyhole full and partial penetrations of welding 
workpiece among so many different controlled parameters of welding 
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process between bode! computation and experimental me 

Set ° f dimensionless parameters are — urements, a 

Figure H -63 with Proposed and illustrated in 

VI^k-A, + ep L) , as “ " ^ h0ri2 ° ntal CO ° rdlnate 

m a PJ-if) ] as the vertical , . 

‘-ical coordinate. Here 
v _ ... e ' 
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k 

At 
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p 

Ls 


ten 


Welding Speed 
Thickness of workpiece 
Thermal diffusivity of workpiece 
Thermal conductivity of workpiece 
Temperature difference between molt 
initial temperatures of workpiece 
Density of workpiece 
Latent heat of molten workpiece 

0 - Total power absorption of worhpiece during th, 

welding process. 

Recall that there is 

“ ere is no available experimental measures * 

values because of the diff icu i t v , suremen 3 

values experimentally q acc “hate measurements of Q 

y* Q values adopted -in r>* 

calculated based on the physical model i lustr „ ^ "" 

3 ) Of this report. ln Ration (7- 

Figure H-63 shows the critic 
Partial penetrations of welding JIT- * eyh ° le 
computation. i„ this figure, the i ef ^ 

tone of keyhole fun DB “PPer region shows the 

1 e rul1 penetration while 

illustrates the zone * . r *<Jht-hand lower region 

° partial penetration. This 

computed based on the Th Curve is 

on the condition of the he^ 

workpiece to att-a • ° m P° r *-i° n of the 

attain " olt en temperature of the met., 
locations of each evne • tal * Precise 

each experimental measurements 

urements with experimental id 

fs ^ 

^ POOR Qu AUT y 


.nd 


numbers 7 and 8 from the Nichols report are also plotted. In these 
experiments, aluminum workpiece with thickness of 1/4 •• ( 6 .35 mm) 
are used. The results show that the model computations predict the 
keyhole full and partial penetrations of welding completely 
agreeable with that of the experimental measurements. 
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"Table. H-A (.a) 


IDi 

CALCULATION CCUFA?£3 WITH 

Ea.-EEINENT 

e i Y/N > ? 

Exp. 

i 

> 

* 1 
r 1 

: (kwo 

j 

i 

QORIC 

O.W) 

(kW) 

KtvhVL 

Cal. 

7A1 

i 

i 

i i r>c ‘ c > 
1 
1 

1 001 

.143 

Y 

Y 

7A2 

! 1 . 306 

J 

1.474 

.163 

Y 

Y 

7 A3 

; 1.4C4 

l 

1.604 

.199 

Y 

Y 

7B1 

! 1.241 

1 
1 

1.390 

.143 

Y 

Y 

7B2 

! 1.337 

1 

J 

1.515 

.173 

Y 

Y 

7S3 

\ 1.444 

i 

i 

1.657 

.212 

Y 

Y 

7C1 

! 1.300 

i 

f 

1 . 467 

.167 

Y 

Y 

7C2 

1.394 

1 

1.590 

.196 

Y 

Y 

7C3 

1 

1 

1.507 

1.740 

.233 

Y 

Y 

7D1 ! 
( 
1 

1 OOn 
i - w/ / / 

1.567 

.191 

Y 

Y 

732 ! 

f 

1 /»C? 

1.703 

oor 

Y 

Y 

7D3 ; 

1 

l 

1 . 594 

1.853 

1 

• 4-W *T 

Y 

Y 

7E1 1 
1 
1 

1 . 223 

1.340 

.115 

N 51.4 % 

N 

7E2 j 

1 

1 

1.347 

1 . 479 

1 ■so 

- 

N 55.6 % 

N 

7 S3 i 

1 

1 

1.456 

1.614 

.153 

N 63.0 % 

N 

7F1 f 

1 

1 

1.257 

1.373 

.121 

N 51.9 % 

N 

7F2 ! 

1 

1 

1.335 

1.526 

.141 

N 59.3 % 

N 

7F3 ! 

1 

1 

1.502 

1.673 

.171 

N 64.3 % 

N 

7G1 

1 

1 

1.321 

1.455 

.134 

N 53.6 % 

N 

7G2 ! 

1 

r 

1.432 

1.535 

.135 

N 61.7 % 

N 

7G3 ! 
• 
1 

1.563 

1.755 

.133 

N 66.2 % 

N 


“Table H - A Cb ) 


7H1 

! 1 35^ 

1 K * i 





\ 

t 


. 145 

N 60. 

0 % N 

7H2 

! 1.516 

1.631 

.174 

N 64. 

3 % N 

7H3 

I 1.665 

1.874 

.209 

N 75. 

0 % N 

711 

J 0.381 

1 .087 

.106 

N 77. 

14% INTERMIT 

712 

1.098 

1.232 

.134 

Y 

INTERMIT 

713 

1 1.231 

1.385 

1 7? 


INTERMIT 




Y 

7J1 

! 1.004 

1.111 

.107 

N 

INTERMIT 

7J2 

j 1 . 133 

1.284 

.151 

Y 

INTERMIT 

7J3 

! 1.246 

1.426 

.180 

Y 

INTERMIT 

7K1 

1.074 

1.202 

.128 

N 

INTERMIT 

7K2 

1.137 

1 

1.364 

.167 

Y 

INTERMIT 

7K3 

1.305 

1.501 

.196 

Y 

INTERMIT 

7T *} l 

f iJ-L i 

1 

f 

1.137 

1.290 

.153 

Y 

INTERMIT 

7f o i 

i i -it- i 

r 

1.261 

1.444 

.133 

V 

1 

INTERMIT 

7L3 ! 

1.380 

1.592 

.212 

Y 

INTERMIT 

7M1 j 

0.961 

1.071 

.110 

N 46.7 % 

N 

7M2 : 

1.111 

1.248 

.137 

N 50.5% 

N 

7M3 ; 

1.228 

1.334 

.156 

N 52.7% 

N 

7N1 J 

1.011 

1.123 

.113 

N 47.0% 

N 

7N2 • 

1.155 

1.299 

.144 

H 51.2% 

N 

7N3 ! 

1.280 

1.444 

.164 

N 53.2% 

N 

701 ! 

1 

r 

1.072 

1.202 

.130 

N 49.4% 

N 

702 ! 

1 

1 

1.216 

1.370 

.154 

N 52.1% 

N 

703 1 

_ 1 
1 

1.345 

1.513 

.152 

N 54 . 1% 

N 


'Table H-A tc) 


7Pi : 

1.183 

1.311 

1 4 O 

N 

C i oo/ 

Ox -j 

7P2 ! 

1.311 

1.473 

.139 

H 

53.6" 

7P3 ! 

1.445 

1 . 623 

.134 

N 

CO 1 o/ 

'O''? . x*> 














ttateral 

OORIG = Power obsorbed by workpiece. 


1/4” ??1 9 ~ TS7 Aluainum C7A-7H) 

1/4 (S.vjffim) A3o iliict Steel (7I-7P) 


Gii - Latent power of workpiece. 

QSQS _ Conduction loss in workpiece ( = QORIG-QM) 
These values are calculated by the model 



Table H-B 

CALCULATION COMPARED WITH WPtr-rvr 
Estimating whether keyhole can 


-r . . • 

~‘-i ».TOt 


ID# ! 


QSGS 
! (kW) 


CALCULATION 

QOc\^a ON Keyhole? 
(kW) (kW) (Y/N) 


EXPERIMENT results 


.'Weld width! Weld height 
.'Crown Root .'Crown Root 


Bead Ap- 
pearance 


' i — ; 1 

t 

8A I 1.103 1.134 
1 
1 

.089 

N 

35.7 

_ l 

J 

I 

1 

% i 9.86 * 

1 

« 

j 

i 

l 

1.24 * 

1 

t 

l 

P 

' PP.N 

9B ! 1.101 1.180 

1 
» 

.070 

N 

85.2 

% ! 9.45 * 
1 

1.30 * 

PP,N 

ec ! 1.1C5 1.135 
1 
1 

.030 

N 

85.9 

% ! 9.40 * 
| 

1.31 * 

PP, N 

8D ! 1.103 1.183 
1 
1 

.080 

N 

85.6 

% ! 9.41 * 

1 

1.08 * 

PP, N 

8E ! 1.103 1.183 
1 
J 

.030 

N 

85.6 

% .' 9.58 * 

I 

1.52 * 

PP.N 

8F i 0.935 0.982 

l 

1 

.046 

N 

61.6 

% ! 5.74 * 

1.70 * 

PP, N 

8G 1 1.164 1.362 

1 
1 

.198 

Y 


! 10.03 * 

1 

1.14 0.75 

PP.PT.N 

8H I 0.979 1.084 

1 
l 

.105 

N 

76.1% 

1 

! 9.00 * 

t 

1.25 * 

PP,N 

81 0.976 1.081 

1 
1 

.105 

N 

75 . 8% 

l 8.38 * 
1 

1.85 * 

PP, D 

8J ! 0.977 1.082 
1 
1 

.105 

N 

75.9% 

1 

! 9.30 * 
1 

.1.17 * 

PP,D 

8K ! 0.977 1.082 

t 

1 

.105 

N 

75.9% 

I 9.32 * 

1.02 * 

PP, D 

8L ! 0.972 1.078 

r 

i 

.104 

N 

76.5% 

I 3.73 * 
1 

1.30 * 

P?,N 

8M ! 0.889 0.973 

1 
f 

.084 

N 

67.4% 

» 

! 7.87 * 

1 

1.11 * 

PP, N 

8N ! 1.059 1.182 
1 
1 

.123 

N 

87.0% 

! 8.23 * 
1 

* * 

PP, D 


Materal . 1/4 (6.35mm) 2219-T87 Aluminum (8A-8G) 

I/'* (6.35mm) A36 Mild Steel (8H-8N) 

QORIG = Power ohsorbed by workpiece. 

CM = Latent power of workpiece. 

&=QS _ Conduction loss in workpiece (= QORIG-QM) . 

These values are calculated by the model 
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